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ABSTRACT 


This  report  details  the  development  and  provides  the 
documentation  for  custom  computer  software  applications  that 
evaluate  antenna  parameters.  The  applications  are  written  in 
Mathcad  3.1  for  the  following  antenna  types:  linear,  planar, 
and  circular  arrays?  folded  dipole;  caged  dipole;  parabolic 
reflectors  with  helical  and  spiral  feeds.  Inputs  to  the 
applications  are  limited  to  the  antenna ' s  physical  dimensions. 
In  some  cases,  ground  parameters  are  required. 

The  chapters  are  structured  to  provide  the  user  with  the 
necessary  information  needed  to  use  and  interpret  the  software 
for  each  antenna  type.  The  software  applications  are  provided 
as  appendices  and  give  examples  of  each  antenna  type. 

Outputs  of  the  applications  provide  various  numerical  and 
performance  predictions,  as  well  as  far-field  radiation 
patterns.  The  results  computed  are  consistent  with 
predictions  provided  in  applicable  publications,  as  well  as 


those  calculated  by  numerical  antenna  analysis  programs  such 


as  ELNEC. 
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I.  BACKGROUND  AND  PURPOSE 


This  thesis  details  the  development  of  user  friendly 
mathematical  applications  capable  of  computing  the  radiation 
pattern  and  other  pertinent  antenna  parameters  of  an  antenna 
or  antenna  system  based  on  available  information,  primarily 
dimensional  information  obtained  from  photographs.  The 
applications  developed  are  compatible  with  any  IBM  personnel 
computer  using  MS-DOS  version  3.2  or  higher  and  a  math 
coprocessor.  The  applications  are  written  in  Mathcad  version 
3.1  engineering  software  and  can  be  run  using  either  Mathcad 
3.1  or  Mathcad  version  4.0.  These  applications  are  intended 
for  use  by  the  Naval  Maritime  Intelligence  Center 
(NAVMARINTCEN) . 

As  already  mentioned,  user  inputs  are  limited  to  antenna 
dimensions  based  on  physical  measurements  primarily  obtained 
from  photographs,  although  other  source  data  can  be  used.  The 
applications  are  flexible  to  the  extent  that  other  parameters 
may  be  estimated  and  used  as  inputs  to  increase  the  accuracy 
of  the  computations  to  gain  better  insight  into  the  antenna's 
performance.  The  Mathcad  applications  provide  various 
performance  predictions  as  well  as  a  graphical  representation 
of  the  antenna's  far-field  radiation  pattern.  The  necessary 
background  information  needed  to  interpret  the  application's 
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formulas  and  displays  are  provided  in  the  corresponding  thesis 
chapters . 

Dietrich  [Ref.  1]  and  Gerry  [Ref.  2]  have  completed  the 
first  and  second  reports  of  this  project,  respectively.  This 
thesis  is  the  last  in  a  series  of  three  reports  intended  to 
fulfill  the  NAVMARINTCEN  statement  of  work. 
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II.  INTRODUCTION 


The  increasingly  sophisticated  design  by  foreign 
countries  of  antennas  warrants  the  development  of  improved 
analytical  methods  to  provide  timely  and  accurate  technical 
assessment  of  these  antennas.  Through  the  process  of  reverse 
engineering,  foreign  communication,  navigation. 
Identification,  Friend  or  Foe  (IFF) ,  and  radar  antennas  are 
evaluated  to  determine  their  capabilities,  limitations,  and 
vulnerabilities.  Having  limited  information  available 
concerning  the  parameters  of  an  antenna,  usually  restricted  to 
photographs  that  provide  only  the  type  and  physical  dimensions 
of  the  antenna,  current  methods  of  antenna  evaluation  are  slow 
and  tedious. 

With  the  availability  of  high-speed  personal  computers 
matched  with  sophisticated  off-the-shelf  engineering  software, 
the  assessment  oi  antennas  using  only  photographic  information 
may  now  be  obtained  rapidly.  The  objective  of  this  report  is 
to  provide  NAVMARINTCEN  with  user  friendly  Mathcad  software  to 
achieve  this  type  of  performance  in  their  evaluation  of 
antenna  systems. 

As  with  any  computer  analysis  of  antennas,  the  accuracy 
of  the  results  depend  upon  the  sophistication  of  software 
used,  the  complexity  of  calculations,  the  reliability  of  input 
data  obtained  from  photographs,  and  the  accuracy  of  any 
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required  input  data  that  must  be  estimated.  This  report  and 
the  Mathcad  applications  are  written  with  existing  engineering 
equations  for  the  type  of  antenna  analyzed  and  inputs  are 
limited  to  physical  dimensions.  Additional  information  such 
as  ground  effects  and  characteristic  impedance  may  be  included 
to  increase  the  accuracy  of  the  analysis.  To  this  end,  this 
report  should  provide  an  initial  understanding  as  to  the 
capabilities  of  the  antenna  system  in  question. 

The  chapters  of  this  report  document  the  equations  and 
assumptions  used  for  each  specific  antenna  type  and  is  written 
as  a  comprehensive  reference  for  the  software.  Copies  of  each 
application  are  included  as  appendices  to  provide  the  user 
with  a  printed  illustration  of  the  software. 
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III.  THE  UNIFORM  LINEAR  ARRAY  AMTEWMA 

Arrays  are  highly  directional  antenna  system  with  narrow, 
steerable  beams  and  low  side  lobes  used  for  long  distance 
communication  and  radar  systems.  The  linear  array  is  one 
geometrical  configuration  that  consists  of  radiating  elements 
lying  along  a  straight  line.  As  a  result  of  the  dependance  of 
inter-element  spacing  (d)  to  the  wavelength  (1) ,  linear  arrays 
are  inherently  narrowband  antenna  systems. 

In  this  Mathcad  application  uniformly  spaced,  equally 
excited,  isotropic  point  sources  aligned  on  the  z-axis  are 
used  as  the  radiating  elements.  In  using  isotropic  point 
sources,  polarization  of  the  array  is  not  calculated  and  an 
assumption  is  made  that  the  antenna  input  resistance  (Rjn)  is 
equal  to  the  radiation  resistance  (Rr)  .  This  assumption 
results  in  calculations  for  gain  (G)  that  are  idealized 
values.  The  following  inputs  are  needed  to  implement  the 
software: 

N  =  number  of  isotropic  radiating  elements 

d  =  inter-element  spacing  between  adjacent  elements 

0„  =  direction  of  main  lobe 

f  =  frequency  of  interest 

I  =  antenna  feed  current 
0 

ZQ  =  characteristic  feed  impedance 
The  first  two  inputs  are  obtained  from  visual  data  such  as 
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photographs,  and  the  *  indicates  input  parameters  that  are 
either  known  or  are  estimated.  The  linear  array  geometry  is 
illustrated  in  Figure  3.1. 


The  total  electric  field  of  a  linear  array  in  the  far- 
field,  neglecting  mutual  coupling  between  adjacent  elements, 
is  determined  by  the  vector  addition  of  the  fields  of 
individual  radiating  elements.  It  can  be  shown  that  the  total 
field  is  equivalent  to  the  product  of  the  field  of  a  single 
element,  selected  at  a  reference  point,  and  a  factor  which  is 
referred  to  as  the  array  factor  (AF) : 

E(  total)  =  [E  (element)  ]  x  [AF]  (3.1) 

In  (3.1),  E  is  the  vector  electric  field  intensity,  while  AF 
is  a  scalar  quantity.  This  result  is  referred  to  as  pattern 
multiplication  for  arrays  of  identical  elements.  The  array 
factor  is  a  function  of  the  geometry  of  the  array  and  the 
element  excitation  amplitude  and  phase.  For  uniform 
excitation,  by  varying  the  inter-element  spacing  (d)  and/or 
the  phase  (/3)  between  the  elements,  the  characteristics  of  the 
array  factor,  and  consequently  the  total  electric  field,  of 
the  array  can  be  controlled.  [Ref.  3:  pp.  204-207] 

This  Mathcad  application  is  written  for  an  N-element 
linear  *ray  with  equally  spaced,  uniformly  excited  isotropic 
radiating  elements.  In  addition,  element  phase  is  assumed  to 
vary  linearly  along  the  length  of  the  array.  The  radiation 
pattern  of  a  specific  element  type  is  neglected  in  (3.1)  since 
in  normal  usage  it  will  have  little  effect  when  the  array 
consists  of  a  large  number  of  elements.  However,  the  user  of 
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this  software  should  be  able  to  use  information  regarding  the 
type  of  element  and  the  orientation  of  the  elements  to  narrow 
the  range  of  possible  0o's.  The  normalized  array  factor  (AFn) 
is  [Ref.  3:  pp.  212-214]: 


af„  ( e )  =  — 
«  N 


Sin|-|i|f  ( 0 ) | 
sin|-|4,(6)j 


( dimensionless )  (3.2) 


where 


(0)  =  kdcos(Q)  +  p  ( radians ) 


(3.3) 


P  =  -kdcosid  )  ( radians ) 


(3.4) 


In  (3.4),  0o  is  the  direction  of  the  maximum  array  factor 
value  (i.e.,  for  a  broadside  array  0O  =  n/2,  for  an  end-fire 
array  with  the  main  beam  directed  at  180*,  0o  *  n ,  and  for  a 
phased  array  with  the  main  beam  directed  to  ±60',  0O  =  tt/3)  , 
P  is  the  phase  excitation  difference  between  the  elements,  k 
is  the  wavenumber  (2jt/A)  ,  and  d  is  the  inter-element  spacing. 

With  the  linear  array  aligned  along  the  z-axis,  the  array 
factor  is  a  function  of  0  but  not  of  <p.  The  visible  region  of 
the  array  in  terms  of  0,  i|r ( 0 )  ,  and  the  inter-element  spacing 
(d/A)  is: 

0  <  0  <  7i  ( radians )  (3.5) 
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P~Jcd  <  i|»  (0)  <  P+Ad  ( radians ) 


(3.6) 


d /A  determines  how  much  of  the  array  factor  appears  in  the 
visible  region  as  defined  by  (3.5)  and  (3.6).  Although  d/A 
provides  guidance  in  determining  the  frequency  to  select  for 
this  Mathcad  application,  it  should  be  noted  that  most  arrays 
are  designed  with  the  inter-element  spacing  less  than  one 
wavelength,  usually  close  to  a  half -wavelength.  Thus,  linear 
arrays  are  characteristic  narrowband  antenna  systems.  [Ref. 
4:  p.  123] 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

i  z  1.6  A  ( meters )  (3.7) 

r  *  5  (Nd)  ( meters )  (3.8) 

r  ;>  2  ( meters )  (3.9) 

The  array  length,  Nd,  in  (3.8)  and  (3.9)  is  the  maximum 
dimension  of  the  antenna  (the  array  length  is  assumed  to 
include  a  distance  of  d/2  beyond  each  end  element  [Ref.  5: 
p . 55 ] ) .  The  minimum  distance  to  the  far-field  is  found  by 
selecting  the  largest  value  of  (3.7),  (3.8),  and  (3.9). 

Radiation  intensity  (U(0))  is  related  to  the  power 
radiated  in  a  given  direction  and  is  independent  of  the 
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distance  to  the  observation  point  [Ref.  3:  p.  27;  Ref.  4;  p. 
33].  For  linear  arrays  U(8)  is  given  by  [Ref.  3:  pp.  229- 
233] ; 

£/{0)  =  (AFn) 2  {w/  solid  angle)  (3.10) 

Radiated  power  (Prad)  is  the  total  power  radiated  by  the 
antenna  and  is  obtained  by  integrating  U(0)  over  a  surface 
surrounding  the  antenna  [Ref.  4:  p.  33].  Since  the  geometry 
of  the  linear  array  is  aligned  on  the  2-axis  and  the  array 
factor  is  independent  of  0,  Prad  is  thus  [Ref.  3:  pp.  229-233]: 

P[ad  =  2%  f*U(6)  sin(6)  dB  (W)  (3.11) 

JO 

Directivity  (Dc)  is  the  maximum  value  of  directive  gain 
where  directive  gain  is  the  ratio  of  U(0)  in  a  specific 
direction  to  the  average  radiation  intensity  [Ref.  4:  pp.  34- 
36].  Hence,  Do  is  given  by  (Ref.  3:  pp.  229-233]: 

D0  =  ( dimensionless )  (3.12) 

Prad 

In  (3.12),  Umax  equals  unity  and  occurs  at  0O  for  the  uniform 
linear  array. 

Effective  isotropic  radiated  power  (EIRP)  is  a  commonly 
used  communications  term  that  is  defined  as  the  product  of 
antenna  gain  and  total  power  accepted  by  the  antenna  from  the 
transmitter.  EIRP  is  determined  as  [Ref.  4:  p.  62]: 
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EIRP  =  PradDa 


( W ) 


(3.13) 


Assuming  that  the  power  delivered  to  the  array  consisting 
of  N-isotropic  point  sources  is  equivalent  to  the  radiated 
power,  we  can  estimate  the  radiation  resistance  (Rr)  of  the 
antenna  as  [Ref.  3:  p.  55]: 

Rr  =  2  (Q)  (3.14) 

Using  the  same  logic,  we  obtain  the  input  resistance  (Rin)  of 
the  array  made  up  of  isotropic  point  sources  as  equivalent  to 
the  radiation  resistance: 

Rin  =  Rr  (Q>  (3.15) 

The  exact  gain  (G)  of  the  array  cannot  be  calculated  with 
the  data  assumed,  but  an  idealized  value  can  be  computed  based 
on  the  assumption  made  to  determine  Rjn  and  the  assumptions 
made  to  calculate  total  antenna  efficiency  (et) .  The  total 
antenna  efficiency  is  given  by  [Ref.  3:  pp.  44-45]: 

tc  =  erecd  =  l  -  Jr |2  ( dimensionless )  (3.16) 


where 


ecd  =  1  ( dimensionless ) 


tz  =  1  -  | r  |2  (dimensionless) 


(3.17) 


(3.18) 
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( dimens i onl ess ) 


(3.19) 


p  _  Rin  ~  Zo 
Rin  +  20 

In  (3.17),  ecd  is  the  conduction  and  dielectric  efficiency 
which  is  unity  for  an  ideal  lossless  antenna.  In  (3.19),  Z0 
is  the  characteristic  impedance  of  the  transmission  line  and 
r  is  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(3.18)  and  the  ideal  gain  for  a  linear  array  is  (Ref.  3:  pp. 
43-44] : 

G  =  ttDa  ( dimensionless )  (3.20) 

G(dB)  =  10log10  (tcD0)  ( dB )  (3.21) 

As  a  result  of  constructing  the  linear  array  with 
isotropic  point  sources,  the  polarization  loss  factor  (PLF) 
and  antenna  polarization  are  not  calculated.  Therefore,  the 
maximum  effective  aperture  (Aw)  of  the  array  when  used  as  a 
receiver  is  calculated  assuming  negligible  polarization 
mismatches.  The  is  then  (Ref.  3:  p  63]: 

Aem  =  —  PLF  t™2)  (3.22) 

e  4  71 


where 


PLF  =  l  ( dimensionless ) 


(3.23) 
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The  assumptions  used  in  finding  Rr  and  Aw  result  in  an 
idealized  value  in  calculating  the  effective  height  (h^)  .  The 
effective  height  is  [Ref.  6:  p.  42]:] 


*ea,  =  2, 


( meters ) 


(3.24) 


where  r?0  is  the  characteristic  impedance  of  free  space  (120ir)  . 

Two  displays  are  produced  in  this  Mathcad  application. 
The  first  display  is  the  rectangular  representation  of  the 
magnitude  of  the  array  factor  from  -it  to  n .  The  second 
display  is  the  polar  representation  of  the  magnitude  of  the 
array  factor. 

The  Hansen-Woodyard  end-fire  array  was  not  addressed  in 
this  application  nor  were  arrays  of  uniform  spacing  but 
nonuniform  amplitude  (i.e.,  Dolph-Tschebyschef f  array).  In 
addition,  nonuni formly  spaced,  uniformly  excited  arrays  were 
not  addressed  by  this  application. 

The  results  obtained  with  this  application  are  identical 
to  those  found  in  [Ref.  3:  pp.  216-240]  with  respect  to  array 
factor  patterns  and  directivity.  Table  3.1  is  a  comparison  of 
data  from  [Ref.  3]  for  a  10-element  equally  excited,  uniformly 
spaced,  end-fire  linear  array  while  Figure  3.2  is  a  comparison 
of  the  array  factor  pattern  of  the  same  array.  Of  note,  the 
array  factor  patterns  from  this  application  and  the  patterns 
from  [Ref.  3]  are  identical  as  expected. 
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Table  3.1  End-fire  Array  Data  Comparison 


ANTENNA 

REFERENCE 

CALCULATED 

DIMENSIONS 

DIRECTIVITY  ( DQ) 

DIRECTIVITY  ( DQ) 

N  =  10 

d  =  1/4  m 

10.0  dB 

10.003  dB 

0=0° 

0 

=  -kd 

14 


1 


FIGURE  3.2  Array  factor  pattern  of  10-element  equally 
excited,  uniformly  spaced  end-fire  linear  array  with  d  = 
1/4  m,  X0  =  O',  and  /3  =  -kd 
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Table  3.2  is  a  comparison  of  data  from  [Ref.  3]  for  a  10- 
element  egually  excited,  uniformly  spaced,  broadside  linear 
array,  and  Figure  3.3  is  a  comparison  of  the  array  factor 
pattern  of  the  same  array.  The  array  factor  patterns  from 
this  application  and  the  patterns  from  [Ref.  3]  are  identical. 


Table  3.2  Broadside  Array  Data  Comparison 


1 — - — — — - - 

ANTENNA 

REFERENCE 

CALCULATED 

DIMENSIONS 

DIRECTIVITY  (D  ) 

DIRECTIVITY  (Dc) 

N  =  10 

d  =  1/4  m 

6.99  dB 

7.14  dB 

6o  =  tt/2 

o 

ll 

1 

16 


FIGURE  3.3  Array  factor  pattern  of  10-element  equally 
excited,  uniformly  spaced  broadside  linear  array  with  d  = 
1/4  m,  0O  =  jt/2,  and  0  =  0 
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IV.  THE  P LAMAR  ARRAY  ANTENNA 


Planar  arrays,  frequently  used  long  distance 
communications  and  radar  systems,  exhibit  characteristics 
analogous  to  those  of  linear  arrays  but  have  additional 
capability  to  control  and  shape  the  radiation  pattern.  Planar 
arrays  are  more  versatile  than  linear  arrays,  providing 
symmetrical  patterns  with  lower  side  lobes,  and  can  be  used  to 
scan  the  main  beam  of  the  antenna  array  toward  any  point  in 
space  [Ref.  3:  p.  261].  As  with  linear  arrays,  planar  arrays 
are  narrowband  antenna  systems  due  to  the  dependance  of  the 
inter-element  spacing  (dx  and  dy)  on  the  wavelength  (1) .  The 
following  inputs  are  required  analyze  a  planar  array: 

M  =  number  of  isotropic  radiating  elements  in  the 
x-direction 

N  =  number  of  isotropic  radiating  elements  in  the 
y-directicn 

dx  =  inter-element  spacing  between  adjacent  elements  in 
the  x-direction 

dy  =  inter-element  spacing  between  adjacent  elements  in 
the  y-direction 
f  =  frequency  of  interest 
0„  =  direction  of  main  lobe  at  0  =  0  * 

<t>0  =  direction  of  main  lobe  at  <t>  -  <t>0  * 

I0  =  antenna  feed  current 
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Z0  =  characteristic  feed  impedance 
The  first  four  inputs  are  data  obtained  primarily  from  visual 
data  such  as  photographs,  and  the  *  indicates  input  parameters 
that  are  either  known  or  estimated.  The  planar  array  geometry 
is  illustrated  in  Figure  4.1. 


Z 


FIGURE  4.1  Geometry  of  M  x  N  element  planar  array 
of  isotropic  point  sources 
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This  Mathcad  application  is  written  for  an  M  x  N  planar 
array  with  uniformly  spaced,  equally  excited,  identical 
radiating  elements  positioned  on  the  x-y  plane.  This 
application  may  also  be  used  to  approximate  the  array  factor 
and  field  patterns  for  arrays  that  have  elliptical  or  circular 
geometries  by  assuming  an  ellipse  is  equivalent  to  a  rectangle 
or  a  circle  is  equivalent  to  a  square,  respectively  [Ref.  6: 
p . 187 ] .  Since  isotropic  point  sources  are  assumed  as  array 
elements,  the  polarization  of  the  array  is  not  calculated,  and 
an  assumption  is  made  that  the  antenna  input  resistance  (Rjn) 
is  equivalent  to  the  radiation  resistance  (Rr) .  This 
assumption  results  in  values  for  gain  (G)  that  are  idealized. 
As  with  the  linear  array,  the  radiation  pattern  of  a  specific 
element  type  is  neglected  in  the  analysis  of  the  planar  array 
since  in  normal  usage  it  will  have  little  effect  for  a  planar 
array  with  a  large  number  of  elements.  However,  as  with  the 
linear  array,  the  user  of  this  software  should  be  able  to  use 
information  regarding  the  type  of  element  and  the  orientation 
of  the  elements  to  narrow  the  range  of  the  possible  0o,s  and 

Vs- 

The  normalized  array  factor  (AFn)  is  calculated  for  far- 
field  observations  neglecting  mutual  coupling  between  adjacent 
elements  with  equal  amplitude  excitation  and  is  [Ref.  3:  pp. 
260-263] : 
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where 


AFn(0,4>)  = 

sin(-j j 

sin| 

(4.1) 

M  sin|  ) 

N  si 

(  4fy(0,<t>)  ] 

\  2  JJ 

ijrx  (  0,  <J) )  = 

Jcdxsin  ( 0  )  cos  ( <J> ) 

-  P* 

( radians ) 

(4.2) 

4ry(0,4>)  = 

Jcdysin  <  0  )  sin  ( <J) ) 

+  Py 

(zadians) 

(4.3) 

Px  =  -Jed  sin(  0O)  cos  (<t>  )  ( radians ) 


(4.4) 


P  =  -Jed  sin  (0O)  sin  (<j>0)  ( radians ) 


(4.5) 


Equation  (4.1)  is  the  normalized  array  factor  that  is  a 
function  of  both  0  and  <t>  where  (4.2)  and  (4.3)  are  the  array 
factor  phase  shifts  in  the  x  and  y  directions,  respectively. 
(4.4)  and  (4.5)  are  the  progressive  phase  shift  between 
adjacent  elements  in  the  x  and  y  directions,  respectively, 
with  the  main  beam  directed  along  0  =  0O  and  <t>  =  <po.  In  (4.2) 
through  (4.5),  k  is  the  wavenumber  (2 n/X) . 

When  the  inter-element  spacing  between  the  elements  (dx 
and  dy)  equals  or  is  greater  than  the  wavelength  (A.)  ,  multiple 
maxima  of  equal  magnitude  are  formed.  To  avoid  creating 
multiple  maxima,  the  inter-element  spacing  must  be  less  than 
a  wavelength  at  the  frequency  of  operation  [Ref.  3:  p.  263]. 
With  the  dependance  of  inter-element  spacing  on  the 
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wavelength,  planar  arrays  are  characteristic  narrowband 
antenna  systems. 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r  2  1.6  1  (meters)  (4.6) 


r  z  5 y/  ( Mdx ) 2  +  ( Ndy)2  ( meters ) 


2{(Mdx)2  +  ( Ndy )2  ) 
r  2:  — i ^ ^ ( meters ) 


(4.7) 

(4.8) 


The  length  of  an  array  is  assumed  to  include  a  distance  d/2 
beyond  each  end  element  [Ref.  5:  p.  55].  Therefore,  the 
maximum  dimension  of  the  planar  array  is: 

^  ( Mdx ) 2  +  ( Ndy ) 2  (meters)  (4.9) 

The  minimum  distance  to  the  far-field  is  found  selecting  the 
largest  value  of  (4.6),  (4.7),  and  (4.8). 

Directivity  (D0)  for  planar  arrays  is  determined 
approximately  by  [Ref.  7:  p.  78]: 


D0  = 


_ _ 

/o7;/Vn^)sin(0)  d8d<|> 


(4.10) 


The  maximum  value  of  the  numerator  in  (4.10)  is  unity  at  0  = 
0O  and  <p  =  <p0  since  the  array  factor  is  normalized.  The 
denominator  is  the  radiated  power  (Prad)  of  the  antenna  where 
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Pr«d  obtained  fay  integrating  over  all  angles  around  the 
antenna  above  the  x-y  plane  (0  <  0  <  it/ 2)  [Ref.  3:  p.  28;  Ref. 
7;  p.  78].  This  assumes  the  array  is  always  pointed  toward 
one  side  or  the  other  of  the  x-y  plane. 

Effective  isotropic  radiated  power  (EIRP)  for  planar 
arrays  is  defined  as  the  product  of  antenna  gain  and  the  total 
power  radiated.  EIRP  is  determined  as  [Ref.  4:  p.  62]: 

EIRP  =  PradD0  iff)  (4.11) 

Assuming  that  the  power  delivered  to  planar  arrays 
constructed  of  M  x  N  isotropic  radiating  elements  is 
equivalent  to  the  radiated  power,  we  estimate  the  radiation 
resistance  (Rr)  of  the  array  as  [Ref.  3:  p.  55]: 

Rr  =  2 \P^~-  (Q>  (4.12) 

Uol2 

Since  the  radiating  elements  are  assumed  to  be  ideal  point 
sources,  the  input  resistance  of  a  planar  array  (Rin)  is 
assumed  to  be  equivalent  to  the  radiation  resistance: 

Rin  =  Rz  (Q)  (4.13) 

The  gain  (G)  of  a  planar  array  cannot  be  calculated  given 
the  information  assumed  to  be  available,  but  an  idealized 
value  can  be  derived.  Another  factor  affecting  the  gain  is 
the  assumption  made  when  determining  the  total  antenna 
efficiency  (et) .  The  total  antenna  efficiency  for  an  antenna 
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is  [Ref.  3:  pp.  44-45]: 


e£  =  erecd  =  1  -  j T |2  ( dimensionless )  (4.14) 


where 


ecd  =  1  ( dimensionless ) 


(4.15) 


er  =  l  -  I r I2  ( dimensionless )  (4.16) 

r  =  — —  (dimensionless)  (4.17) 

°in  + 

In  (4.14),  ccd  is  the  conduction  and  dielectric  efficiency, 

which  is  unity  for  an  ideal  lossless  antenna.  In  (4.17),  Z0 

is  the  characteristic  impedance  of  the  transmission  line  and 
r  is  the  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(4.15)  when  ideal  lossless  antenna  is  assumed.  The  ideal  gain 
of  a  planar  array  in  dimensionless  and  decibel  quantities  is 
then  [Ref.  3:  pp.  43-44]: 

G  =  ttDa  ( dimensionless )  (4.18) 

G(dB)  =  10  log  10(ttDo)  (dB)  (4.19) 

The  maximum  effective  aperture  (A^)  for  a  planar  array 
is  derived  by  assuming  polarization  mismatches  are  negligible 
(this  application  is  constructed  using  isotropic  radiating 
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elements  for  which  polarization  loss  factor  (PLF)  and  antenna 
polarization  are  not  calculated)  [Ref.  3:  p.  63]: 


A 


em 


G(X)2 

4  n 


PLF 


(m2  ) 


(4.20) 


PLF  =  1  ( dimensionless ) 


(4.21) 


The  effective  height  (hm)  of  the  planar  array  is  now 
estimated  by  [Ref.  6:  p.  42]: 


hem  =  2 


N 


(m) 


(4.22) 


In  (4.22),  rjQ  is  the  characteristic  impedance  of  free  space 
(1207 r)  . 

Four  two-dimensional  array  pattern  displays  are  produced 
in  this  Mathcad  application:  the  first  is  in  the  x-z  plane; 
the  second  is  in  the  y-z  plane;  the  third  display  generated  is 
a  plane  perpendicular  to  the  x-y  plane  when  <f>  =  n/4;  the  last 
display  is  a  plane  perpendicular  to  the  x-y  plane  when  <p  =  <pa. 
The  last  display  is  redundant  when  <p  is  either  0,  tt/4,  or  n/2. 

Results  obtained  with  this  application  are  the  same  as 
those  found  in  [Ref.  3:  pp.  260-274].  Directivity  was 
calculated  using  the  approximate  equation  (4.10).  An 
alternative  approximation  for  directivity  is  given  by  either 
[Ref.  3:  p.  272-273]: 
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( dimens i onl ess) 


(4.23) 


D0  =  n  COS  {  eo)  DxDy 


or 

D  -  - H -  ( dimensionless )  (4.24) 

Qa ( rads2 ) 

In  (4.23),  Dx  and  Dy  are  the  directivities  of  broadside  linear 
arrays  of  length  and  number  of  elements  Mdx,  M,  and  Ndy,  N  in 
the  x  and  y  directions,  respectively.  In  (4.24),  nA  is  the 
beam  solid  angle. 

Table  4.1  is  a  comparison  of  directivity  computed  for  a 
10x10  planar  array.  The  reference  directivity  is  derived 
using  (4.24)  while  the  calculated  directivity  is  from  (4.10). 


Table  4.1  Planar  Array  Data  Comparison 


ANTENNA 

REFERENCE 

CALCULATED 

DIMENSIONS 

DIRECTIVITY  (D0) 

DIRECTIVITY  (D0) 

M  *  N  =  10 

dx  =  dy  =  1/2  m 

eQ  =  30 0 

23.67  dB 

24.07  dB 

<Po  =  45" 

0,  =  /3y  =  -*7(272) 

Figure  4.2  and  Figure  4.3  are  comparisons  of  array  factor 
patterns  for  a  5x5  planar  array  with  patterns  from  [Ref.  3] 
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and  the  calculated  patterns.  Figure  4.2  is  a  representation 
of  the  array  factor  patterns  in  the  x-z  and  y-z  planes,  and 
the  scale  of  this  figure  has  been  increased  in  order  to 
increase  the  clarity  of  the  pattern.  Figure  4.3  is  the  array 
pattern  in  the  plane  at  <p  =  45°,  where  the  maximum  beam  is 
directed  at  <pa  =  45°.  The  array  patterns  produced  by  this 
application  are  identical  to  the  patterns  given  in  [Ref.  3]. 
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FIGURE  4.2  Array  factor  pattern  of  a  5x5  element,  equally 
excited,  planar  array  with  dx  =  dy  =  X/2  m,  0O  =  30*, <pa  - 
45*,  in  the  x-z  and  y-z  planes 
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FIGURE  4.3  Array  factor  pattern  of  a  5x5  element,  equally 
excited,  planar  array  with  dx  =  dy  =  1/2  m,  0O  =  30°,  0O  - 
45°,  in  the  plane  0  =  45° 
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V.  THE  CIRCULAR  ARRAY  AWTBWMA 


Circular  arrays,  or  ring  arrays,  are  antennas  in  which 
the  radiating  elements  are  placed  on  a  circle  with  no  elements 
positioned  inside  the  circle.  Circular  antenna  arrays  with 
elements  positioned  inside  the  circle  are  planar  arrays,  and 
these  antennas  are  addressed  in  Chapter  iv.  Circular  array 
antenna  are  used  for  radio  direction  finding,  air  and  space 
navigation,  underground  propagation,  radar,  and  sonar  systems 
[Ref.  3:  p.  274].  As  with  linear  and  planar  arrays,  circular 
arrays  with  uniformly  spaced,  equally  excited  elements  are 
narrowband  antenna  systems  due  to  the  dependance  of  the  number 
of  elements  (N)  and  the  circumference  of  the  circle  to  the 
wavelength  (A.)  [Ref.  3:  p.  277]. 

The  following  inputs  are  required  to  analyze  a  circular 
array: 

N  =  number  of  isotropic  radiating  elements 
a  =  radius  of  the  circle 
f  =  frequency  of  interest 
0„  =  direction  of  main  lobe  at  0  =  0  * 

<p0  =  direction  of  main  lobe  at  <p  —  <f>0 
I„  =  antenna  feed  current 
Z0  =  characteristic  feed  impedance 
The  first  two  inputs  are  data  primarily  obtained  from  visual 
data  such  as  photographs,  and  the  *  indicates  input  parameters 
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that  are  either  known  or  estimated.  The  planar  array  geometry 
is  illustrated  in  Figure  5.1. 


FIGURE  5.1  Geometry  of  N-element  circular  array  of 
isotropic  point  sources 


31 


This  Mathcad  application  is  written  for  N-element 
circular  arrays  with  uniformly  spaced,  equally  excited, 
isotropic  radiating  elements  on  the  x-y  plane.  The 
polarization  of  the  array  is  not  calculated  since  isotropic 
point  sources  are  used,  and  the  antenna  input  resistance  (Rjn) 
is  assumed  to  be  equivalent  to  the  radiation  resistance  (Rr)  . 
This  assumption  results  in  the  computation  of  ideal  gain  (G) . 
As  with  the  linear  and  planar  arrays,  the  radiation  pattern  of 
a  specific  element  type  is  neglected  in  the  analysis  of  the 
circular  array  since  in  normal  usage  it  will  have  little 
effect  for  a  circular  array  with  a  large  number  of  elements. 
However,  as  with  both  the  linear  and  planar  arrays,  the  user 
of  this  application  should  be  able  to  use  information 
regarding  the  type  of  element  and  the  orientation  of  the 
elements  to  narrow  the  range  of  possible  0o's  and  <P0's- 

The  normalized  array  factor  (AF)  is  calculated  for  far- 
field  observations  neglecting  mutual  coupling  between  adjacent 
elements  with  equal  amplitude  excitation  and  is  [Ref.  3:  pp. 
274-278;  Ref.  4:  350-354]; 

AF ( 6 , 4> )  =  A  y  ejikaain*con*-*o)  *«.]  ( dimensionless )  (5.1) 
N 


where 

Qn  =  n  =  1,  2,  .  .  .  N  { radians )  (5.2) 


an  =  -icasin80cos  (<|>0-Gn)  ( radians )  (5.3) 
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Equation  (5.1)  is  the  normalized  array  factor.  Equation  (5.2) 
is  the  anqular  position  of  the  nth  element  on  the  x-y  plane. 
Equation  (5.3)  is  the  phase  excitation  of  the  nth  element  with 
the  position  of  the  main  beam  directed  at  0  =  0O  and  <p  =  <p0, 
and  k  is  the  wavenumber  (2n/X)  and  a  is  the  radius  of  the 
array. 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r  *  1.6  A.  ( meters )  (5.4) 

r  i.  5  (2a)  ( meters )  (5.5) 

i  ( meters )  (5.6) 

The  minimum  distance  to  the  far-field  is  found  by  taking  the 
maximum  of  (5.4),  (5.4),  and  (5.6),  where  the  diameter  of  the 
circle  is  taken  to  be  the  maximum  dimension  of  the  antenna. 

To  determine  the  directivity  (D0)  for  a  circular  array, 
the  approximate  equations  of  radiation  intensity  (U(0,0))  for 
arrays  and  radiated  power  (Prad)  are  used.  Radiation  intensity 
is  [Ref.  3:  pp.  229-233]: 

U{  0,<J>)  =  ( AF) 2  (  W/  solid  angle)  (5.7) 

The  radiated  power  is  [Ref.  3:  pp.  28]: 

Prad  =  f2*fn/2U(B'4>)  sin(0)  d8d<t>  (V)  (5.8) 

Jo  Jo 
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In  (5.8),  Prad  is  computed  by  integrating  over  all  angles  above 
the  x-y  plane  (0  <  0  <  n / 2)  [Ref.  7:  p.  78].  This  assumes 
that  the  array  has  a  pattern  maxima  above  the  x-y  plane  with 
no  radiation  below  the  x-y  plane.  The  directivity  for  the 
circular  array  is  [Ref.  3:  pp.  229-233]: 

Da  _  .I { dimensionless )  (5.9) 

R iad 

In  (5.9),  Uraax  is  unity  and  occurs  at  0  =  0O  and  <p  =  <pc  since 
the  array  factor  is  normalized. 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4 :  p.  62 ] : 

EIRP  =  PradD0  (W)  (5.10) 

Assuming  that  the  power  delivered  to  a  circular  array  of 
N  isotropic,  uniformly  spaced  radiating  elements  is  equivalent 
to  the  radiated  power,  we  can  estimate  then  the  radiation 
resistance  (Rr)  of  the  array  as  [Ref.  3:  p.  55]: 

Rt  =  ^-T^TT  (£J)  (5.11) 

Since  the  radiating  elements  are  assumed  to  be  ideal  point 
sources,  the  input  resistance  of  a  circular  array  (Rjn)  is 
estimated  to  be  equivalent  to  the  radiation  resistance: 

Rin  =  R,  (£5)  (5.12) 
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Since  the  input  resistance  is  estimated,  the  gain  (G)  of 
a  circular  array  cannot  be  calculated,  but  an  ideal  value  can 
be  derived.  As  with  linear  and  planar  arrays,  addressed  in 
Chapter  III  and  Chapter  IV,  respectively,  the  total  antenna 
efficiency  (et)  is  [Ref.  3:  pp.  44-45]: 

£c  =  tztcd  ~  1  -  | r I2  ( dimensionless )  (5.13) 


where 


ccd  =  1  ( dimensionless ) 


(5.14) 


er  =  1  -  | r |2  ( dimensionless )  (5.15) 

r  =  — —  (dimensionless)  (5.16) 

In  (5.14),  t cd  is  the  conduction  and  dielectric  efficiency, 

which  is  unity  for  an  ideal  lossless  antenna.  In  (5.16),  ZD 

is  the  characteristic  impedance  of  a  transmission  line,  and  r 
is  the  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(5.15)  when  an  ideal  lossless  antenna  is  assumed,  and  the  gain 
of  a  circular  array  is  [Ref.  3:  pp.  43-44]: 

G  =  ttDa  ( dimensionless )  (5.17) 

G(dB)  =  10 log10  (tcD0)  ( dB )  (5.18) 
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The  maximum  effective  aperture  (km)  of  a  circular  array 
is  derived  assuming  polarization  mismatches  are  negligible 
(this  Mathcad  application  is  constructed  using  isotropic 
radiating  elements  for  which  polarization  loss  factor  (PLF) 
and  antenna  polarization  are  not  calculated)  [Ref.  3:  p.  63]: 

A  =  GJA21  PLF  ( m2 )  (5.19) 

em  4  n 

PLF  =  1  (dimensionless)  (5.20) 


The  effective  height  (h^)  of  a  circular  array  is 
estimated  by  [Ref.  6:  p.  42]: 


K.  =  2 


(m) 


(5.21) 


In  (5.21),  t}0  is  the  characteristic  impedance  of  free  space 
( 120ir)  . 

Four  two-dimensional  array  pattern  displays  are  produced 
in  this  Mathcad  application:  the  first  is  the  x-z  plane;  the 
second  is  the  y-z  plane;  the  third  is  a  plane  perpendicular  to 
the  x-y  plane  when  0  =  tt/4;  the  last  is  a  plane  perpendicular 
to  the  x-y  plane  when  <p  =  <pQ.  The  last  display  is  redundant 
when  0  is  either  0,  n/4,  or  n/2. 

The  directivity  of  a  uniform  circular  array  approaches 
the  number  of  elements  (N)  as  the  radius  (a)  of  the  array 
becomes  very  large  [Ref.  3:  p.  277].  Table  5.1  provides  the 
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directivity  of  a  uniform  circular  array  as  a  function  of  the 
radius  of  the  array. 


Table  5.1  Circular  Array  Directivity  Data 


RADIUS  (m) 

NUMBER  OF 

DIRECTIVITY  (D0)  1 

ELEMENTS 

1.59 

10 

13.70  dB 

10 

10 

12.83  dB 

20 

10 

12.89  dB 

50 

10 

11.14  dB 

100 

10 

10.33  dB 

Figures  5.2  and  5.3  provide  comparisons  of  the  array 
radiation  patterns  computed  for  a  10-element,  uniformly 
spaced,  equally  excited  circular  array  with  a  radius  of  1.59 
meters  with  those  in  [Ref.  3].  Figure  5.2  is  the  radiation 
pattern  in  the  x-z  plane,  and  the  radiation  pattern  in  the  y-z 
plane  is  shown  in  Figure  5.3.  For  these  patterns,  8a  =  <pa  = 
0° .  The  array  radiation  patterns  generated  by  this 
application  are  identical  to  the  patterns  given  in  [Ref.  3]. 
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1 


FIGURE  5.2  Array  factor  pattern  of  a  10-element,  uniformly 
spaced,  equally  excited  circular  array  with  a  =  1.59  m, 
and  0O  =  <pQ  =  0°,  in  the  x-z  plane 
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and  0O  =  <p0  =  0°,  in  the  y-z  plane 
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The  folded  dipole  is  a  center-fed  wire  antenna  that 
performs  in  a  manner  similar  to  that  of  a  linear,  center-fed 
thin-wire  antenna  with  the  added  advantage  of  a  more  easily 
matched  input  impedance  [Ref.  8:  p.  7-36].  The  folded  dipole 
antenna  consists  of  two  parallel  dipoles  connected  at  the  ends 
to  form  a  narrow  wire  loop.  The  folded  dipole  antenna  is 
illustrated  in  Figure  6.1. 


FIGURE  6.1  Geometry  of  folded  dipole  ' 

antenna  I 
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The  folded  dipole  can  be  assumed  to  behave  as  a  thin-wire 
dipole  when  the  following  are  satisfied:  the  length  of  the 
antenna  is  much  greater  than  the  diameter  of  the  feed 
conductor  (L  >  2a)  ;  the  distance  between  the  conductors  is 
greater  than  the  radius  of  the  feed  conductor  (d  >  a)  ; 
( 2wd/k )J  <  1;  all  conductors  are  assumed  to  be  perfect  [Ref. 
9:  p.  172]; 

This  Mathcad  application  calculates  the  parameters  of  a 
folded  dipole  by  assuming  the  antenna  behaves  as  a  thin-wire 
dipole  with  an  equivalent  radius  (ae)  .  The  current  maximum  is 
normalized  to  one  amp.  The  equivalent  radius  is  [Ref.  8:  p. 
7-38]  : 


a 


e 


exp 


a2ln(a)  +  b2ln(i?)  +  2a±>ln(d) 
(a  +b)  2 


( meters )  (6.1) 


This  application  provides  the  characteristics  of  the  antenna 
both  in  free  space  and  positioned  horizontally  over  the  earth. 
Although  the  folded  dipole  positioned  vertically  over  the 
earth  occurs  frequently  in  VHF/UHF  communications,  it  is  not 
addressed  in  this  research. 

The  following  inputs  are  required  to  analysis  a  folded 
dipole: 

a  =  radius  of  feed  conductor 
b  =  radius  of  second  conductor 
d  =  spacing  between  conductors 
L  =  length  of  antenna 
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h  =  height  of  antenna  above  ground  plane 
f  =  frequency  of  interest 

ecd  =  conduction/dielectric  efficiency  of  conductors 
er  =  relative  dielectric  constant  of  ground  plane 
a  -  conductivity  of  ground  plane  * 

rf£  =  distance  for  far-field  parameter  calculations 
a*  =  incoming  wave  electric  field  unit  vector  for  antenna 
in  free  space 

aa  =  unit  polarization  vector  for  antenna  in  free  space* 
ah*  =  incoming  wave  electric  field  unit  vector  for 
horizontally  positioned  antenna  * 

0hp  =  antenna  polarization  direction  for  horizontally 
positioned  antenna  * 

The  first  five  inputs  are  physical  dimensions  obtained  through 
photographs,  the  *  indicates  input  parameters  that  are  either 
known  or  estimated. 


A.  FOLDED  DIPOLE  IN  FREE  SPACE 

Assuming  the  folded  dipole  behaves  as  a  thin-wire  dipole 
with  a,  «  d  <c  L,  we  have  the  far-field  electric  field 
intensity  of  the  antenna  aligned  with  the  z-axis  as  [Ref.  3: 
p.  120]: 


E(B)  =jr\, 


r  e  ~jkI“ 


27i  r 


tt 


cos|~cos(0)|  -  COS^^Y 
sin  (0) 


( v/m ) 


(6.2) 


In  (6.2),  r)0  is  the  intrinsic  impedance  of  free  space,  equal 
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to  120tt,  and  k  =  2n/X  is  the  wavenumber.  The  normalized 
electric  field  is  found  by  dividing  (6.2)  by  the  maximum  value 
of  the  electric  field  which  occurs  at  8  =  ir/2  for  L  less  than 
about  1.25A..  Clearly,  the  radiation  pattern  of  the  folded 
dipole  is  identical  to  that  of  a  conventional  dipole  of  length 
L. 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r  *  1.6X  ( meters )  (6.3) 

r  i  5L  (meters)  (6.4) 

r  z.  (meters)  (6.5) 

The  minimum  distance  to  the  far-field  is  found  by  comparing 
the  values  of  (6.3),  (6.4),  and  (6.5)  and  selecting  the 

maximum  value.  The  distance  between  the  conductors  is  not 
taken  into  account  in  determining  the  minimum  distance  to  the 
far-field  since  it  is  assumed  that  the  folded  dipole  has  an 
over  all  effective  radius  (ae)  . 

Directivity  (D0)  for  the  folded  dipole  is  found  by 
determining  the  radiation  intensity  (U(0))  and  radiated  power 
(Pr,d)  .  Since  the  radiation  pattern  is  symmetric  in  the  x-y 
plane,  U(6)  and  Prad  are  functions  of  0  only.  Radiation 
intensity  is  [Ref.  3:  p.  28]: 
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£7(0)  =  -^-|£(0)|2  (  w/ solid  angle)  (6.6) 

2tl0 

The  radiated  power  is  [Ref.  3:  pp.  28]: 

PraCj  ~  2it  f*  u(  Q )  sin(0)  dB  (W)  (6.7) 

J  0 

Directivity  for  a  folded  dipole  is  [Ref.  3:  pp.  29): 

4  * 

Dg  =  — ~ —  ( dimensionless )  (6.8) 

° rad 

where  U,,.,  =  U(ir/2)  . 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4 :  p.  62 ] : 

EIRP  =  P radD0  (W)  (6.9) 

The  folded  dipole  is  unique  in  its  capability  to  act  as 
a  step  up  impedance  transformer  when  the  antenna  is  resonant. 
The  antenna  effectively  operates  as  an  unbalanced  transmission 
line,  and  it  can  be  analyzed  by  decomposing  the  current  into 
two  distinct  modes:  a  transmission  line  mode  [Figure  6.2]  and 
an  antenna  mode  [Figure  6.3].  [Ref.  3:  p.  341;  Ref.  4:  pp. 
205-207] 
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The  input  impedance  (Zjn)  at  the  feed  for  the  folded 
dipole  is  given  in  (6.10),  and  the  radiation  resistance  (Rin) 
for  the  antenna  is  the  real  component  of  ( Z  in) )  [Ref.  3: 
p.342] : 


= 


42tZd 


2Zd  + 


(Q) 


(6.10) 
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In  (6.10),  Zt  is  the  input  impedance  of  the  transmission  line 
mode  and  is  obtained  from  the  impedance  transfer  equation 

Zc  =  jZQ  tan|-^)  (0)  (6.11) 


where  Z0  in  (6.11)  is  the  characteristic  impedance  of  a  two- 
wire  transmission  line  and  is  given  in  (6.12)  [Ref.  3:  pp. 
341-342;  Ref.  8:  p.  7-39]: 

20  =  —  cosh^f^— )  (Q)  (6.12) 

The  input  impedance  of  a  linear  dipole  of  length  L  and 
equivalent  radius  ae  is  calculated  by  [Ref.  3:  pp.  124,  127, 
294,  and  342] : 

Zd  -  Rd  *  jXd  (Q>  (6.13) 


= 


(Q) 


(6.14) 


= 


(Q) 


(6.15) 


Rr 


~[y  +  ln(lc.L)  -  Ci  ( kL)  ] 


271 


sin  {kL)  [  5i  (  2kL) 


2  7t 


~  cos  ( kL )  (y  -*•  ln(-y) 


2Si(*L)]j 
+  Ci  (  2  kL)  -2  Ci  (  kL) 


(6.16) 
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jJpSjUcL)  +  cos  ( kL)  [2  Si(kL)  -  Si  (  2  k  L)  ]] 


(6.17) 


= 


4  n 


f  (  2kaJ))l 

sin ( kL) 

2Ci(kL) -Ci(2kL) -C^  8  Jj 

In  (6.16)  and  (6.17),  C^x)  and  St(x)  are  the  cosine  and  sine 
integrals,  respectively,  [Ref.  3:  pp.  124  and  743-746): 


Cj  (x) 


=  f'=2£iEldt  =  y.lnlx)  Y, 

Jm  x  n? 1 


2n{2n)  ! 


(6.18) 


<x>  =  L 


*  sin (x) 


*  =  Y 


n-i 


(  j  n-1  <2n- 1) 

(2/3-1)  (  2  n  - 1 )  ! 


(6.19) 


and  y  =  0.57721  (Euler's  constant). 

Mathcad  computes  summations  more  quickly  than  integrals. 
In  this  application,  the  summation  forms  of  (6.18)  and  (6.19) 
are  used  with  the  index  n  ranging  from  1  to  50  to  ensure 
proper  convergence.  The  maximum  value  of  the  argument  of 
either  C^x)  or  Sj(x)  is  x  =  4irL/X,  and  for  L  =  A,  the 
succeeding  terms  in  the  summation  rapidly  decrease  for  n  >  20. 

The  gain  (G)  of  a  folded  dipole  is  obtained  as  the 
product  of  the  antenna  reflection  efficiency  (et)  and  the 
directivity  of  (6.8)  where  the  antenna  reflection  efficiency 
is  [Ref.  3:  pp.  43-45): 

et  =  trecd  =  ecd(l-|r|2)  ( dimensionless )  (6.20) 
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where 


er  =  l  -  | r |2  ( dimensionless )  (6.21) 

r  =  — —  ( dimensionless )  (6.22) 

Rm*Zo 

In  (6.21)  er  is  the  mismatch  efficiency  and  ecd  is  the 
conduction  and  dielectric  efficiency  in  (6.22).  In  (6.22),  ZD 
is  the  characteristic  impedance  of  the  transmission  line  given 
in  (6.12).  Therefore,  the  gain  of  the  antenna  is: 

G  =  tcD0  ( dimensionless )  (6.23) 

G(dB)  =  10  log  l0(tcDo)  ( dB )  (6.24) 

The  folded  dipole  in  free  space  aligned  on  the  z-axis  is 
horizontally  polarized  with  respect  to  the  x-y  plane. 
Therefore,  given  the  direction  of  the  incoming  wave's  unit 
polarization  vector  (aj  ,  the  polarization  loss  factor  (PLF) 
is  [Ref.  3:  pp.  48-53]: 

PLF  =  |ow- 5a|2  (dimensionless)  (6.25) 

The  maximum  effective  aperture  (A^)  of  a  folded  dipole 
is  [Ref.  3 :  p. 63 ] : 

A  =  g 'l\)lpLF  (W2)  (6.26) 

^  4  n 
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The  effective  height  (hem)  of  the  folded  dipole  is  [Ref. 
6 :  p.  42  ]  : 


The  bandwidth  of  the  folded  dipole  antenna  is  dependant 
on  the  length  of  the  wires  and  has  better  bandwidth 
characteristics  than  a  conventional  dipole  of  the  same  length. 
It  can  be  assumed  that  the  bandwidth  of  the  folded  dipole  is 
essentially  the  same  as  that  of  a  conventional  dipole  of  an 
radius  a,  (a  <  ae  <  d/2)  [Ref.  3:  p.  346].  The  bandwidth  is 
increased  when  the  distance  between  the  conductors  (d)  is 
increased  and/or  when  the  conductor  radii  a  and  b  are 
increased  [Ref.  10:  p.  190]. 

Tables  6.1  and  6.2  are  comparisons  of  measured  and 
calculated  input  impedance  data  for  the  folded  dipole  as  a 
function  of  its  electrical  length.  The  radii  of  the 
conductors  (a,  b)  are  the  same  (0.0005  m)  while  the  spacing 
between  the  conductors  (d)  is  held  constant  at  0.00625  m  and 
X  =  1  m  [Ref.  8:  pp.  7-40  -  7-41]. 
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TABLE  6.1  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

LENGTH  (L) 

Re(Zjn)  n 

Re  (Zjn)  n 

k/lf 

2000 

6533  ** 

2  k/n 

2400 

2592 

3k/n 

0.0 

0.1 

Ak/n 

300 

324 

5  k/n 

1700 

1685 

6  k/n 

0.5 

1.0 

7  k/lT 

100 

100 

50 


TABLE  6.2  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

CALCULATED 

LENGTH  (L) 

Im(Zin)  n 

Im(Zin)  n 

k/n 

5000 

5437 

2k/n 

-650 

-1204  ** 

3k/n 

-90 

-85 

4k/n 

1000 

1061 

5  k/ir 

700 

793 

6k/lT 

-190 

-174 

lk/TT 

600 

633 

To  further  test  the  validity  of  the  Mathcad  application, 
a  numerical  antenna  analysis  program  was  also  used.  The 
program  used  was  ELNEC,  a  powerful,  easy-to-use  program  for 
modeling  and  analyzing  a  wide  variety  of  antenna  types 
including  ground  effects  and  parasitic  structures.  The 
fundamental  computation  portion  of  ELNEC  is  the  same  as 
MININEC  (Version  3)  ,  that  was  developed  by  the  Naval  Ocean 
Systems  Center  [Ref.  11:  pp.  5-7]. 

Table  6.3  is  a  comparison  of  the  gain  for  a  folded  dipole 
in  free  space  with  different  dimensions  [Ref.  8:  p.  3-29]. 
For  each  case,  the  radii  of  the  conductors  (0.0005  m)  and  the 


distance  between  the  conductors  (12.5a  m)  are  the  same.  As 
can  be  seen,  the  comparison  between  the  gains  obtained  from 
various  sources  is  very  good. 


TABLE  6.3  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

ELNEC 

Mathcad  1 

DIMENSION 

Gain  (G) 

Gain  (G) 

Gain  (G)  1 

1.71  dB 

1.70  dB 

1.64  dB 

keb 

I  L  =  1/2  m 

2.14  dB 

2.20  dB 

2.15  dB 

The  computed  parameters  in  Tables  6.1,  6.2,  and  6.3  are 
essentially  the  same  as  those  measured  or  computed  with  ELNEC 
with  exception  of  the  ones  marked  with  **.  These  discrepancies 
result  in  using  the  transmission  line  model  and  equivalent 
radius  to  analysis  the  folded  dipole  in  this  application. 

As  previously  mentioned,  the  folded  dipole  is  a  center- 
fed  wire  antenna  that  performs  in  a  manner  similar  to  a 
linear,  center-fed,  thin-wire  antenna.  Therefore,  the 
radiation  pattern  of  a  folded  dipole  is  also  similar  to  that 
of  a  thin-wire  antenna.  The  folded  dipole  radiation  pattern 
in  free  space  is  illustrated  in  Figure  6.4. 
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FIGURE  6.4  Radiation  Pattern  folded  dipole  in  free  space 
with  L  =  1/2  m,  a  =  b  =  0.0005  m,  and  ae  =  0.004  m 


Figure  6.5  is  a  comparison  of  the  radiation  pattern  of  a 
folded  dipole  in  free  space  obtained  using  ELNEC  and  the 
Mathcad  application.  The  patterns  are  not  identical,  but  are 
similar  and  oriented  in  the  same  direction  as  expected.  The 
difference  can  be  attributed  to  the  fact  that  ELNEC  takes  into 
account  mutual  coupling  between  the  conductors  while  the 
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Mathcad  application  models  the  folded  dipole  as  a  thin-wire, 
center-fed  dipole. 
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B.  FOLDED  DIPOLE  POSITIONED  HORIZONTALLY  OVER  THE  EARTH 

Through  the  use  of  image  theory,  the  far-field  electric 
field  intensity  of  a  folded  dipole  oriented  horizontally  above 
a  flat  earth  can  be  obtained.  The  6  component  of  the  electric 
field  of  a  horizontal  dipole  (parallel  to  the  y-axis)  in  free 
space  is  [Ref.  3:  pp.  143-144]: 


£0(0, 4>)  =jr\( 


re  ~Jkz“ 


2itr 


ft 


cos|-^sin  (6)  sin(<t>)j  -cosj-^j 
yjl  -  sin2  (0)  sin2(<|>) 


(V/m)  (6*28) 


From  image  theory,  the  0  component  of  the  electric  field 
over  ground  is  [Ref.  4:  pp.  229-235]: 

£0^(0)  =£0(0, <j>)(e^os(6)  -  rve-^cos(e))  (V/m)  (6.29) 


In  (6.28)  and  (6.29),  Ee(0,<£)  is  valid  only  above  the  ground 
(0  <  0  <  n/2)  [Ref.  3:  pp.  135-142].  In  (6.29),  h  is  the 
distance  from  the  earth  to  the  antenna  and  the  vertical 
reflection  coefficient  (Tv)  of  the  ground  is  [Ref.  4:  pp.  231- 
233]  : 


e^cosfO)  -  \j*-' z~  sin2  ( 0 ) 
e/rcos(0)  +  yjt'  x  -  sin2  ( 0 ) 


(dimensionless)  (6.30) 


where 

t'  =  e  -  j  — - ■ —  (dimensionless)  (6.31) 

r  r  2nfen 
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In  (6.30)  and  (6.31),  (e'r)  is  the  relative  complex  effective 
dielectric  constant  of  the  ground,  (er)  is  the  relative 
dielectric  constant  of  the  ground,  and  (a)  is  the  conductivity 
of  the  ground.  For  a  perfect  ground  plane,  a  -*  «  and  rv  =  1 . 

The  *  component  of  the  electric  field  of  a  horizontal 
dipole  (parallel  to  the  y-axis)  in  free  space  is  [Ref.  3:  pp. 
135-142] : 

I  e 

£<t>(0,<|>)  sJT)o  "|-y - 

2nzff 

Now  from  image  theory,  the  #  component  of  the  electric  field 
over  earth  is  [Ref.  4:  pp.  229-235]: 

£$^(0)  =E4>(e,<t>)(eJkhcoste>  +  r/1e-J'Wlcos<e))  (V/m)  (6.33) 


cos|-^sin(0)  sin(4>)j -cos|-^ 

v/l  -  sin2  (0)  sin2(<|>) 


( V/m ) 


(6.32) 


In  (6.33),  (Th)  is  the  horizontal  reflection  coefficient  of 
the  ground  and  is  [Ref.  4:  pp.  231-232]: 


cos(0)  -  yjz* ;  -  sin2  ( 0 ) 
cos(0)  +  \]z' r~  sin2  (0) 


( dimensionless )  (6.34) 


When  a  -  oo,  rh  =  -1. 

The  radiation  intensity  (Uh(0,0))  of  a  folded  dipole  over 
a  earth  must  take  into  account  both  the  0  and  #  components  of 
the  electric  field.  Therefore,  the  radiation  intensity  is 
[Ref.  3:  p.  28] : 
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Uh(6)  =  ^-[lEQgpid)  I2  +  IftyjpfG)  I2]  (w/ solid  angle)  (6.35) 

2T1o 

The  radiated  power  (Phrad)  follows  (6.7)  except  that  the  field 
is  present  only  above  the  earth.  Thus, 

Phiad  =  2nf  $  Uh(B)  sin(0)  d0  (W)  (6.36) 

J  o 


Directivity  (Dh0)  and  effective  isotropic  radiated  power 
(EIRPh)  for  a  horizontal  folded  dipole  over  the  earth  are 
identical  to  (6.8)  and  (6.9),  respectively,  with  the  results 
of  (6.36)  inserted. 

Input  resistance  (Rhin)  for  the  folded  dipole  over  the 
earth  plane  is  found  by  (where  I0  in  (6.32)  is  normalized  to 
one  amp)  [Ref.  3:  p.  124]: 


Rhin  = 


(Q> 


(6.37) 


where  the  radiation  resistance  (Rhr)  is. 


Rh 


X 


2Phgmd 

\l2o\ 


(0) 


(6.38) 


The  voltage  reflection  coefficient  (TH) ,  gain  (Gh)  , 
maximum  effective  aperture  (Ah^,)  ,  and  maximum  effective  height 
(hh.m)  for  a  horizontal  folded  dipole  over  the  earth  are 
determined  using  the  equations  for  those  parameters  for  a 
folded  dipole  in  free  space,  as  listed  above. 
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The  Polarization  Loss  Factor  ( PLFh)  for  the  antenna  over 
the  earth  can  be  determined  when  the  direction  of  the  main 
beam  of  the  polarized  antenna  (0hp)  is  estimated.  The  antenna 
polarization  vector  (aha)  is  then  calculated,  and  estimating 
the  unit  polarization  vector  of  the  incoming  wave  (crh*)  ,  the 
PLF  is  found  from  (6.25). 

Table  6.4  is  a  comparison  of  measured  [Ref.  8:  p.  3-29), 
ELNEC,  and  Mathcad  software  gains  for  a  horizontal  folded 
dipole  over  a  reflecting  sheet.  The  radii  of  the  conductors 
(0.0067  m)  are  the  same.  As  can  be  seen,  the  Mathcad  software 
gain  is  similar  to  the  measured  and  ELNEC  gains. 


TABLE  6.4  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

ELNEC 

Mathcad 

DIMENSION 

Gain  (G) 

Gain  (G) 

Gain  (G) 

d  =  L/ 25.5  m 

7.14  dB 

6.98  dB 

7.19  dB 

■bsb 

No  input  impedance  data  has  been  found  to  compare  a 
folded  dipole  orie-  1  horizontally  over  the  earth  with  that 
computed  by  this  Mathcad  application. 

Figures  6.6  and  6.7  are  the  H-plane  and  E-plane  radiation 
patterns  of  a  folded  dipole  oriented  horizontally  over  the 
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earth,  respectively.  Since  the  folded  dipole  is  similar  to  a 
thin-wire  center-fed  dipole  of  length  L,  the  illustrations 
below  are  similar  to  those  of  a  dipole  of  length  L. 


FIGURE  6.6  H-plane  radiation  pattern  of  a  folded  dipole 
with  L  =  k/2  m,  a  =  b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5 
m,  a  =  0.01  S/m,  er  =  15 
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FIGURE  6,7  E-plane  radiation  pattern  of  a  folded  dipole 
with  L  =  X/2  m,  a  =  b  =  0.0005  m,  d  =  0.00625  in,  h  =  0.5 
m,  a  =  0.01  S/m,  er  =  15 


Figures  6.8  and  6.9  are  comparisons  of  the  H-plane  and  E- 
plane  radiation  patterns  between  those  obtained  with  ELNEC  and 
Mathcad  software  of  a  folded  dipole  oriented  horizontally  over 
the  earth,  respectively.  As  expected,  the  patterns  are 
similar.  As  previosly  mentioned,  the  difference  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
coupling  between  the  conductors  while  the  Mathcad  application 
is  modeled  after  a  thin-wire,  center-fed  dipole. 
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1.2 


Figure  6.8  Comparison  between  ELNEC  and  Mathcad  of  H-plane 
radiation  patterns  of  a  folded  dipole  with  L  =  1/2  m,  a  = 
b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5  m,  a  =  0.01  S/m,  er  = 
15 


Figure  6.9  Comparison  between  ELNEC  and  Mathcad  of  E-plane 
radiation  patterns  of  a  folded  dipole  with  L  =  1/2  m,  a  = 
b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5  m,  a  =  0.01  S/m,  er  = 
15 
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VII.  THE  CAGED  DIPOLE  ANTENNA 

The  caged  dipole  antenna  is  used  in  communications  for 
its  broadband  antenna  characteristics.  Its  performance  is 
similar  to  that  of  a  thick  cylindrical  dipole.  Since  thick 
cylindrical  dipoles  are  usually  too  clumsy  to  use  except  at 
very  short  wavelengths,  thick  dipoles  are  replaced  with  a 
caged  dipole  which  consists  of  a  number  of  thin  wire  dipoles 
arranged  to  form  a  circular  cross  section.  The  antenna  is 
usually  supported  by  a  central  conductor  and  support  members 
that  provide  rigidity  but  contribute  little  to  the  electric 
properties  of  the  caged  dipole.  The  caged  dipole  (without  the 
central  conductor  and  the  support  members)  is  illustrated  in 
Figure  7.1  [Ref.  12:  pp.  172-174]. 


FIGURE  7.1  Geometry  of  caged  dipole  antenna 


This  Mathcad  application  models  the  caged  dipole  antenna 
as  a  circular  array  of  center-fed,  equally  excited,  uniformly 
spaced,  thin-wire  dipoles.  This  approach  is  taken  to  account 
for  the  number  of  conductors  that  are  present  in  the  caged 
dipole.  As  with  the  arrays  previously  considered,  mutual 
coupling  between  the  dipoles  is  neglected.  In  this  manner, 
this  application  simulates  the  characteristics  of  a  thick 
cylindrical  dipole.  The  thin-wire  dipoles  are  modeled  after 
finite  length  dipoles  in  free  space  that  have  their  length 
much  greater  than  their  diameter  (L  >  wire  diameter) ,  and  the 
dipoles  are  assumed  to  be  perfect  conductors.  Since  the 
central  conductor  and  the  support  wires  have  little  effect  on 
the  electrical  properties  of  a  caged  dipole,  mutual  coupling 
from  these  components  are  not  taken  into  account  in  this 
application.  As  seen  in  Figure  7.1,  the  inner  ends  of  a  caged 
dipole,  the  feed  points,  maybe  coned  in  order  to  reduce  base 
capacity  [Ref.  12:  pp.  172-174].  This  effect  is  also  not 
addressed. 

This  application  provides  the  characteristics  of  the 
antenna  in  free  space  as  well  as  and  oriented  both  vertically 
and  horizontally  over  the  earth. 


A.  CAGED  DIPOLE  IN  FREE  SPACE 

The  following  inputs  are  required  to  analyze  a  caged 
dipole  in  free  space: 
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N  =  number  of  conductors 


a  =  radius  of  antenna 
L  =  length  of  antenna 
f  =  frequency  of  interest 
I0  =  antenna  feed  current 
ZQ  =  input  impedance 

ecd  =  conduction/dielectric  efficiency  of  conductors 
rff  =  distance  for  far-field  parameter  calculations 
a„  =  incoming  wave  electric  field  unit  vector  for 
antenna 

aa  =  unit  polarization  vector  for  antenna 
The  first  three  inputs  are  physical  dimensions  obtained 
primarily  from  photographs,  and  the  *  indicates  input 
parameters  that  are  either  known  or  estimated. 

Assuming  the  caged  dipole  behaves  as  a  circular  array 
with  equally  excited  (IJ  ,  uniformly  spaced  elements  in  the  x- 
y  plane  with  a  radius  of  (a)  ,  we  have  the  array  factor 
(AF(0 ,0) )  [Ref.  3:  pp.  274-278]: 

AF  ( 0 ,  <t> )  =  £  lnejl*asin8cos(4~*fl>1  ( dimensionless )  (7.1) 

n  - 1 

where  the  angular  position  in  the  nth  element  in  the  x-y  plane 
is: 

=  2  it (-2),  n  =  1,  2,  .  .  .  N  ( radians )  (7.2) 
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Equation  (7.1)  is  derived  from: 

AF  (0, 4>)  =  £  jneJt*«inecos  <♦-♦„>♦«,,]  ( dimensionless )  (7.3) 

n*l 

where  the  phase  excitation  of  the  nth  element  is: 

an  =  0  =  -kasinB0cos  ($0-$n)  ( radians )  (7.4) 


In  (7.3),  and  the  phase  excitation  of  the  nth  element  (an)  is 
zero  as  the  result  of  desiring  the  main  beam  of  the  array  to 
be  directed  perpendicular  to  the  length  of  the  antenna  (0O  = 
ir/2)  .  In  (7.1),  k  is  the  wavenumber  (2n/X)  . 

Assuming  each  conductor  of  a  caged  dipole  behaves  as  a 
thin-wire  dipole,  we  have  the  far-field  electric  field 
intensity  of  a  single  conductor  aligned  with  the  z-axis  as 
[Ref.  3:  p.  120]: 


E(d)  —  jt\  InG 


2nz 


tf 


cos/— cos(0)|  -  cos/— 

\  2  I  \  2  ) 


sin  (0) 


( V/m) 


(7.5) 


In  (7.5),  r) o  =  120tt  and  is  the  intrinsic  impedance  of  free 
space. 

As  previously  discussed  in  Chapter  III,  the  total 
electric  field  of  an  array  constructed  with  identical  elements 
is  the  product  of  the  field  of  a  single  element  and  the  array 
factor  as  a  consequence  of  the  principal  of  pattern 
multiplication  [Ref.  3:  p.  207].  Therefore,  the  total 
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electric  field  of  a  caged  dipole  in  free  space  is  approximated 
as: 

Ee(0,4)  *  E(0)AF(6,4)  (V/m)  (7.6) 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25): 

r  *  l . 6  X  (meters)  (7.7) 


r  *  5 \JL2  +  (2a)2 

2  (L2  +  (2a)2) 

X 


( meters ) 


( meters ) 


(7.8) 


(7.9) 


The  minimum  distance  to  the  far-field  is  found  by  taking  the 
maximum  of  (7.7)  through  (7.9).  In  (7.8)  and  (7.9),  the 
maximum  dimension  of  the  antenna  is  the  diagonal  length  of  the 
array . 

Directivity  (D0)  for  a  caged  dipole  is  found  by 
determining  the  radiation  intensity  (U(0,0))  and  radiated 
power  (Pr,d)  .  Radiation  intensity  is  [Ref.  3:  p.  28): 

r  ^ 

U( 6,<J>)  =  — —  \Et  (0, <|> )  |2  (  w/  solid  angle)  (7.10) 

2  T|  o  1  C 

The  radiated  power  is  [Ref.  3:  pp.  28): 

Prad  =  f2*  fnU(  6,  <t>)  sin(6)  d0d<|>  (V)  (7.11) 

JO  Jo 
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Directivity  for  a  caged  dipole  is  [Ref.  3:  pp.  29]: 

D  _  ,4  n  { dimensionless )  (7.12) 

Pzad 

where  U,,,.,  =  U(tt/2,0). 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4 :  p.  62 ] : 

EIRP  =  PiadD0  (W)  (7.13) 


Input  resistance  (Rin)  for  a  caged  dipole  in  free  space 
is  found  from  [Ref.  3:  p.  124]: 

R, 


Rin  « 


sm 


2|  kL 


kL\ 

2  I 


(Q) 


(7.14) 


where  the  radiation  resistance  (Rr)  is: 

2Pr 


= 


rad 


(Q) 


(7.15) 


The  gain  (G)  of  a  caged  dipole  is  obtained  as  the  product 
of  the  antenna  reflection  efficiency  (et)  and  the  directivity, 
where  the  antenna  reflection  efficiency  is  ([Ref.  3:  pp.  43- 
45]  : 

et  =  erecd  =  ecd(l-|r|2)  ( dimensionless )  (7.16) 


where 


er  =  l  -  | r  [2  ( dimensionless )  (7.17) 
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( dimensionless ) 


(7.18) 


r  = 


Rjn-Zc 
Rin  +  Zo 


The  parameter  er  is  the  mismatch  efficiency,  £ed  is  the 
conduction  and  dielectric  efficiency  which  is  unity  for  a 
lossless  conductor,  and  Z0  is  the  characteristic  impedance  of 
a  transmission  line.  Therefore,  the  gain  of  the  antenna  is: 

G  =  ttDa  ( dimensionless )  (7.19) 


G(dB)  =  lOlog  10(tcDo)  ( dB )  (7.20) 


A  caged  dipole  in  free  space  aligned  parallel  to  the  z- 
axis  is  horizontally  polarized  with  respect  to  the  x-y  plane. 
Therefore,  given  the  direction  of  the  incoming  wave's  unit 
polarization  vector  (aw)  ,  the  polarization  loss  factor  (PLF) 
is  [Ref.  3:  pp.  48-53]: 

PLF  =  1 3 u"  od | 2  ( dimensionless )  (7*21) 

The  maximum  effective  aperture  (A^)  of  a  caged  dipole  is 
[Ref.  3 :  p. 63 ] : 

A  =  X  PLF  ( m 2)  (7.22) 

A 


The  effective  height  (h*,,)  of  the  caged  dipole  is  [Ref. 
6 :  p.  42  ]  : 


h 


em 


Rr^em 

tie 


(/n) 


(7.23) 
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The  bandwidth  of  a  thick  cylindrical  dipole  antenna  is 
dependent  on  the  ratio  of  the  length  of  the  antenna  to  the 
diameter  of  the  antenna  (L/2a)  .  When  L/2a  »  5000  the 

acceptable  bandwidth  of  a  thick  cylindrical  dipole  is  about  3% 
of  the  center  frequency,  and  when  L/2a  «  260  the  bandwidth  is 
about  30%.  Hence,  if  we  assume  a  linear  relationship  between 
bandwidth  and  L/2a,  and  since  the  caged  dipole's  performance 
is  similar  to  that  of  a  thick  cylindrical  dipole,  then  the 
bandwidth  of  a  caged  dipole  is  [Ref.  3:  p.  333]: 

BW  =  fhigh~  flou  (Hz)  (7.24) 

where  the  upper  and  lower  frequencies  are: 

■  tc  *  <«*>  (’-25) 


■ low 


=  f 


(Hz) 


(7.26) 


In  (7.25)  and  (7.25),  fc  is  the  center  frequency  and  F(L/2a) 
is  the  linear  relationship  function  between  antenna  length  and 
diameter.  This  relationship  is  obtained  by  solving  y  =  mx  + 
fc  for  m  using  xx  =  5000,  Yi  ~  3%  and  x2  =  260,  y2  =  30%.  The 
relationship  becomes: 

^jy-5.696.10-  (^)  +  0.3148  (dimensionless)  (7.27) 
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Table  7.1  is  a  comparison  of  the  gain  of  a  thick 
cylindrical  dipole  with  that  of  a  caged  dipole  (N  =  8)  in  free 
space  computed  both  by  ELNEC  and  Mathcad  software  [Ref.  8:  p. 
2-38].  As  previously  mentioned  in  Chapter  VII,  ELNEC  is  an 
easy-to-use  numerical  analysis  program  for  modeling  and 
analyzing  different  antenna  types  [Ref.  11:  pp.  5-7].  As  can 
be  seen,  the  comparison  between  the  measured,  ELNEC,  and 
Mathcad  application  values  are  similar. 


TABLE  7.1  Caged  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

ELNEC 

Mathcad  | 

DIMENSIONS 

Gain  (G) 

Gain  (G) 

Gain  (G)  | 

L  =  X/2  m 

a  =  0.025  m 

2.14  dB 

1.87  dB 

1.88  dB 

N  =  8 

Figure  7.2  is  the  radiation  pattern  of  a  caged  dipole  (N 
=  8)  in  free  space.  It  is  essentially  equivalent  to  that  of 
a  center-fed,  thick  cylindrical  dipole  having  the  same 
conductor  length  and  diameter  [Ref.  8:  p.  2-38]: 
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Figure  7.3  is  a  comparison  of  the  radiation  pattern  of  a 
caged  dipole  (N  =  8)  obtained  using  ELNEC  and  Mathcad.  The 
patterns  are  not  identical,  but  are  similar  and  oriented  in 
the  same  direction  as  expected.  The  difference  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
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coupling  between  elements  while  the  Mathcad  application 
ignores  these  effects. 


FIGURE  7.3  Comparison  of  radiation  patterns  using  ELNEC 
and  Mathcad  of  a  caged  dipole  in  free  space  with  L  =  A/2 
m,  N  =  8,  and  a  =  0.125  m 


l 


72 


B.  CAGED  DIPOLE  ORIENTED  VERTICALLY  OVER  EARTH 


The  following  known  or  estimated  inputs  are  additional 
data  required  to  analyze  a  caged  dipole  oriented  vertically 
over  the  earth.  The  height  of  the  antenna  above  ground  is  a 
physical  dimension  obtained  through  photographs: 
h  =  antenna  height  above  ground 
er  =  relative  dielectric  constant  of  earth 
a  =  conductivity  of  earth 

ov„  =  incoming  wave  electric  field  unit  vector  for 
vertical  antenna 

0vp  =  antenna  polarization  direction  for  vertically 
orientated  antenna 

As  with  the  caged  dipole  in  free  space,  the  analysis  of 
a  caged  dipole  positioned  vertically  above  the  ground  does  not 
take  into  account  the  central  conductor  and  support  members 
since  they  produce  negligible  effects  on  the  antenna's 
electrical  properties. 

From  image  theory  and  pattern  multiplication  for  arrays, 
the  total  far-field  electric  field  intensity  of  a  caged  dipole 
oriented  vertically  above  a  flat  earth  is  [Ref.  4:  pp.  229- 
235]  : 

£vc(0,<t>)  =  £v(0)  AFv(e,4>)(ejkhcos'e]  +  rve-jkhcos{e))  ( V/m )  (7.28) 

In  (7.28),  Ev(0)  is  the  far*-field  electric  field  intensity  of 
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a  finite  length  dipole  in  free  space  which  is  given  by  (7.5), 
and  AFv(0,0)  is  the  array  factor  for  a  circular  array  which  is 
(7.1).  The  total  electric  field  intensity  is  valid  only  above 
the  ground  (0  <  0  <  it/ 2)  [Ref.  3:  pp.  135-142).  The  distance 
from  the  earth  to  the  center  of  the  antenna  is  (h) ,  and  the 
vertical  reflection  coefficient  (rv)  of  the  ground  is  [Ref.  4: 
pp.  231-  233]  : 


rv(0) 


e^cost©)  -  y/e ;r  -  sin2  (0) 
c/rcos(0)  +  yjt' z  -  sin2  (0) 


(dimensionless)  (7.29) 


where 


=  Cr" 


2  n  ft. 


( dimens i onl ess) 


(7.30) 


In  (7.29)  and  (7.30),  (e'r)  is  the  relative  complex  effective 
dielectric  constant  of  the  ground,  (er)  is  the  relative 
dielectric  constant  of  the  ground,  (o)  is  the  conductivity  of 
the  ground,  and  eQ  is  the  permittivity  of  free  space  (l/(36x) 
x  10~9)  .  For  a  perfect  ground  plane,  o  -  ®>  and  rv  =  1. 

The  radiation  intensity  (Uv(0,0))  of  a  caged  dipole  over 
earth  is  [Ref.  3:  p.  28): 

jr  2 

C/v(0,  <J>)  =  — —\\Ev  (0, 4>)  |2|  (  W/  solid  angle)  (7.31) 

2t)0  j 


The  radiated  power  (Pvrad)  is  present  only  above  the  earth. 
Thus,  [Ref.  3:  p.  28): 


(7.32) 


Pvzad  =  f  2  C/v(0,4>)  sin(0)  d0d$  (f*0 
Jo  Jo 

Directivity  (Dv0)  and  effective  isotropic  radiated  power 
(EIRPv)  for  a  vertically  oriented  caged  dipole  over  the  earth 
are  given  by  (7.12)  and  (7.13),  respectively,  with  the  results 
of  (7.32)  inserted.  Uvmax  in  (7.12)  is  simply  the  maximum 
value  of  (7.31)  which,  in  Mathcad,  is  found  by  constructing 
(7.31)  as  a  vector  and  then  finding  the  maximum  value  of  the 
vector. 

Input  resistance  (Rvin)  ,  voltage  reflection  coefficient 
(TV)  ,  gain  (Gh) ,  maximum  effective  aperture  (AvM)  ,  and  maximum 
effective  height  (hvTO)  for  a  vertical  caged  dipole  over  the 
earth  are  determined  using  the  equations  for  those  parameters 
for  a  caged  dipole  in  free  space. 

The  Polarization  Loss  Factor  (PLFv)  for  the  antenna  over 
the  earth  can  be  determined  when  the  direction  of  the  main 
beam  of  the  polarized  antenna  (9vp)  is  estimated.  The  antenna 
polarization  vector  (ava)  is  then  calculated  and,  estimating 
the  unit  polarization  vector  of  the  incoming  wave  (avw)  ,  PLF 
is  found  from  (7.21). 

Table  7.2  is  a  comparison  of  the  gain  of  a  vertical  caged 
dipole  over  a  flat  earth  obtained  with  both  ELNEC  and  the 
Mathcad  application.  As  can  be  seen,  the  gains  are  similar. 
As  previously  mentioned,  the  difference  in  the  gains  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
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coupling  between  the  elements  as  well  as  ground  effects.  The 
Mathcad  application  ignores  mutual  coupling  between  elements 
and  approximates  ground  effects. 


TABLE  7.2  Caged  Dipole  Antenna  Gain  Comparison 


ANTENNA 

ELNEC 

Mathcad 

DIMENSIONS 

Gain  (G) 

Gain  (G) 

L  =  X/2  ra 

a  =  0.125  m 

N  =  8 

4.25  dB 

4.50  dB 

h  =  0.25  m 

a  =  0.01  S/m 

er  =  15 

Since  the  caged  dipole  behaves  in  a  manner  similar  to 
that  of  a  thin  dipole  except  for  bandwidth  characteristics, 
the  radiation  pattern  of  a  vertically  oriented  caged  dipole 
over  ground  will  resemble  that  of  a  thin  vertical  dipole  over 
ground.  Figure  7.4  is  a  plot  of  the  radiation  pattern  of  a 
caged  dipole  over  earth.  There  is  a  very  close  resemblance  to 
the  radiation  pattern  of  a  thin  vertical  dipole  over  ground 
[Ref.  4:  p.  234].  Figure  7.5  is  a  comparison  of  radiation 
patterns  computed  by  ELNEC  and  Mathcad  software.  As  expected, 
the  patterns  are  similar. 
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FIGURE  7.4  E -plane  radiation  pattern  of  a  caged  dipole 
with  L  =  A/2  m,  N  =  8,  a  =  0.125  m,  h  =  0.25  m,  a  =  0.01 

S/m,  and  er  =  15 


FIGURE  7.5  E-plane  radiation  pattern  comparison  between 
ELNEC  and  Mathcad  of  a  caged  dipole  with  L  =  A/2  m,  N  = 
8,  a  =  0.125  m,  h  =  0.25  m,  a  =  0.01  S/m,  and  er  =  15 
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C.  CASED  DIPOLE  ORIENTED  HORIZONTALLY  OVER  EARTH 

The  following  known  or  estimated  inputs  are  additional 
data  required  to  analyze  a  caged  dipole  oriented  horizontally 
over  the  earth: 

h  =  antenna  height  above  ground 

er  =  relative  dielectric  constant  of  earth 

o  =  conductivity  of  earth 

uhw  =  incoming  wave  electric  field  unit  vector  for 
vertical  antenna 

0hp  =  antenna  polarization  direction  for  vertically 
orientated  antenna 

As  with  the  evaluations  of  a  caged  dipole  both  in  free 
space  and  oriented  vertically  over  the  earth,  the  analysis  of 
a  caged  dipole  positioned  horizontally  above  the  ground  does 
not  take  into  account  the  central  conductor  and  support 
members  since  they  produce  negligible  effects  on  the  antenna's 
electrical  properties. 

Through  the  use  of  image  theory,  the  far-field  electric 
field  intensity  of  a  caged  dipole  oriented  horizontally  above 
a  flat  earth  can  be  obtained.  The  0  component  of  the  electric 
field  of  a  horizontal  dipole  (parallel  to  the  y-axis)  in  free 
space  is  [Ref.  3:  pp.  143-144]: 
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Ie-m„ 

cosj 

|-^sin  (0)  sin(<t>)j 

|  -  cosj 

(f)l 

J  \Q  O  ^ 

2nr  ff 

Vl  -  sin2  (0)  sin2(<j>)  j 

(V/m) 


(7.33) 


The  0  component  of  the  total  electric  field  over  ground 
is  [Ref.  4:  pp.  229-235]: 

EQhc  (0,  <J>)  =E6h(6,4>)  AFh(e,4>)(ejkhcosi»  -  rve^khcosW)  ( V/m)  (7.34) 


The  array  factor  (AFh(0,0))  in  (7.34)  for  a  horizontal  array 
is  determined  by  [Ref.  3:  pp.  274  and  776]: 

AFhid,  4>)  =  £  Inejkas  in6sin*sin*«  ( dimensionless )  (7.35) 

n  *  1 

In  (7.33)  through  (7.35),  the  electric  field  (E0h(0,0)),  the 
total  electric  field  (E0ht(0,0)),  and  array  factor  (AFh(0,0)) 
are  valid  only  above  the  ground  (0  <  0  <  n/2)  [Ref.  3:  pp. 
135-142].  In  (7.34),  (h)  is  the  distance  from  the  earth  to 
the  center  of  the  antenna,  and  the  vertical  reflection 
coefficient  (Tv)  of  the  ground  is  given  by  (7.29). 

The  *  component  of  the  electric  field  of  a  horizontal 
dipole  (parallel  to  the  y-axis)  in  free  space  is  [Ref.  3:  pp. 
135-142]  : 

I  e'jkz“ 

E4>h(8,<b)  =jr\o~!L— - 

2nr{{ 

From  image  theory,  the  *  component  of  the  total  electric  field 
over  earth  is  [Ref.  4:  pp.  229-235]: 


cosi^sin  (0)  sin(4>)j- 


cos 


■(f) 


\J  1  -  sin2  (0)  sin2(4>) 


(V/m) 


(7.36) 
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£4>/it(0,<t>)  =  £$A(0,<t>)  AF*(0,<|))(e^co*(9)  ♦  r^e-^0*'6*)  iv/m)  (7.37) 


In  (7.37),  ( AFh (0,0) )  is  given  by  (7.35),  and  (Th)  is  the 

horizontal  reflection  coefficient  of  the  ground  and  is  [Ref. 
4:  pp.  231-232]: 


cos(0)  -  ^e'r-  sin2  (0) 
cos(0)  +  Je'z  -  sin2  ( 0 ) 


( dimensionless )  (7.38) 


When  a  -  <»,  rh  =  -1. 

The  radiation  intensity  (Uh(0,0))  of  a  caged  dipole  over 
a  earth  must  take  into  account  both  the  0  and  *  components  of 
the  electric  field.  Therefore,  the  radiation  intensity  is 
[Ref.  3 :  p.  28 ] : 

t/A(0,<t>)  =3—  [|£0ht  (0,<1>)  |2  +  |£4>/2c  <0 ,  <|>)  |2]  {W/ solid  angle)  (7.39) 

O 


The  radiated  power  (Phrad)  is  present  only  above  the  earth. 
Thus, 

Phiad  =  f*' f*/2Uh(Q,$)  sin(e)  d0d<t>  (W)  (7.40) 


The  remaining  parameters  (i.e.,  Dh0,  EIRPh,  Rhin,  etc.) 
for  a  horizontal  caged  dipole  over  the  earth  are  calculated 
using  the  equations  for  these  parameters  listed  previously. 

Table  7.3  is  a  comparison  of  the  gain  of  a  horizontal 
caged  dipole  over  a  flat  earth  obtained  with  ELNEC  and  the 
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Mathcad  application.  As  can  be  seen,  the  gains  are  similar. 
As  previously  mentioned,  the  difference  in  gains  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
coupling  between  the  elements  as  well  as  ground  effects.  The 
Mathcad  application  ignores  mutual  coupling  between  elements 
and  approximates  ground  effects. 


TABLE  7.3  Caged  Dipole  Antenna  Gain  Comparison 


As  previously  mentioned  with  regard  to  the  vertically 
oriented  caged  dipole  over  the  earth,  the  caged  dipole  behaves 
in  a  manner  similar  to  that  of  a  thin  cylindrical  dipole. 
Therefore,  the  radiation  patterns  of  a  horizontally  oriented 
caged  dipole  are  similar  to  those  of  horizontal  thin 
cylindrical  dipole  which  in  turn  resembles  that  of  a  thin 
dipole.  Figures  7.6  and  7.7  illustrate  the  H-plane  and  E- 
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plane  patterns,  respectively,  of  a  caged  dipole  over  the  earth 
[Ref.  4:  pp.  231-237].  The  radiation  patterns  of  the 
horizontally  orientated  caged  dipole  over  a  flat  earth  are 
very  similar  to  those  of  a  horizontal  thin  dipole  over  earth. 
Figures  7.8  and  7.9  compare  the  radiation  patterns  between 
those  computed  with  ELNEC  and  those  computed  by  Mathcad 
software.  As  expected,  the  patterns  are  similar. 


FIGURE  7.6  H-plane  radiation  pattern  of  a  caged  dipole 
with  L  =  1/2  m,  N  =  8,  a  =  0.125  m,  h  =  1.0  m,  a  =  0.01 
S/m,  and  er  =  15 
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S/m,  and  cr  =  15 


FIGURE  7.8  Comparison  between  ELNEC  and  Mathcad  of  H-plane 


radiation  patterns  of  a  caged  dipole  with  L  =  1/2  m,  N  = 


8,  a  =  0.125  m,  h  =  1.0  m,  a  =  0.01  S/m,  and  er  =  15 
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FIGURE  7.9  Comparison  between  ELNEC  and  Mathcad  of  E-plane 
radiation  patterns  of  a  caged  dipole  with  L  =  1/2  m,  N  = 
8,  a  =  0.125  m,  h  =  1.0  m,  a  -  0.01  S/m,  and  er  =  15 
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VIII.  PARABOLIC  REFLECTORS 


Many  communications  and  radar  operations  require  very 
large  values  of  antenna  gain  (G)  that  are  difficult  to  achieve 
with  a  single  device.  One  of  the  most  popular  ways  to 
construct  a  very  high  gain  antenna  is  to  reflect  the 
electromagnetic  energy  of  a  relatively  small  antenna  off  a 
metallic  parabolic  dish  in  the  far-field  of  the  feed  device 
[Ref.  13:  p.  12-2].  If  the  source  antenna  is  located  at  the 
focus  of  the  parabola  and  if  the  radius  (a)  of  the  parabolic 
dish  periphery  is  large  with  respect  to  wavelength  (k) ,  then 
the  reflector  may  be  approximated  as  a  uniform  phase  aperture 
antenna.  Gain  in  excess  of  30  dB  is  common  for  this  type  of 
structure  [Ref.  4:  pp.  423-424].  The  geometry  of  a  parabolic 
reflector  is  illustrated  in  Figure  8.1,  where  the  parameters 
include  the  focal  distance  (f),  the  distance  (p)  from  the 
origin  to  any  point  (P)  on  the  reflector  surface,  and  the 
distance  (r')  from  the  origin  to  the  projection  of  point  (P) 
onto  the  z  =  0  plane.  The  primed  angles  in  Figure  8.1 
correspond  to  the  parameters  associated  with  the  feed  horn. 
The  Mathcad  application  user  should  note  that  the  focal 
distance  should  exceed  the  minimum  distance  to  far-field  for 
the  feed  antenna. 
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FIGURE  8.1  Parabolic  Reflector  Geometry 


The  Mathcad  parabolic  reflector  applications  are  unique 
in  that  the  software  is  not  written  as  a  stand  alone  package. 
Rather,  the  reflector  software  is  run  as  an  addendum  to  the 
programs  that  already  exist  for  various  antennas  [Ref.  2:  pp. 
4-15,  67-114].  Four  types  of  feeds  are  analyzed:  helical, 
planar  equiangular  spiral,  pyramidal  horn,  and  conical  horn. 
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The  following  additional  inputs  are  required  to  analyze 
parabolic  reflectors: 

foe  =  reflector  focal  length 

are£  =  radius  of  the  mouth  of  the  reflector 

E0  =  electric  field  scale  factor 

eb  =  blockage  efficiency 

esp  =  spar  efficiency 

r££2  =  secondary  far  field  observation  distance 
x,y,z  =  polarization  loss  factor  coordinates 
aw  =  incoming  wave  unit  polarization  vector 
tl,  t2  =  radiated  power  increments  for  <p  and  0 
il  =  secondary  field  increments 
N,M  =  summation  increments 

The  first  two  inputs  are  physical  dimensions  obtained 
primarily  through  photographs,  inputs  with  the  *  are 
parameters  that  are  either  known  or  estimated,  and  the  last 
inputs  are  only  used  to  affect  the  processing  time  of  the 
application. 

In  order  to  calculate  the  far-field  patterns  of  a 
parabolic  reflector,  the  electric  field  at  the  aperture  of  the 
device  must  be  known.  Since  the  electric  field  patterns  of 
interest  are  in  the  far-field,  it  is  assumed  that  the  circular 
aperture  created  by  the  reflector  is  located  at  the  origin  and 
is  parallel  to  the  z  =  0  plane.  The  Cartesian  components  of 
the  aperture  electric  field  are  expressed  as  [Ref.  14:  pp. 
121-123;  Ref.  15:  pp.  29-31]: 
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(V/m) 


(8.1) 


E«(r',0')  =  [-  Eef(0'(r')  ,0')cos(0') 

+  E,f(0'  (r')  ,0')sin(<*>')  ]eJk(p(r  )'2£oc) 


Eay(r',0')  =  [EM(0/(r/)  ,0/)sin(0/) 

(8.2) 

+  E*f(0'(r')  ,0')cos(0')  ]eJk('>(r')-2£oc)  (V/m) 

In  (8.1)  and  (8.2),  (E„ (r'  ,  0' ) )  and  (Eay(r'  ,0' ))  are  the 

aperture  electric  fields  and  k  is  the  wavenumber  (2n/X) . 

The  relationships  between  (0'),  (p) ,  and  (r'  )  needed  to 
evaluate  (8.1)  and  (8.2)  are  expressed  as  [Ref.  4:  p.  426]: 

0'  =  2 tarr1^  2§oc]  ( radians )  (8-3) 

p(r')  ■  -i lj°£)2 +{:')*  (jn)  (8.4) 


The  electric  field  from  the  entire  parabolic  assembly  is 
commonly  referred  to  as  the  secondary  pattern.  Referring  to 
Figure  8.1,  the  electric  field  in  the  -z  half-plane  is 
considered  the  secondary  pattern.  The  secondary  field  of  a 
parabolic  reflector  is  calculated  using  the  electric  field 
integral  equation  solution  and  computing  the  component  vector 
potentials  (Px(0,0)),  (Py(0,0))  as  follows  [Ref.  4:  p.  426]: 


*«.♦>« 


E  (r'  0/)  (0)(cos  ($')  cos  (t> 


sin  (♦')  sin 


<♦>)] 


*  r'dr'dfy' 


{V/m) 


(8.5) 
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Double  integrals  in  Mathcad  take  significantly  longer  to 
converge  to  an  answer  than  do  single  integrals  [Ref.  15:  p. 
200].  To  reduce  the  processing  time  of  the  applications, 
(8.5)  and  (8.6)  were  each  separated  into  two  single  integrals. 
Therefore,  for  the  applications,  (8.5)  and  (8.6)  are  each 
computed  as: 

P(0,4>)  =  f2*A  (0 , 4> » 4>')  ^4/  ( V/m )  (8.7) 
J  0 


where 

A(0,<|>,4>')  ={*'•'  Ea{r',$')  eJ*',sin(e)[cos(«/)cos(*>  <•  Ein(V)sin(*)]r/dr/  (8.8) 
J  o 

Any  contribution  from  feed  system  back  lobes  is  ignored, 
and  the  component  vector  potentials  calculated  in  (8.7)  and 
(8.8)  are  applied  to  the  subsequent  formulas  to  compute  the 
secondary  field  [Ref.  4:  p.  383]: 

-j  h-a  -jkr  [C2 

Eq(Q,Q>)  =  - [Px(0,<t>)  cos  (<j>)  +  Py (0, <t>)  sin (4>) ]  (V/m)  (8.9) 

t  /**>  ~  jkl f  £  2 

£^(0,40  =~^ - f-Px(0,<J>)  sin(4>)  +  Py(0, 4>)  cos  (<|>)  jcos  (0)  (8.10) 
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The  radiation  intensity  (U(0,0)),  radiated  power  (Prad)  , 
and  directivity  (D0)  of  the  secondary  field  pattern  are 
expressed  using  the  intrinsic  impedance  of  free  space  ( rj0  = 
120jt)  as  follows  [Ref.  3:  pp.  27-30]: 

£/(0,  <J>)  =  [|ge  (6,  <t>)  |2  (0, <J>)  |2]  (W/ solid  angle)  (8.11) 

2r\0 

Prad  =  (2n  f'/2U  (6, 4>)  sin  (Q)cBd4>  (w)  (8.12) 

Jo  Jo 

D0  =  —  ( dimensionless )  (8.13) 

Piad 

In  (8.13),  (Uroax)  is  the  maximum  value  of  radiation  intensity, 
obtained  by  creating  (Umax)  as  a  matrix  and  finding  the  maximum 
value  of  the  matrix. 

The  computer  processing  time  required  to  evaluate  (8.12) 
is  prohibitively  long.  Thus,  the  parabolic  reflector  Mathcad 
applications  use  the  following  numerical  approximation  to 
evaluate  radiated  power  [Ref.  3:  pp.  37-42]: 

-  (^)(-^)^  (s  (8-14) 

In  (8.14),  (N  and  M)  are  user  selected  integers.  This 
approximation  for  Prad  reduces  the  processing  time  of  the 
applications  substantially.  As  N  and  M  increase,  the  results 
obtained  with  (8.14)  converge  to  that  obtained  with  (8.12). 
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When  N  =  M  =  3,  the  processing  time  for  only  Pr4d  of  the 
helical  feed  reflector  option  is  approximately  four  hours  on 
a  33  MHz,  386  machine,  and  when  N  and  M  is  increased  to  5,  the 
time  to  calculate  Pr4d  for  the  same  application  is  ten  hours. 
For  the  applications,  we  chose  N  =  M  =  10  to  evaluate  Pr4d  with 
reasonable  accuracy.  With  this  choice,  it  takes  approximately 
one  day  for  the  application  to  run,  but  this  is  much  less  than 
the  thirty-six  hours  it  takes  if  (8.12)  is  used. 

Parabolic  reflectors  are  members  of  the  aperture  antenna 
family.  The  gain  of  an  aperture  antenna  is  the  product  of  its 
aperture  efficiency  (eap)  and  directivity.  Therefore,  the  gain 
of  a  parabolic  reflector  is  written  as  [Ref.  6:  p.  575]: 

G  =  tapDa  ( dimensionless )  (8.15) 

Aperture  efficiency  is  the  product  of  several  separate 
terms,  including  spillover  efficiency  (es) ,  taper  efficiency 
(et)  ,  phase  (or  random  surface  error)  efficiency  (ep)  , 
polarization  efficiency  (ex) ,  blockage  efficiency  (eb) ,  spar 
efficiency  (esp)  ,  and  ohmic  efficiency  (eohmic)  . 

Spillover  efficiency  (es)  is  the  fraction  of  power 
radiated  by  the  feed  system  that  is  intercepted  and  reflected 
by  the  parabolic  dish.  Spillover  efficiency  is  the  most 
important  efficiency  term  for  any  feed  system  and  is  defined 
as  [Ref.  6:  pp.  583]: 
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f  Kf  °Uf  (0')  sin  (0')  dQ'dfy' 

es=  — - — - -  ( dimensionless )  (8.16) 

f  Tuf(0')Sin<e')d6'd*' 

Jo  Jo 


In  (8.16),  (Uf(9'  ))  is  the  radiation  intensity  function  of  the 
feed  antenna  and  (0O)  is  the  value  of  (8.3)  when  (r' )  equals 
the  radius  of  the  mouth  of  the  reflector  (aref)  . 

Taper  efficiency  (et)  accounts  for  the  lack  of  amplitude 
uniformity  of  the  feed  pattern  on  the  surface  of  the 
reflector.  Taper  efficiency,  which  is  a  dimensionless 
quantity,  is  expressed  as  [Ref.  3:  pp.  626-627]: 


et  =  2cot‘ 


2 


)r 

Jo  \j 

tanl 

(1) 

2 

d0'd<t>' 

/ 

2*f0o4 nUf(B'  .  /  /  , 

/  - £ — sin(0O  dd'dfy' 

3  J°  pr  ad 

(8.17) 


Phase  efficiency  (ep)  is  a  function  of  the  power  loss 
that  occurs  if  the  field  at  the  mouth  of  the  dish  is  not  in 
phase  at  every  point  in  the  aperture.  Phase  efficiency  is 
approximately  given  by  [Ref.  4:  pp.  434-435]: 


-««610'sa„f 

cp  ~  e  ( dimensionless ) 


(8.18) 


The  power  that  is  lost  due  to  cross-polarized  electric 
fields  in  the  aperture  of  the  reflector  determines 
polarization  efficiency  (ex) .  Polarization  efficiency  is 
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difficult  to  precisely  calculate,  but  can  be  reasonably 
approximated  as  0.98  [Ref,  4:  p.  435]. 


Blockage  (eb)  and  spar  (esp)  efficiencies  are  determined 
by  the  presence  of  feed  structures,  support  struts,  and 
possible  signal  processing  equipment  in  front  of  the 
reflector's  aperture.  Precise  computation  of  blockage  and 
spar  efficiencies  requires  an  extensive  method  of  moments 
simulation  and  is  not  attempted  by  the  Mathcad  parabolic 
reflector  applications.  However,  if  the  radius  (af)  of  the 
blocking  structure  at  the  focal  point  and  the  number  of  spars 
(N)  that  are  (A/2)  thick  is  known,  blockage  efficiency  and 
spar  efficiency  can  be  estimated  using  Tables  10.1  and  10.2, 
respectively  [Ref.  4:  p.436]. 


TABLE  8.1  Blockage  Efficiency 


af/ aref 

0.05 

0.10 

0.20 

0.99 

0.96 

0.84 

TABLE  8.2  Spar  Efficiency 


aref  —  5/ A 

aref  =  100/A 

N  =  3 

0.95 

0.99 

t-* 

« 

o 

o 

N  =  4 

0.94 

0.99 

1.00 
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Ohmic  efficiency  (eohmc)  for  a  reflector  is  usually  very 
high  but  difficult  to  accurately  calculate.  The  Mathcad 
parabolic  reflector  applications  assume  ohmic  efficiency  is 
0.98. 

With  all  reflector  efficiency  terms  known  or  estimated, 
the  aperture  efficiency  ( e ap)  ,  effective  aperture  (A,)  ,  and 
isotropic  radiated  power  (EIRP)  are  given  by  [Ref.  3:  pp.  623- 
630] : 

eap=  tstctptxtbtSpt  ohmic  ( dimensionless )  (8.19) 

Cap^xef)2  ^2)  (8.20) 

EIRP  =  PzadD0  ( W )  (8.21) 

The  magnitude  of  the  feed  current  (IQ)  for  the  source 
antenna  cannot  be  estimated  based  on  measured  geometry. 
However,  if  the  feed  current  is  normalized  to  1  amp,  radiation 
resistance  (Rr)  for  the  entire  reflector  structure  is: 

Rr  =  <Q)  (8.22) 

Uol 

In  (8.22),  (Prad)  is  the  radiated  power  from  the  aperture  of 
the  reflector,  not  of  the  feed  antenna. 

The  Mathcad  parabolic  reflector  applications  use  the 
radiation  resistance  computed  for  the  parabolic  reflector 
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antenna  system  to  determine  effective  height  (h«„)  by  [Ref.  6: 
p.  42 ] : 


(8.23) 


The  bandwidth  of  a  parabolic  reflector  is  assumed  to  be 
the  same  as  its  feed  antenna. 

The  polarization  loss  factor  (PLF)  of  the  parabolic 
reflector  for  a  user  selected  point  in  the  secondary  far-field 
and  a  specific  wave  unit  polarization  vector  (aw)  is  found 
using  the  electric  field  components  of  (8.9)  and  (8.10)  as 
follows  [Ref.  3:  51-53]: 

PLF  =  |o,/da|2  (dimensionless)  (8.24) 


where  the  unit  polarization  vector  of  the  antenna  is: 


-  /  X  aKEx+  ayEy  +  azEz 

oa(x,y,z)  =  -  ,  y  -y  , 

yj\E(x,y,  z)  I 


(dimensionless)  (8.25) 


and  the  Cartesian  components  of  the  electric  field  are  [Ref. 
3:  p.  776]: 

Ex  =  E9 cos  (0)  cos  (4>)  -  fi^sin  (<J>)  ( v/m )  (8.26) 


Ey-  Eecos  (0)  sin  (<J>)  +  E^c  os(<|>)  (v/m)  (8.27) 


Ez=  -£‘esin(0)  (v/m) 


(8.28) 
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All  of  the  secondary  field  parameters  calculated  by  the 
Mathcad  parabolic  reflector  applications  assume  that  the 
observation  point  is  in  the  far-field.  Thus,  all  of  the 
subsequent  formulas  must  hold  for  the  application  calculations 
to  be  valid  [Ref.  4:  pp.  24-25]: 

r  ^  1 .6  A.  (m)  (8.29) 


r  z  10  aref  (m) 


(8.30) 


r  * 


8(*r  eff 

k 


(dimensionless) 


(8.31) 


Table  8.3  is  a  comparison  of  reflector  directivity 
computed  by  the  Mathcad  applications  to  predictions  made  for 
an  identical  antenna  with  unity  aperture  efficiency  [Ref.  3: 
p.  634].  As  can  be  seen,  the  comparisons  vary  by 
approximately  8  dB.  This  can  be  attributed  to  the  fact  that 
Prad  was  approximated  using  (8.14)  instead  of  (8.12),  and  also 
to  the  fact  that  the  numerical  tolerance  variable  (TOL  in 
Mathcad' s  math  menu)  for  the  application  was  increased  to  0.5 
verses  0.001  in  order  to  reduce  processing  time.  This 
increase  in  TOL  reduces  both  the  accuracy  of  (8.14)  and  the 
processing  time  of  the  application. 
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TABLE  8.3  Parabolic  Reflector  Data  Comparison 


ANTENNA 

PREDICTED 

CALCULATED  1 

PARAMETER 

DATA 

DATA 

DIRECTIVITY 

42  dB 

50.2  dB 

Figure  8.2  is  a  comparison  of  measured  data  fed  by  an 
optimum  horn  with  that  calculated  by  the  Mathcad  application 
fed  by  an  equiangular  spiral  for  a  parabolic  reflector  (aref  = 
0.320  m,  f  =  0.147  m)  operating  at  35  GHz  [Ref.  14:  p.  126]. 
As  can  be  seen,  the  comparison  is  good.  The  Mathcad  pattern 
of  Figure  8.2  is  not  symmetric  which  can  be  attributed  to  the 
number  of  increments  used  to  increase  the  speed  of  processing. 

The  reflectors  fed  by  pyramidal  and  conical  horn  antennas 
are  not  included  in  this  report.  The  applications  have  been 
developed  and  are  included  on  the  diskette  containing  the 
software.  The  complexity  of  the  software  results  in 
prohibitively  long  processing  time  for  the  two  applications 
(i.e.,  the  pyramidal  horn  fed  reflector  application  was  still 
calculating  Prad2  after  seven  days  when  run  on  a  50  MHz,  486 
personal  computer  using  Mathcad  4.0). 


FIGURE  8.2  Comparison  of  Parabolic  Reflector  E-Plane 
Electric  Field  Patterns 


IX.  REMARKS  AND  CONCLUSIONS 


This  thesis  and  the  accompanying  software  applications 
are  intended  to  provide  the  Naval  Maritime  Intelligence  Center 
(NAVMARINTCEN)  with  software  for  the  analysis  of  various 
antenna  types  where  information  regarding  the  antenna  is 
restricted  primarily  to  the  antenna  type  and  the  physical 
dimensions  of  the  antenna.  The  goal  of  this  research  was  to 
provide  user  friendly  software  that  an  engineer  familiar  with 
basic  antenna  types  can  easily  use  and  interpret.  Although 
the  software  was  developed  for  that  purpose,  the  complexity  of 
formulas  that  apply  to  many  of  the  antenna  types  of  necessity 
reduced  the  simplicity  of  the  applications. 

The  software  applications  developed  are  written  in 
Mathcad  3.1  and  are  compatible  with  any  IBM  personal  computer 
that  supports  Mathcad  3 . 1  or  Mathcad  4.0  for  Windows.  Many  of 
the  applications  require  extensive  processing  time.  In  order 
to  reduce  processing  time,  which  results  from  the  complexity 
of  various  integrals,  approximate  numerical  techniques  were 
implemented  to  evaluate  the  integrals  and  the  numerical 
tolerance  of  the  more  complicated  applications  was  reduced. 
Even  with  the  processing  time  reduced,  some  programs  require 
up  to  five  days  to  complete  when  run  with  Mathcad  4.0  on  a  50 
MHz,  486  personal  computer.  Applications  that  did  not  exceed 
a  day  to  complete  were  run  on  a  33  MHz,  386  personal  computer 
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with  Mathcad  4.0.  All  applications  were  run  using  Mathcad  4.0 
since  Mathcad  4.0  is  approximately  twice  as  fast  as  Mathcad 
3.1.  As  previously  mentioned,  the  applications  were  written 
in  Mathcad  3.1  due  to  sponsor  requirements.  Mathcad  3.1  is 
compatible  with  Mathcad  4.0,  but  Mathcad  4.0  is  not  compatible 
with  Mathcad  3.1. 

For  this  thesis,  a  number  of  the  antennas  could  not  be 
adequately  researched  through  the  use  of  texts.  In  some 
cases,  such  as  the  caged  dipole,  very  little  information  on 
the  antenna  could  be  found  in  either  texts  or  in  professional 
journals.  For  the  caged  dipole,  the  antenna  is  modeled  as  a 
circular  array  of  center-fed,  equally  excited,  uniformly 
spaced,  thin-wire  dipoles  to  take  into  account  the  number  of 
conductors  that  are  present  in  the  caged  dipole.  This  model 
provides  good  results  when  compared  to  predicted  (the  caged 
dipole's  performance  is  similar  to  that  of  a  thick  cylindrical 
dipole)  and  ELNEC  results. 

It  was  found  that  Mathcad  is  not  the  ideal  tool  to 
analyze  antennas  due  to  its  long  processing  time.  An 
additional  drawback  of  Mathcad  is  that  when  a  function  is  in 
the  beginning  of  an  application  and  the  same  function  is 
placed  later  in  the  application,  Mathcad  does  not  store  the 
values  of  the  function,  but  re-evaluates  the  function.  This 
re-evaluation  increases  processing  time  significantly, 
especially  when  the  function  itself  takes  a  very  long  time 
(hours  vice  minutes)  to  reach  a  numerical  solution. 
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Nevertheless,  the  Mathcad  applications  developed  in  this 
research  allow  the  user  to  analyze  several  antenna  types 
(i.e.,  arrays  with  isotropic  point  sources,  parabolic 
reflectors  with  the  various  feed  antennas)  that  could  not 
easily  be  evaluated  using  some  antenna  numerical  analysis 
programs  such  as  NEC,  ELNEC ,  or  MININEC.  For  the  other 
antenna  types  researched,  the  antenna  numerical  analysis 
programs  mentioned  would  be  easier  to  setup,  run,  and  provide 
parameters  closer  to  those  of  a  real  antenna.  The 
disadvantage  of  this  approach  is  that  it  requires  an  ability 
to  evaluate  antennas  using  a  program  such  as  ELNEC,  while  the 
Mathcad  applications  dc  not  require  any  advanced  programming 
knowledge  on  the  part  of  the  user. 
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LINEAR  ARRAY 

MATHCAD  SOFTWARE  ARRAY_LN . MCD 

This  linear  array  application  is  based  on  N  isotropic  point  sources 
uniformly  spaced  and  aligned  on  the  z-axis.  Mutual  interference  between 
adjacent  sources  are  ignored  when  calculating  antenna  parameters.  Since 
isotropic  point  sources  are  used  as  the  radiating  elements,  polarization 
of  the  antenna  can  not  be  determined.  The  gam  can  only  be  idealized 
when  it  is  assumed  that  the  antenna  input  resistance  is  equivalent  to 
the  radiation  resistance.  Bandwith  calculations  are  not  carried  out  as 
a  result  of  the  narrowband  characteristics  of  arrays,  and  in  large  part 
to  the  dependance  lr.ter-element  spacing  (d)  has  with  the  wavelength  (1)  . 
The  following  limits  should  be  adhered  to  when  analyzing  Broadside, 
End-Fire,  and  Phased  (scanning)  antennas: 

0  <  0O  <  *  (Scanning  Antennas) 

0O  =  «/ 2  (Broadside  Antennas) 

0O  =  0  or  7i  (Ordinary  End-Fire  Antennas) 

(Note:  Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore, 

parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn ) . ) 

The  linear  array  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  »  indicate  parameters  that  can  be  computed  but 
their  results  are  idealized  as  the  result  of  using  isotropic  point 
sources  as  the  radiating  elements): 


1  =  Wavelength 

rmm  =  Minimum  Distance  to  the  Far-Field 

P  =  Phase  Difference  Between  Adjacent  Elements 
Hi  =  Array  Factor  Phase  Shift 
AFn  =  Normalized  Array  Factor 

U(0)  =  Radiation  Intensity 

prad  “  Radiated  Power 
D0  =  Directivity 

EIRP  =  Effective  Isotropic  Radiated  Power 

Rr  =  Radiation  Resistance 

Rin  =  Antenna  Input  Resistance  * 

T  =  Voltage  Reflection  Coefficient  * 

Et  =  Antenna  Efficiency  * 

G  =  Gam  * 

Aem  =  Maximum  Effective  Aperture  * 
hem  =  Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 


N  =  Number  of  Isotropic  Radiating  Elements 

f  =  Frequency  of  Interest 

d  =  Distance  Between  Ad}acent  Elements 

0O  =  Direction  of  Main  Lobe 

I0  =  Antenna  Feed  Current 

ZQ  *  Characteristic  Feed  Impedence 


N  :=  10 


Enter  input  data  here 


(elements ) 


d:=0  25 


(meters) 


f:=3  108  (Hz) 


0o  :=  0 


( radians) 


Io  :=  1  0  (A) 


Zo  :=75  0  (Q) 


Define  constants  and  calculate  wavelength: 

(Q) 


(m_1) 


(dimensionless) 
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2.9979  10*  (meters/ sec) 


ia-c 


1=0  999 


(meters/cycle) 


(meters/cycle) 


7)0  —  120-71 


k:=l-n 


ccd  :=  1 


( radians ) 


Calculate  linear  array  parameters: 
Define  angular  offset  9: 

0:*  (o-  icr “0  (2  n-  10'^ 

Minimum  Distance  to  Far-Field  rm1n'. 

r  (meters) 


r,:»5  Nd 

.  2  (N  d)2 


(meters ) 


(meters ) 


rmin :«  max(  r) 


rmin  =  12  509 


(meters) 


(meters ) 


Array  Factor  AFp, : 
(3  :=  -  kdcos(0o) 


<K0)  :=  k  d  cos(0)-f  p 


( radians) 


( radians) 


(dimensionless) 
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Radiation  Intensity  U: 
U(0)  :=(AFn(6»2 


Radiated  Power  P 


Prad  :*2  it 


rad-L 


U(0)  sin(0)d0 


Prad  =  1.256 


Directivity  Dn: 

Umax  •=  1 


Do:= 


4  n  Umax 
Prad 


Do  =  10  007 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :■  Prad  Do 


EIRP  =  12  566 


(W/solid  angle) 


(W) 


(W) 


(W/solid  angle) 


(dimensionless) 


(dimensionless ) 


(W) 


(W) 
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Radiation  Resistance  Rr: 


n  .  2  Prad 
Rr  :* - 

(|lo|  )2 


(Q) 


Rr  =  2  512 


<Q) 


Input  Resistance  Rln: 

Rin  :*  Rr  (0) 


Rin  =2.512 


<Q) 


Voltage  Rejection  Coeffecient  [: 

T  :r  (dimensionless) 

Rin  +  Zo 


r=~0  935 


(dimensionless) 


Reflection  Efficiency  St- : 

rt  =ecd  [  1  -  ( |r|  )2]  (dimensionless) 


d  =0  125 


(dimensionless) 
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Gain  G: 


i  G  r  El  Do 


(dimensionless) 


GdB  =  10  log(G)  (dB) 


G  =  1  255 


(dimensionless) 


GdB  =0.986  (dB) 


Maximum  Effective  Aperture  Aom: 

.  .  G  J.2  ,  2> 

Aem  := (m^) 

4  k 


Aem=0  1  (m2) 


Maximum  Effective  Height  h^,m: 


hem  :=  2 


J 


Rr  Aem 

no 


(m) 


hem  =  0  052 


(m) 


) 
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N  =  10 
0o=  0 


X(0)  :*  |AFn(0)|  cos/  0  ♦  — ] 


NORMALIZED  ARRAY  TACTOR  POLAR  PLOT 

Number  of  Elements: 

Direction  of  Main  Beam  (radians): 


Y(0)  | AFn( 0)|  sinld 


PLANAR  ARRAY 

MATH CAD  SOFTWARE  ARRAY  PL.MCD 


This  planar  array  application  is  based  on  M  and  N  isotropic  point 
sources  uniformly  spaced  and  aligned  on  the  x-axis  and  y-axis 
respectively.  Mutual  interference  between  adjacent  elements  are 
neglected  when  calculating  antenna  parameters.  Since  isotropic  point 
sources  are  used  as  the  radiating  elements,  polarization  of  the  antenna 
can  not  be  determined.  The  gain  can  only  be  idealized  when  it  is 
assumed  that  the  antenna  input  resistance  is  equivalent  to  the  radiation 
resistance  and  the  total  antenna  efficiency  is  calculated  for  a  ideal 
lossless  antenna. 

(Note:  Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore, 

parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 

The  planar  array  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  can  be  computed  but 
their  results  are  idealized  as  the  result  of  using  isotropic  point 
sources  as  the  radiating  elements): 

1  =  Wavelength 

Lx  =  Antenna  Dimension  in  x-direction 
Ly  =  Antenna  Dimension  in  y-direction 
rmin  =  Minimum  Distance  to  the  Far-Field 

Px  =  Progressive  Phase  Shift  Between  Adjacent  Elements  in  x-direction 

py  =  Progressive  Phase  Shift  Between  Adjacent  Elements  in  y-direction 

i(ix  =  Array  Factor  Phase  Shift  in  x-direction 

ijiy  =  Array  Factor  Phase  Shift  in  x-direction 

AFn(0, 4)  =  Normalized  Array  Factor 

U<0»  =  Radiation  Intensity 

prad  =  Radiated  Power 
D0  =  Directivity 

EIRP  =  Effective  Isotropic  Radiated  Power 
Rr  =  Radiation  Resistance 

Rin  =  Antenna  Input  Resistance  * 

T  =  Voltage  Reflection  Coefficient  * 

=  Antenna  Efficiency  * 

G  =  Gain  * 

Aem  =  Maximum  Effective  Aperture  * 
hem  =  Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 

M  =  Number  of  Isotropic  Radiating  Elements  in  x-direction 
N  =  Number  of  Isotropic  Radiating  Elements  in  y-direction 
f  =  Frequency  of  Interest 

dx  =  Distance  Between  Adjacent  Elements  in  x-direction 

dy  =  Distance  Between  Adjacent  Elements  in  y-direction 

0O  =  Direction  of  Main  Beam  at  0  =  0o 

=  Direction  of  Main  Beam  at  ♦  =  4° 

I0  =  Antenna  Feed  Current 

ZQ  =  Characteristic  Feed  Impedence 

3  =  Number  of  Increments  in  Degrees  for  Far-Field  Radiation  Pattern 


Enter  input  data 

here 

M  :=5 

(elements ) 

dx  :=0.S 

(meters) 

II 

(elements ) 

dy  =0.5 

(meters) 

f:=3  10® 

(Hz) 

0o:=- 

6 

( radians) 

lo  =  10 

(A) 

.  _  * 

*°-4 

( radians ) 

Zo  :=  75  0 

( Q) 

j  :=  180 

(degrees ) 
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Define  constants  and  calculate  wavelength: 


t)0  =  120  x 

k  — 

1 

ocd  -  1 

Calculate  linear  array  parameters: 

Define  angular  offset  8 andj: 

e  =  (o-  10-^. 10-4+i^  (2  n-  1(T4) 

^  :*  (0-  ia^  lff4+  — J  (2  k-  1CT4) 


c  *29979  108  (meters/sec) 


c 

l'-3-  (meters/cycle) 


1=0  999  (meters/cycle) 


Minimum  Distance  to  Far-Field  rml n : 
rQ:=16  1  (meters) 

ri  :=5  >/(M  dx)2+  (N  dy)2  (meters) 

2 


U 


(Mdx)24(Ndy)2J 


l  ^  r  *5  ) 


rmin  :=  max(r)  (meters) 


rmin  =  25  018 


( Q ) 


tm~  1 ) 


(dimensionless ) 


( radians ) 


( radians) 


(meters ) 
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Array  Factor  AFn  : 


|3x  :=-k  dx  sin(0o)  cos(*o) 


( radians ) 


Py  --k  dy  sin(8o)  sm($0) 


( radians ) 


♦>  ;  =  k  dx  sin(0)  cos{$)-»-  px 


( radians ) 


<l»y(0,4)  =k  dy  sin(0)  sin(^)  +  py 


(radians) 


AFn(0 


/M  \  /w  ' 

sml— i|ix(0.*)j  sinl—  i^yCO.4) 


(dimensionless) 


Radiation  Intensity  U: 


U(0,*)  -=AFn(M)  AFn(0,*) 


(W/solid  angle) 


Radiated  Power  P^: 


r  2 «  I 


Prad  :=  I  U(8J)  sin(0)  d0d* 

•'0  0 


Prad  =  0.206 


Directivity  Dn: 


Umax  :*  1 


(W/solid  angle) 


Do  ■- 


4*  Umax 


Prad 


( dimensionless ) 


Do  =61  09 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRP: 

EIRP  =  Prad  Do  (W) 

EIRP  =  12  566  (W) 


Radiation  Resistance  Rr: 

R r:  =  i^l 
(llol  )2 


(D) 


Rr  =0411 


(Q) 


t 


Input  Resistance  R^r,: 

Rin-Rr  (Q) 


Rin  =0411 


(Q) 
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Voltage  Rejection  Coef fecient  f : 
Rin  -  Zo 


r  •=■ 


Rin  +  Zo 


r  =-0  989 


Reflection  Efficiency  c,. : 

cl  :* ecd- [  1  -  ( |r|)2j 


ct  =0  022 


Gain  G: 


Maximum  Effective  Aperture  Aom: 

,2 


Aem  :■ 


G  1 
4 n 


(m2) 


Maximum  Effective  Height  horn : 


•*  2 

IJO 


.  —  i  *  v*  Aem 

hem  *2  | -  (m) 


(dimensionless) 


( dimensionless ) 


(dimensionless) 


(dimensionless ) 


G:  =  ct  Do  (dimensionless)  GdB  :*  10  log(G)  (dB) 


G  =  1.326  (dimensionless)  GdB  =  1  225  (dB) 


Aem  =0  105  (m^) 


hem  =0021  (m) 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  (x~z)  plana 


Number  of  Elements  along  x-axis: 


Number  of  Elements  along  y-axis: 


Direction  of  Main  Beam: 


00  =  0  524 


*o=0  785 


( radians ) 


( radians ) 


X1(0)  =  |Afn(0,O)|  cosjo-  Y  1( 0)  =  |AFn(0,O)|  sm|0- 


0  =  0C 


'  > _ i _ _  '  /  ✓  v  f 

'  /  *  \  , 

'  '  ■  '  ,  *  „  , 

,'v  '  V  ,  ' 

\  \  •  ;  ,  •>  *  >  > 

V  .  .’ft  '+  -  1  » 

\  1  (  '  '  '  <  I 

^  \  s'  \  I 

l  l  t  (  '  .  s  ^  , 

'  *  ^  ,  s<  \  .-*• 

'  V  *  ,  ,  ■*.  '  ^  >■  '  ,  " 

"  -  "  •  „  '  Y 

'  .  4  S  ,  ,  ,  /  ^  I  -  "  * 


TWO-DIMENSIONAL  ANTENNA  FATTEKN:  *  ■  »/4 


Number  of 

Elements 

along  x-axis: 

M  =5 

Number  of 

Elements 

along  y-axis: 

N  =5 

Direction 

of  Main 

Beam: 

0o=O524 

( radians) 

*>=0  785 

( radians ) 

X3(0)  :=  |AFn|e,^J  |  cos|e-  Y3(9) :■  AFn|e,-^J  |  smjo- 

6=0° 


xref  ,X3(8) 
m 
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M  =5 


TWO-DIMENSIONAL  ANTENNA  PATTERN:  ± 

Number  of  elements  along  x-axis: 
Number  of  Elements  along  y-axis: 
Direction  of  Main  Beam: 


N  =5 

0o=O.524  (radians) 

^0=0  785  (radians) 


X4(6)  :=  |AFn(0,4o)|  cos|e-  Y4(0)  :=  |AFn(e,*>)|  sinje- 


=  oc 
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CIRCULAR  ARRAY 

MATH CAD  SOFTWARE  ARRAY  CI.MCD 


This  circular  array  application  is  based  on  N  isotropic  point  sources 
uniformly  spaced  on  the  x-y  plane  along  a  circular  ring  of  constant 
radius.  Mutual  interference  between  adjacent  elements  are  neglected 
when  calculating  antenna  parameters.  Since  isotropic  point  sources  are 
used  as  the  radiating  elements,  polarization  of  the  antenna  can  not  be 
determined.  The  gain  can  only  be  idealized  when  it  is  assumed  that  the 
antenna  input  resistance  is  equivalent  to  the  radiation  resistance  and 
the  total  antenna  efficiency  is  calculated  for  a  ideal  lossless  antenna 


(Note:  Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore, 

parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 


The  circular  array  antenna  Mathcad  application  will  compute  thefollowing 
parameters  (Items  with  *  indicate  parameters  that  can  becomputed  but 
their  results  are  idealized  as  the  result  of  usinqisotropic  point 
sources  as  the  radiating  elements) : 


1  =  Wavelength 

rmin  =  Minimum  Distance  to  the  Far- Field 
un  =  Phase  Excitation  of  the  nth  Element 
AF (0, 4)  =  Normalized  Array  Factor 
U(0>)  =  Radiation  Intensity 

Pracj  =  Radiated  Power 

D0  =  Directivity 

EIRP  =  Effective  Isotropic  Radiated  Power 

Rr  =  Radiation  Resistance 

R^n  =  Antenna  Input  Resistance  * 

T  =  Voltage  Reflection  Coefficient  * 

£t  =  Antenna  Efficiency  * 

G  —  Gain  * 

Aem  =  Maximum  Effective  Aperture  * 
hem  =  Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 

N  =  Number  of  Isotropic  Radiating  Elements 
f  =  Frequency  of  Interest 
a  =  Radius  of  circle 
0O  =  Direction  of  Main  Beam  at  0  =  0o 

40  =  Direction  of  Main  Beam  at  ♦  = 

I0  =  Antenna  Feed  Current 

ZQ  =  Characteristic  Feed  Impedence 

3  =  Number  of  Increments  in  Degrees  for  Far-Field  Radiation  Pattern 


Enter  input  data  here 


N  :=  10  (elements) 


a  •=  ■ 


JO 

2-n 


(meters ) 


f:=3  10*  (Hz) 


0o~O  (radians) 


lo  :=  1 .0  (A) 


Zo  :=75  0  (Q) 


♦©:=  0  (radians) 


j  := 360  (degrees) 


Define  constants  and  calculate  wavelength: 

c:=  2.9979  10*  (meters /sec)  ijo:sI20n  (D) 


1 


c 

f 


(meters/cycle) 


k 


2  n 


X 


(m  1 ) 
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Define  constants  and  calculate  wavelength: 

X  =  0 999  (meters/cycle)  scd  :*  1  (dimensionless) 

Calculate  linear  array  parameters: 

Define  angular  offset  9_and_±: 


0:*(o-  Id4) , 

Kt) 

(2  n-  id4) 

( radians ) 

♦  =(o-  ior4), 

io-4+i-?j 

(2  a-  id4) 

( radians ) 

Minimum  Distance  to  Far-Field  r^in: 

r0  =1.6-1  (meters)  ^=5  (2  8)  (meters) 


V 


■  _2(2a)z 


(meters ) 


rmin  :=  max(r)  (meters) 


rmin  =  20.278  (meters) 


Array  Factor  AF: 
n  = 1  N  (elements) 


*  -  t  n 

0  =  2n — 
"  N 


( radians ) 


a 

n 


k  a  sin(0o)  cos^o-  ^ 


(radians) 
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Array  Factor  AF: 


AF(M) 


(k »iin(9) 


a 


(dimensionless) 


Radiation  Intensity  U: 


U(0,*)  :=  AF(0,$)  AF(0,*)  (W/ solid  angle) 


Radiated  Power  Prari: 


Prad  := 


'0 


U(0,4)  Sin(0)d0c4 

0 


Prad  =0  536 


(W) 


(W) 


Directivity  D^: 


Umax  '•=  1 


(W/solid  angle) 


Do  •= 


4  n  Umax 
Prad 


(dimensionless) 


Do  =23.427 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :*  Prad  Do 


(W) 


EIRP  =  12.566 


(W) 
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Radiation  Resistance  Rr: 


R r:» 


2  Prad 

( llol  )2 


(Q) 


Rr  =  I  073 


(Q) 


Input  Resistance  Rin: 


Rin  :*  Rr 


(Q) 


Rin  -1.073 


(Q) 


Voltage  Relection  Coeffecient  [: 

~  (dimensionless) 

Rin  +  Zo 


r  =-0  972 


(dimensionless ) 


Reflection  Efficiency  g» : 

d  =ocd  [l  -  ( |r|)2]  (dimensionless) 


d  =0.056 


(dimensionless) 
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Gci in  G : 


G-  =  ctDo  (dimensionless)  GdB  •=  10  log(G)  (dB) 


G  =  1.303  (dimensionless)  GdB  =  1.149  (dB) 


Maximum  Effective  Aperture  Aom: 

r  ,2 

Aem :« -  (m^)  Aem  =0.104  (m^) 

4* 


Maximum  Effective  Height  hQm: 

hem  •=  2  [fi1 .  ^*0  (m)  hem  =0034  (m) 

V  no 
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TMO-DIMEN3IOMJU.  ANTENNA  PATTERN:  (X-X)  plan* 


Number  of  Elements  in  x-y  plane: 

N  =  10 

Direction  of  Main  Beam: 

00=0 

( radians ) 

*o  =  0 

( radians ) 

X1(0)  :=  |aF(0,O)|  cos|e- 

y  i(0) 

:=  |AF(0,O)|  sin^-  -ij 

0 

o 

II 

<E> 
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TWO-DIMENSIONAL  ANTENNA  PATTERN :  (y-«)  plan* 


Number  of 

Elements  in  x-y  plane: 

N  =  10 

Direction 

of  Main  Beam: 

6o=0 

( radians ) 

♦o  =0 

(radians) 
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TWO-DIMENSIONAL  ANTENNA  PATTERN :  j  -  t/4 


Number  of  Elements  in  x-y  plane: 
Direction  of  Main  Beam: 


X3(8) 


N  =  10 


8o=0  (radians) 

4o=0  (radians) 


Y3(8)  := 


6=0° 


TWO-DIMENSIONAL  ANTENNA  PATTERN:  $  -  ja 


Number  of  Elements  in  x-y  plane: 

N  =  10 

Direction  of  Main  Beam: 

O 

II 

& 

( radians ) 

o 

ii 

z 

( radians ) 

X4(0)  :=  |AF(6,*o)|  cosje- 

Y4(0) 

:*  |AF(0.*o)|  sm|0-  -I 

\  2/ 

e  =  o° 
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FOLDED  DIPOLE 

MATHCAD  SOFTWARE  FOLDED. MOD 


To  provide  impedance  matching  to  coaxial  input  cables  and  maintain  good 
directional  pattern  characteristics,  a  folded  dipole  is  used  instead  of 
a  single  dipole  element.  The  geometry  of  a  folded  dipole  is  such  that 
the  spacing  between  dipole  elements  (d)  is  much  less  than  the  wavelength 
(1)  and  the  length  of  an  element  (L) . 

This  Mathcad  application  calculates  the  parameters  of  a  folded  dipole 
aligned  vertically  on  the  z-axis  and  a  folded  dipole  positioned 
horizontally  above  the  earth.  For  the  calculations,  the  antenna  is 
assumed  to  behave  as  a  cylindrical  dipole  with  equivalent  radius  (ae) . 

(Note:  Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore, 

parameters  do  not  contain  subscripts  (i.e.,  Rin  will  be  written  Rin) . ) 

The  folded  dipole  antenna  Mathcad  application  will  compute  the  following 
parameters  (*  indicates  antenna  parameters  that  are  calculated  for  a 
folded  dipole  in  free  space  and  positioned  horizontally  over  the  earth) : 


1  =  Wavelength 

ae  =  Equivalent  Radius 

Rr  =  Radiation  Resistance 
Rd  =  Input  Resistance  of  Linear  Dipole 
Xr  =  Reactance  of  Linear  Dipole 
Xd  =  Input  Reactance  of  Linear  Dipole 

Zd  =  Input  Impedance  of  Linear  Dipole 

Zq  =  Characteristic  Impedance  of  Two-wire  Transmission  Line 
Zt  =  Characteristic  Impedance  Transform  of  Folded  Dipole 
zin  =  Characteristic  Impedance  of  Folded  Dipole 
ZQ  =  Maximum  Current  of  Folded  Dipole 

rmin  =  Minimum  Distance  to  the  Far-Field 
E  (9)  =  Electric  Field  * 

En  (0)  =  Normalized  Electric  Field  * 

U (0)  =  Radiation  Intensity  * 

UmaX  =  Maximum  Radiation  Intensify  * 
prad  =  Radiated  Power  * 

D0  =  Directivity  * 

EIRP  =  Effective  Isotropic  Radiated  Power  * 

R^n  =  Antenna  Input  Resistance  * 

T  =  Voltage  Reflection  Coefficient  * 
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Antenna  Efficiency  * 

Antenna  Polarization  Vector 


Et  = 

°a  = 

PLF  =  Polarization  Loss  Factor  * 

G  =  Gain  * 

Aem  =  Maximum  Effective  Aperture  * 
hem  -  Maximum  Effective  Height  * 

crp  =  Relative  Complex  Effective  Dielectric  Constant  of  Ground 
rv  =  Vertical  Reflection  Coefficient  of  Ground  Plane 
rh  =  Horizontal  Reflection  Coefficient  of  Ground  Plane 


The  following  known  or  estimated  data  must  be  entered: 

f  =  Frequency  of  Interest 

a  =  Radius  of  Conductor 

b  =  Radius  of  Conductor 

d  =  Distance  Between  Conductors 

L  =  Length  of  Antenna 

h  =  Height  of  Antenna  Above  Earth 

Ecd  =  Conduction/Dielectric  Efficiency 

er  =  Relative  Dielectic  Constant  of  Ground 

o  =  Conductivity  of  Ground 
rff  =  Far  F’ielcl  Distance 

i  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 
ow  =  Incoming  Wave  Electric  Field  Unit  Vector  For  Free  Space  Antenna 

°wh  =  INcoming  Wave  Electric  Field  Unit  Vector  For  Horizontal  Antenna 
8ph  =  Antenna  Polarization  Direction  For  Horizontal  Antenna 
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Enter  input  data  here 


ow  :*  0 1 

o/ 

(dimensionless ) 

M 

awh  0 

1, 

f:=3  10* 

(Hz) 

Oph  :*0 

a  :*  0  0005 

(meters) 

rff  :=  1  105 

b  :=  a 

(meters ) 

i  :=  90 

d  :=  12  5  a 

(meters) 

ecd  :=  1 

L  :*  0  5 

(meters) 

o  :*25 

h  :=  5 

(meters ) 

a  :=  15 

( dimensionless ) 


( radians ) 


(meters) 


( radians ) 


(dimensionless) 


(siemens/meter) 


(dimensionless) 
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Define  constants: 


oa  •= 


(dimensionless) 


c  :*  2  9979  10*  (meters/sec ) 


Tp  :=  1 20  n 


1 


(meters/cycle) 


1  »  0.999 


(meters/cycle) 


7:=  0.57721 


eo . 


—  10' 9 
36  n 


(F/m) 


z  =0 . 2 i 


Define  angular  offset  $: 


e  :* 0,—  .  2  JI 

> 


Define  cosine  [Ci(x)}  and  sine  [Si(x)]  integrals: 


n:=  1.  50 


Ci(x)  .=  0.57721  + 


^)-S 


(- 1 )"~  1  x2n 
(2  n)  (2  n)! 


n 


Si(x)  =  ]T 


(-D 


(2  n- 


n 


(Q) 


(m  1 ) 


(dimensionless ) 


( radians ) 


(radians) 


n-  1  x(2n-  1) 

1 )  (2  n-  1)! 
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FOLDED  DIPOLE  ALIGNED  ON  THE  Z-AXIS  IN  FREE  SPACE: 


Equivalent  radius  aa: 


ae  =  exp 


a2  ln(a)  +  b2  ln(b)  +  2  a  b  ln(d) 


(a+  b) 


ae  =0  002 


Electric  field  E(0): 


e-j  kfff 

cos 

cos(0)j  -  cos 

m 

•  J 

2  *  rff 

1  sin(0)  f 

En(0)  : 


..  E(0) 


Minimum  Distance  to  Far-Field  r. 


minj. 


v=161 


(meters ) 


r,:=5L 


V* 


2d  r 


(meters ) 


rmin  :=  max(r)  (meters) 


rmin  =  2.5 


(meters ) 


(meters) 


(meters) 


(V/m) 


(V/m) 


(meters) 
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Radiation  Intensity  U (0)  : 


U(8)  :*——•(  |E(0)|  )2 


Radiated  Power  Pran= 


Prad  :*  2  a  U(0)sin(0)d0 
•  0 


Prad  =  36.64 


Directivity  D^: 


41 


4-ff' Umax 


Do  =  1  641 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP Prad  Do 


EIRP  =60.132 


(W/solid  angle) 


(W/solid  angle) 


(dimensionless) 


(dimensionless) 
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Radiation  Resistance  Rr: 


Rr :-iE. 
2  n 


ln(kL) -  Ci(k  L)  +  —  sin(k  L)  (Si(2  k  L)-  2  Si(kL)) 
2 


+  -  cos(k  L) 
2 


Mt) 


+■  Ci(2  k  L)  -  2  Ci(kL) 


Rr  =  73  281 


Input  Radiation  Resistance  of  linear  dipole  R^ : 
Rr 


Rd  •'=  • 


Rd  =73  281 


Reactance  of  linear  dipole  Xr: 


Xr:=J32 

4i 


2  Si(kL)+ cosfkL)  (2  Si(kL)-  Si(2  kL»  . 

- 1  sin(k  L)  (2  Ci(kL)-  C^2  kL)-  cJ2  k 


r| 


Xr  =  43  142 


Input  reactance  of  linear  dipole  X^: 
Xr 


Xd:=- 


Ht)) 


Xd  =43  142 


(Q) 


(Q) 


(Q) 


(Q) 


<Q) 


(Q) 


<Q) 


(Q) 


136 


Zd  =Rd  +  j  Xd 


(Q) 


Zd  =  73.281  +43.142i 


(Q) 


Characteristic  impedance  of  two-wire  transmission  line  Zn: 
Zo  =—acosh[ — (Q) 

*  \2  jTb/ 


Zo  =  302.312 


(Q) 


Characteristic  impedance  transform  for  folded  dipole  Z«- : 
Zt:=j  Zotan/^j  (Q) 


Zt  =-2.747  105i 


(Q) 


Characteristic  impedance  of  folded  dipole  Zjr,: 

.  4  Zt  Zd 

Zin  := - 

2  Zd+Zt 


(Q) 


Zin  =293.307  +172.465i 


(Q) 


Input  Resistance  Rin: 

Rin  =Re(Zin)  (Q)  Rm  =  293.307  (Q) 
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Voltage  Rejection  Coeffecient  [: 

r  _  Rm-  Zo  (dimensionless)  r=-0015  (dimensionless) 

Rin  ♦  Zo 

Reflection  Efficiency  z*  : 

a  :*ecd[l-  (|r|)2]  (dimensionless)  ct  =  1  (dimensionless) 

Gain  G: 

G:*£t  Do  (dimensionless)  GdB  :=  10  Iog(G)  (dB) 

G  =  1641  (dimensionless)  GdB  =2  15  (dB) 

Polarization  loss  factor  PLF: 

PLF  =  flow  oa|)  (dimensionless)  PLF  =  I  (dimensionless) 


Maximum  Effective  Aperture  flpm: 
Cl2 

Aem  :■ - PLF  (m2) 

4  * 


Aem  =  0  13  (m2) 


Maximum  Effective  Height  hom: 

.  .  /Rr  Aem  ,  , 

hem  »  2  I -  (m) 


hem  =0  318 


(m) 


T|0 
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HORIZONTALLY  POSITIONED  FOLDED  DIPOLE  OVER  THE  EARTH: 


Define  angular  offsets  81  and  82: 

( radians ) 


.  n  ji  n  n 
61  > — , —  + 


2  2  2  i  2 


02  j 

2  2  2  i 


Relative  complex  effective  dielectric  constant  £rp: 


£rp:=er- j 


2  n  f  so 


4  component  of  E- field  Ej: 


E4(01)  :*j  vf>- 


e'  krfT 
2  it  rff 


sm(0l )  sm(0)j  -  cos|^-tj 


J 1  -  sin(01 ) 


2  sin(0)2 


E41  :»j  no- 


e-J  krfr 
2  n  rff 


cos|-~  s|n(02^  sin(0)j  -  cos|—'j 


J 1  -  H*2)  2  sin(O)2 


0  component  of  E- field  E0: 

s|i±  su*8l )  sin(^)  -  cos(li) 


E0(01 )  :=  j  jf>- 


e-J  k  rff 

2*rff 


r 

1 


I  -  sin(01 )  sin 


E01z  :*j  ^ 


e-j  krff 

cosj 

[kL 

i  2 

sin  ^02^  sinj 

2  n  rff 

J 

1 1  -  sin^02^ 

HiT 

( radians ) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 
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Vertical  plane  wave  reflection  coefficient 


rvi 


,_apcos(0I )-  1 

Jap-  sin(01 ) 

ap  cos(01 )  4  < 

erp  cos(02  J  -  J 

Jap  -  sin(01 ) 

xp-  sin  (02  J  2 

ap  cos(02^  4  Ji 

rp-  sin  (02^ 2 

2 


2 


(dimensionless ) 


(dimensionless) 


Horizontal  plane  wave  reflection  coefficient 


rh(01 )  -C0^61)-  yop-  sm(81)Z 

cos(01 )  4  -Jap -  sin(01  )2 


Fhl 


(02  J  -  Jap  -sin  (02  J‘ 
(02 }  4  Jap-  sin (02 


(dimensionless ) 


(dimensionless) 


E-field  (4  component)  with  ground  plane  reflection  Ji*2£± 
E4gp(01).=E*(01)  (e1  khc“<ei>4  Fh(01 )  e*j  khc0,(9l>)  (v/m) 


E*gplz's 


kh-coi 


W 


4  rhl  e 


-j  khcoa(92^ 


m 


(V/m) 


E4n(01 )  ■=  lL 

max(E*gpl) 


(V/m) 
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(V/m) 


E8gp(01):*E8(ei)  (e*  k  h “*<9l) -  rv(01 )  e'J  khco‘<8l)) 


E0gplz:: 


j  k  h 


-  Fvl  e 

z 


■j  k  h  co*(02^J 


E8n(81 )  :=  ll 

max(  E0gp  1 ) 


Radiation  Intensity  Uh(6l): 

Uh(01 )  =  — [( |E*gp(01  )|  )2+  ( |E0gp(01  )|  )2] 
2 VP 


Radiated  Power  Ph 


radi. 


Phrad  —  2  « 


Uh(01)  sin(01)d01 


Phrad  =  70  415 


(V/m) 


(V/m) 


(W/solid  angle) 


(W) 


(W) 


Directivity  Dh,-,: 

Ugp_max^ :=  ~~  [  +  (|E0gp]J)  J 


(W/solid  angle) 


Uhmax  =  max(  Ugpmax ) 


(W/solid  angle) 


Directivity  Dh^: 


Dho  = 


4  *  Uhmax 
Phrad 


(dimensionless) 


Dho  =  3.938 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRPh: 


EIRPh  :=  Phrad  Dho 


(W) 


EIRPh  =  277  3 


(W) 


Input  Resistance  Rh,„: 


Rhin 


2  Phrad 


(Q) 


Rhin  =  140  83 


(Q) 


Voltage  Relection  Coeffecient  TH: 

Rhin-  Zo 
!H  :* - 

Rhin  Zo 


(dimensionless ) 


rH  =-0.364 


(dimensionless) 


Reflection  Efficiency  ghj- : 

Eht:=8cd  [l-(|m|  )2]  (dimensionless)  dit  =  0.867  (dimensionless) 
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G3 1  n  Gti  i 


Gh;=ditDho  (dimensionless)  GhdB  =  10  log(Gh)  (dB) 


Gh=3415  (dimensionless)  GhdB  =  5.334  (dB) 


Polarization  loss  factor  PLFv: 


Ehmax  =  |  E0n(  0ph )  |  (V/m) 


Lobes  =  floor|—  +  lj 


oah 


0 

0 

Ehmax 


Lobes 


(dimensionless ) 


PLFh : 


2 


(dimensionless ) 


PLFh  =7  736  10"4 


(dimensionless ) 


Maximum  Effective  Aperture  Avpm: 

Gh  l2 

Ahem  — - PLFh  (m2)  Ahem  =  2  099"  10  (m2) 

4  « 


Maximum  Effective  Height  hvPTt,: 


hhem  :* 


JRhinAhem 


(m) 


hhem  =0018  (m) 
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E-Plane  Radiation  Pattern  (x-z  plane)  of  Horizontal  Folded  Dipole  Over 
The  Earth: 
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CAGED  DIPOLE  ANTENNA 
MATHCAD  SOFTWARE  CAGED_DI.MCD 

This  caged  dipole  application  is  modeled  after  a  circular  array  with  the 
radiating  elements  as  dipoles.  The  dipoles  are  modeled  after  conductors 
with  infinitely  small  radii  and  antenna  length  greater  than  1/4.  Mutual 
interference  between  adjacent  elements  are  neglected  when  calculating 
antenna  parameters. 

(Note:  Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore, 

parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 


The  caged  dipole  antenna  Mathcad  application  will  compute  the  followin' 


VERTICAL  CAGED  DIPOLE  IN  FREE  SPACE: 

1  =  Wavelength 

rmin  =  Minimum  Distance  to  the  Far-Field 

AF (0,4)  =  Normalized  Array  Factor 
E  (0,^)  =  E-Field  of  Finite  Length  Dipole 

£t(0,4)  =  Total  E-Field 

Etn<*W  =  Total  Normalized  E-Field 
U  (0,4)  =  Radiation  Intensity 

prad  ~  Radiated  Power 

Dc  =  Directivity 

EIRP  =  Effective  Isotropic  Radiated  Power 

R r  =  Radiation  Resistance 
R^n  =  Input  Resistance 

T  =  Voltage  Reflection  Coefficient 
ct  =  Antenna  Effeciency 
G  =  Gain 

PLF  =  Polarization  Loss  Factor 

Aem  =  Maximum  Effective  Aperture 

hem  =  Maximum  Effective  Height 

flow  =  Low  Frequency 

fhigh  =  High  Frequency 
BW  =  Bandwidth 

VERTICAL  CAGED  DIPOLE  OVER  EARTH: 

erp  =  Relative  Complex  Effective  Dielectric  Constant 

AFv (01,4)  =  Normalized  Array  Factor 
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Ev (81, 4)  =  E-Field  of  Finite  Length  Dipole 

rv  <81 )  =  Vertical  Plane  Wave  Reflection  Coefficient 

Evt(81,4)  =  Total  E-Field 

Evtn*®1'^)  =  Total  Normalized  E-Field 

Uv(01,4)  =  Radiation  Intensity 

Pvrad  =  Radiated  Power 

Dvc  =  Directivity 

EIRPv  =  Effective  Isotropic  Radiated  Power 
Rvr  =  Radiation  Resistance 

Rvin  =  Input  Resistance 

IV  =  Voltage  Reflection  Coefficient 
evt  =  Antenna  Efficiency 
Gv  =  Gain 

PLFv  =  Polarization  Loss  Factor 
Avem  =  Maximum  Effective  Aperture 

hvem  =  Maximum  Effective  Height 
HORIZONTAL  CAGED  DIPOLE  OVER  EARTH: 

AFh(01,41)  =  Normalized  Array  Factor 

rh(01)  =  Horizontal  Plane  Wave  Reflection  Coefficient 

E0h  (01 )  =  0  Component  of  E-Field 

E0ht(01,4H  =  0  Component  of  Total  E- Field 

E0htn (01,41)  =  0  Component  of  Total  Normalized  E-Field 

E4h(01)  =  4  Component  of  E-Field 

E4ht (01,41)  =  4  Component  of  Total  E-Field 

E4htn(01'41>  =  ♦  Component  of  Total  Normalized  E-Field 

Uh(01,41)  =  Radiation  Intensity 

Phrad  =  Radiated  Power 
Dh0  =  Directivity 

EIRPh  =  Effective  Isotropic  Radiated  Power 
Rhr  =  Radiation  Resistance 
Rhin  =  Input  Resistance 

TH  =  Voltage  Reflection  Coefficient 
tht  =  Antenna  Efficiency 
Gh  =  Gain 

PLFh  =  Polarization  Loss  Factor 
Ahem  =  Maximum  Effective  Aperture 

hhem  =  Maximum  Effective  Height 
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The  following  known  or  estimated  data  must  be  entered: 

N  =  Number  of  Dipoles 

a  =  Radius  of  Antenna 

L  =  Length  of  Antenna 

h  =  Antenna  Height  Above  Ground 
f  =  Frequency  of  Interest 
IQ  =  Antenna  Feed  Current 
ZQ  =  Input  Impedance 

£r  =  Relative  Dielectric  Constant  of  Ground 
a  =  Conductivity  of  Ground 

ow  =  Incoming  Wave  E-Field  Unit  Vector  for  Antenna  (Free  Space) 
oa  =  Antenna  Unit  Polarization  Vector  (Free  Space) 
ovw  =  Incoming  Wave  E-Field  Unit  Vector  for  Vertical  Antenna 

0Vp  =  Direction  of  Incoming  Unit  Vector  (Vertical  Antenna) 

<jhw  =  Incoming  Wave  E-Field  Unit  Vector  for  Horizontal  Antenna 

0hp  =  Direction  of  Incoming  Unit  Vector  (Horizontal  Antenna) 

j  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 


Enter  input  data  here: 


I1' 

'l\ 

aw  := 

°i 

(dimensionless) 

oa  := 

0 

ioj 

A 

(dimensionless ) 


ovw 


(dimensionless)  0vp  = 


80-  a 

Tio 


(dimensionless ) 


ohw  :=  0 


(dimensionless)  0hp:=O  (dimensionless) 


N  =8 


(elements ) 


lo  =  10  (A) 


a  :=- 
8 


(meters) 


Zo  '.=75.0  (Q) 
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L  *0  5 


(meters ) 


h  •=  0.25 


(meters) 


j  :*  90 


o  ■=  102 


f  :=  3  10 


(Hz) 


fx  :=  15 


Define  constants  and  calculate  wavelength: 


c:=2  9979  10*  (meters/sec) 


T)Q  :=  1 20  h 


(meters/cycle) 


l  •  2-n 

k  :» — 


i  -  Io 

In  ■=  — 
N 


(A) 


rff  •=  105 


1  m'9 

EO  = -  10 

36  * 


( farads/m) 


ted  :=  1 


i  =  0  999 


(meters/cycle) 


z  =0  2  j 


Calculate  caged  dipole  parameters  parameters: 


Define  angular  offset  8 : 


0  =0,—  ..  2  ji 


( radians) 


4  =0,—  .  2  n 


( increments) 


( siemens/m) 


(dimensionless) 


< 


(Q) 


(m  x) 


(meters) 


(dimensionless ) 


( increments ) 


( radians) 
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Minimum  Distance  to  Far-Field  rm1n: 
r>  1.6-31  (meters) 


r,  =5  VI 


L24(2«)2 


(meters ) 


_  2  [l24-  (2  a)2] 


(meters) 


rmin  :=  max(r)  (meters) 


rmin  =  2.795 


VERTICAL  CAGED  DIPOLE  IN  FREE  SPACE: 


Array  Factor  AF(0): 

n-l.N  (elements)  <J>  :=2n  — 

"  N 


AF(«,*)  sJV* 


j  k»«in(9) co«(t-  <!>,) 


E-Field  of  finite  dipole  E(8): 


E(0)  :=j  T)o. 


2 nrff 


fcos 

fk  L 

—  cos(0) 
12 

-  cos 

m 

\  sin(0)  / 

(meters) 


( radians) 


(dimensionless) 


(V/m) 
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Et(8,*)  :»E(0)  AF(0,*) 


(V/m) 


Etn(0,4) 


(V/m) 


Radiation  Intensity  U  (0)  : 

U(M)  ■=— (|Et(0.*)|  )2 

2 tip 


(W/solid  angle) 


Radiated  Power 


f  2  * 


Prad  :* 


x 

U(0»-sin(0)d0d* 


0  '0 


(W) 


Prad  =0441 


(W) 


Directivity  D^: 


Umax  :=  U 


(W/solid  angle) 


Do  := 


4  x  Umax 
Prad 


(dimensionless) 


Do  =  1.543 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :=  Prad  Do  (W) 


EIRP  =0.681  (W) 
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Radiation  Resistance  Rr: 

Rr  :=  — (Q)  Rr  =  56  487  (Q) 

( |ln|  )2 

Input  Resistance  Rin: 

Rin  := — —  < Q ) 

2 

Voltage  Rejection  Coeffecient  T: 

r  Rin  -  2o  ,  . .  . 

1  .* -  (dimensionless) 

Rm  +  Zo 

Reflection  Efficiency  g«- : 

d  seed  [  1  -  ( |r|)2]  (dimensionless)  d=0  98  (dimensionless) 

Ga in  G : 

G:  =  et  Do  (dimensionless)  GdB  :=  10  log(G)  (dB) 

G  =  1.513  (dimensionless)  GdB  =  1.798  (dB) 

Polarization  Loss  Factor  PLF : 

PLF  *  (low  oal)  (dimensionless)  PLF  =  1  (dimensionless) 


Rm=  56.487  (Q) 


r  =  —0  1 4 1  (dimensionless) 
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Maximum  Effective  Aperture  Acm: 

Cl2 

Aem  :* - PLF  (m^)  Aem  =0  12  (m^l 

4  * 


Maximum  Effective  Height 


hem  :=  2  [—  (m)  hem  =0  268  (m) 

V  rp 

Define  Frequency  Relationship  F(x): 

F(x)  =- 5.696  10'4x+ 0.31481  (dimensionless) 


Frequency 

Range 

flow  and  fhighl 

flow  =  f-  f|- 

(Hz) 

flow  =2.056  10® 

(Hz) 

fhigh  :■  f-*-  f|- 

rjf 

(Hz) 

fhigh  =  3  944  10® 

(HZ) 

Bandwidth  BW : 

BW  =  fhigh  -  flow 

(Hz) 

BW  =  1  888-10® 

(Hz) 
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£ -FIELD  PATTERN:  This  far-field  radiation  pattern  is  for  the  vertical 
caged  dipole  in  free  space: 


Number  of  Dipole  Conductors : 


N  =8 


Radius  of  Antenna: 


a  =0  125 


(meters ) 


X(0).=  |Etn(0,*)|  cosjo-  — | 


Y(0)  •=  |Etn(0»|  sin  -  -j 


0  =  0° 


xref  ,X(8) 
m 


0  =  180° 
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VERTICAL  CASED  DIPOLE  OVER  EARTH: 

Define  angular  offset  81  and  82: 

( radians ) 


01  z.l,.±  +  JL  .1 


2  2  2  j  2 


«•  *  7 2 


92 


( radians ) 


2  2  2  j 

Relative  complex  effective  dielectric  constant  crpj 


erp  =er- j 


2  n  f  oo 


Array  Factor  AFv(0): 

Ej  tMimen 

Inc  v  * 

n 

AFvgp,  =£,„<'  ‘-N 


E-Field  of  finite  dipole  Ev(81): 


Ev(01 )  =  j  tjo 


I  -J  k  fff 
In  e  1 

^cosj 

fkl  \ 

- cos(01 ) 

\  2  y 

|  -  cosj 

m 

2-n-rff 

i  sin(  0 1 )  f 

.  -j  krff 

_  In  e 

Evgp  =j  tjo - — 

2  » rff 


(k  l  /k  Li 

cosl - cos^02  H  -  cosl — j 


sin^92j 


Vertical  plane  wave  reflection  coefficient  rv(81): 


rvc»n  :s«ES£g*12 


-  •!&?- 


sin(01 ) 


ap  cos(01 ) 


+•  -Jap-" 


sm(0I ) 


( radians) 


(dimensionless ) 


(dimensionless ) 


(dimensionless) 


(V/m) 


(V/m) 


(dimensionless ) 
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fvl  :=- 


cos^02^  -  Jerp-  sin  ^02  J  ‘ 
•  cos^02^  +  Jerp-  sm^02  y 


Total  E-Field  Ev t (01)  &  Normalized  E-Field  Ev t n (61  i 
Evt(ei,4)  =  Ev(01)  AFv(01,*)  (e>  kh“*<9,,4  rv(01 )  e‘J  khc“<8,)) 


Evtgp^ :« 


Evgp  AFvgp 


2  2 


|ej  k  hcc(92)  +  fvi  e-j  k-h-co.(fe}j 


Evtn(61,*) 

max(Evtgp) 


Radiation  Intensity  Uv(91_) 

rff2 

Uv(01 ,4)  :=  — — ( |  Evt( 01  ,f)|  )2 

2  ip 


Uvl  :=-i2l  (|Evtgp|) 
z  2  T|0 


Radiated  Power  PvraH: 

■2n 


Pvrad  := 


0 


Uv(01,*)sin(01)d01  d* 


0 


Pvrad  =  0  483 


Directivity  Dv^: 
Uvmax  :=  Uv(0vp,O) 


(dimensionless) 


(V/m) 


(V/m) 


(V/m) 


(W/solid  angle) 


(W/solid  angle) 


(W) 


(W) 


(W/solid  angle) 
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Directivity  Dv0: 
4  n  Uvmax 


Dvo 


Pvrad 


Dvo  =  ^415 


Effective  Isotropic  Radiated  Power  EIRPv: 

EIRPv  :» Pvrad  Dvo  (W)  EIRPv  =2  133 


Radiation  Resistance  Rvr : 

Rvr;_?Pvrad  (Q)  Rvr  =61  857 

(|ln|)2 


Input  Resistance  Rvln: 

Rvin  - — — —  (D)  Rvin  =  61  857 


Voltage  Rejection  Coeffecient  rv : 

fV  — —  (dimensionless)  TV  =-0  096 

Rvin  +•  Zo 

Reflection  Efficiency  cv,.  : 

cvt  :=ocd  [  1  -  ( |rv|  )2]  (dimensionless )  cvl=0  99l 


(dimensionless ) 


( dimensionless ) 


(W) 


(0) 


<Q) 


(dimensionless ) 


(dimensionless) 
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Gam  Gv: 


Gv  =EVtDvo  (dimensionless)  GvdB  =10  1og(Gv)  (dB) 


Gv  =4.374  (dimensionless)  GvdB  =6409  (dB) 


Polarization  Loss  Factor  PLFv: 
Evmax  :=  Evtn(5vp,0)  (V/m) 


ova 


Evmax 

lobesv 

0 

0 


(dimensionless) 


lobesv  :=  floor 


PLFv 


:  ( |ow  ova|) 


(dimensionless)  PLFv  =0  986  (dimensionless) 


Maximum  Effective  Aperture  Apm: 

Avem  :=  —  PLFv  ( m2 ) 

4n 


Maximum  Effective  Height  hvcnl: 


Avem  =0.343  (m2) 


f 


,  .  .  Rvr  Avem  (  , 

hvem  =  2  1 -  (m) 


hvem  =0  474  (m) 


TP 
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E-FIELD  PATTERN :  This  far-field  radiation  pattern  is  for  the  vertical 
caged  dipole  over  earth: 


Number 

of 

Dipole  Conductors: 

N  =8 

(elements ) 

Radius 

of 

Antenna : 

a  =0.125 

(meters) 

Height 

of 

Antenna  Above  Earth: 

h  =0.25 

(meters ) 

Xl(01 )  :=  |E'  tn(01  ,*)|  cosL  ■+ 


Y  1( 01 )  :=  |Evtn(01  ,|)|  sinfoi  4--) 


HORIZONTAL  CAGED  DIPOLE  OVER  EARTH: 


Define  angular  offset  $1  and  : 


41  :=0,-..2n 

j 


;  radians) 


Array  Factor  AFh  (91, 41 )  : 


n  =  1  N 


(elements ) 


9  '*  2  n  — 
"  N 


( radians) 


Ej  k  1 (nn(»l  )  sin<6l  yun(@h  4-co«<  01 )  eo»fai) 

In  e  v  ™  v  / 


(dimensionless) 


'z:*EIneJ 


ka(..n(*i)  *in  tin  ■+0.(62)  COC  &> 


(dimensionless) 


0  component  of  E-Field  of  finite  dipole  E8h(81) 


E0h(01 )  =j  rp- 


-I  k  rfr 
i  e  1 

2  7i  rff 


cosf - sin(01 )]  -  cosl 


F 1  -  sin(01 ) 


(V/m) 


E0hl  '=  j  1)0 - 

z  2  n  rff 


[I  -  sin^02^ 


(V/m) 


Total  E-Field  E8ht  (01, jl )  : 


E0ht(01,*l):=E0h(01)  AFh(01,*l)  (e>  khcos(9l)-  rv(01)  e'j  khc“<e,>)  (V/m) 


E0htl  :=  E0hl  AFhl 

z  z  z 


-  f  v  I  e 


-j  k  h  01(62)  J 


(V/m) 
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Normalized  E-Field  Etjh t n  (81 ,  41 )  : 

E6htn(91  ,41 )  :=  ^ht(ei'*!^ 
max(  E6hi  1 ) 


Horizontal  plane  wave  reflection  coefficient  Eh  (91) 


rh(0i ) cos(91 )-  Jerp- sin(01 ) 


cos(01 


)-nJap- 


sin(0l ) 


rhl 


(02  J  -  Jcrp-  sin  ^02  J " 
— .  = 

(02 J  ^  op-  sinfey 


4  component  of  E-Field  of  finite  dipole  Ejh(01): 

a 


In  e  J  krfr 

E^h(01 )  =j  rp - 1  I  -  cos 


2  ti  rff 


riL1  .  In  c  J  krff  /,  / k  L\\ 

z  2  7i  rff  l  A  2  }} 


Total  E-Field  Ejjhf  (01 ,  $1 )  &  Normalized  E-Field  E$hrn  (6 

E*ht(01  >1 )  :=  E*h(01 )  AFh(01 ,41 )  (e1  kl,  co‘<ei>+  fh(01 )  t>  khcos<ei)) 


E4htlz:= 


(j  khcotffe^ -j  khco»(()2^\  1 

e  W4rhlzeJ 


(V/m) 


(dimensionless) 

(dimensionless) 

(V/m) 

(V/m) 

>4D  : 

(V/m) 

(V/m) 


E4htn(81 ,41 )  jE4ht(0Ml_)| 
max(E4htl ) 


(V/m) 


Radiation  Intensity  Uh(81,jl): 

Uh(01  ,*l )  •=—  [( | E6hl( 9 1 ,41  )|  )24  ( i E4ht< Q 1  ,*l  )|  )2] 
2  T)0 


(W/solid  angle) 


(W/solid  angle) 


Radiated  Power  Pvrar<: 


2  * 


f! 


Phrad  - 


JO 


Uh(61  >1  )  sm(01 )  <301  d*l 

0 


(W) 


Phrad  =1059  (W) 

Directivity  Dh^: 

Uhmax  •=  Uh(0hp,O)  (W/solid  angle)  Dho  •=  4  *  (dimensionless) 

Phrad 


Dho  =  5  051 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRPh: 


EIRPh  :=  Pvrad  Dvo  ( W ) 


EIRPh  =2.133  (W) 


Radiation  Resistance  Rvr: 

_2_Phrad  (Q)  Rhr  =  135  58  (Q) 

(|ln|)2 
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Rhin  =  135  58 


Input  Resistance  Rhln: 
Rhin  := — — —  (Q) 


Voltage  Relection  Coeffecient  fH : 

. .  Rhin  -  Zo  ,  ,  . 

fH  •= -  (dimensionless) 

Rhm+  Zo 


(Q> 


rH=0  288  (dimensionless) 


Reflection  Efficiency  cvt : 

cht  :  =  8cd-[  I  -  ( |fH|  )2]  (dimensionless)  cht  =  0.917  (dimensionless) 


Gain  Gh: 

Gh  =ditDho  (dimensionless)  GhdB  =10log(Gh)  (dB) 

Gh  =4.633  (dimensionless)  GhdB  =6659  (dB) 


Polarization  Loss  Factor  PLFh: 
Ehmax:=E4htn(0hp,Oj  (V/m) 


aha  := 


0 

0 

Ehmax 


lobesh 


(dimensionless ) 


PLFh 


(dimensionless ) 


lobesh  :=  floor 


PLFh =  (dimensionless) 
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Maximum  Effective  Aperture  Ahain: 


Ahem  PLFh  (m2) 


Ahem  =  35.536  (m2) 


Maximum  Effective  Height  hham: 


.  .  RhrAhem  , 
hhem:*2  i -  (m) 


hhem  =7  15 
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E- PLANT  RADIATION  PATTERN :  This  far-field  radiation  pattern  in  the 
(x-z)  plane  is  for  the  horizontal  caged  dipole  over  earth: 


Number 

of 

Dipole  Conductors: 

Z 

II 

00 

(elements) 

Radius 

of 

Antenna : 

a  =0  125 

(meters ) 

Height 

of 

Antenna  Above  Earth: 

h  =0  25 

(meters ) 

-i  o  l 

xref  ,X3<  61) 
m 
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I  a 


H- PLANE  RADIATION  PATTERN :  This  far-field  radiation  pattern  m  the 
(y-z)  plane  is  for  the  horizontal  caged  dipole  over  earth: 


Number 

of 

Dipole  Conductors: 

Z 

M 

00 

(elements) 

Radius 

of 

Antenna : 

a  =0  125 

(meters) 

Height 

of 

Antenna  Above  Earth: 

h  =  0.25 

(meters ) 

X2(01 )  :■  |E*hln(61 ,0)|  cosfoi  +  Y2(01 )  >  |E*htn(01 ,0)|  sinjoi  +  ~ 


THE  HELICAL  ANTENNA  (REFLECTOR  OPTION) 

MATHCAD  SOFTWARE -DISH_HEL . MCD 

When  built  to  the  proper  specifications,  the  helical  antenna  possesses 
many  qualities  which  make  it  suitable  for  a  wide  variety  of 
communications  applications.  If  the  following  conditions  are  satisfied 
the  helix  will  exhibit  a  highly  directional  axial  mam  lobe,  low  side 
lobe  level,  negligible  mutual  interference  with  adjacent  antennas,  low 
voltage  standing  wave  ratio  (VSWR) ,  and  resistive  input  impedance  over  a 
wide  frequency  band: 


.8  <  CX  <  1.15 

n  >  3 

12  <  a  <  14 

(Note:  1  in  a  subscript  indicates  the  dimension  is  in  wavelengths. 

Mathcad  equations  can  not  use  symbolic  subscripts.  Therefore,  the  symbol 
1  will  immediately  follow  the  parameter  in  equations  (i.e..  Cl)  in  lieu 
of  subscripts.) 

The  helical  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  are  calculated  for 
both  axial  and  peripheral  feed  geometries) : 

C  =  Circumference  of  Helix 

1  =  Wavelength 

a  =  Pitch  Angle 

Do  =  Directivity 

p  =  Relative  Phase  Velocity 

Eg  ^  =  Electric  Field  Components 

¥  =  Array  Factor  Phase  Shift 

U  =  Radiation  Intensity 

Prad  =  Radiated  Power 

R  =  Antenna  Input  Resistance* 

T  =  Voltage  Reflection  Coefficient* 

Er  =  Reflection  Efficiency* 
hem  =  Maximum  Effective  Height* 

G  =  Gain* 

EIRP  =  Effective  Isotropic  Radiated  Power* 

Aem  =  Maximum  Effective  Aperture* 

AR  =  Axial  Ratio 

PLF  =  Polarization  Loss  Factor 

BW  =  Bandwidth 

f high  =  Upper  Frequency  Limit 
flow  =  Lower  Frequency  Limit 
Acceptable  Conductor  Diameter 

Ex, yt z  =  Electric  Field  Cartesian  Coordinates 
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6p'4p  =  Unit  Polarization  Vector  Coordinate  Angles 

oa  =  Antenna  Unit  Polarization  Vector 
rrain  =  Mininum  Distance  to  the  Far-Field 


The  following  data  must  be  input  based  on  known  or  estimated  data: 

D  =  Diameter  of  Helix  (Center  to  Center) 

S  =  Spacing  Between  Turns  (Center  to  Center) 

L  =  Length  Along  Conductor  of  One  Turn 
n  =  Number  of  Turns 
d  =  Diameter  of  Helical  Conductor 
f  =  Frequency  of  Interest 
m  =  Desired  Mode 

i  =  Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 

IQ  =  Antenna  Feed  Current 

ZQ  =  Characteristic  Feed  Impedance 

ow  =  Incoming  Wave  Electric  Field  Unit  Vector 


J5 

* 

0 


Enter  input  data  here: 


(dimensionless ) 


D 


J! 


(meters ) 


L  :=  0.097  (meters) 


n  :=  10 


(turns ) 


f  '.=  3  1010  (Hz) 


m  '=  1 


(dimensionless) 


Z0=150  (Q) 


169 


S  =  0  023 

(meters ) 

i  =360  (degrees) 

d  =  0.005 

(meters ) 

lo  =  1  (A) 

x  :=  1 

(meters ) 

z  :*  1000  (meters) 

y  =1 

(meters) 

Calculate 

helical  qeometnc 

parameters  and  define  constants: 

c  :=  2  9979  10® 

(meters/sec) 

:=  1 20  Jt  (Q) 

C  ::»  D 

(meters ) 

,  c  , 

X=-  (meters/cycle) 

f 

C  =0  1 

(meters) 

X=9993‘I03  (meters/cycle) 

cx^ 

1 

(dimensionless ) 

a:=atanf— ]  (radians) 

\c ) 

Cl  =  10  007 

(dimensionless ) 

a  =  0  22607  (radians) 

1 

(dimensionless ) 

s 

SX  :=—  (dimensionless) 

1 

LI  =  9  70679 

(dimensionless ) 

SX  =  2.30161  (dimensionless) 

..  180 
ad  •= - a 

K 

(degrees) 

k  :=  2  -  (m_1) 

X 

ad  =  12.95276 

(degrees) 
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Minimum  Distance  to  the  Far-Field 


Calculate  helical  antenna  parameters: 


Define  angular  offset  8 
from  helical  axis:  £mmi 


0:SO,~..2  * 

i 


( radians) 


*:-0, 


2  7t 


( radians ) 


r0=!6i 

Tj  :=  5  n  S 

2  (n  S)2 
Fj  '  1 

rmin  :=  max(r) 


rmin  =  10  58741 


Relative  Phase  Velocity  p: 


P  :s 


Li 


S1+  nu 


p  =2.89616 


Array  Factor  Phase  Shift  yr 


Hi(8)  :=  2  7i 


|si  cos(0)- 


(meters) 


(meters) 


imeters ) 


(meters) 


(meters) 


(dimensionless ) 


(dimensionless) 


( radians ) 
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Electric  Field  Components  Efl,E^: 


cos(0) 


(V/m) 


E0(0)  :=  E(0) 


(V/m) 


E*(0)  :=j  E(0)  (V/m) 


Radiation  Intensity  U(0): 

U(0)  •»—■(  |E(0)|  )2  (W/solid  angle) 

% 


Radiated  Pov.er  PraH: 


f2« 


Prad  = 


n 

U(0)  sin(0)d0d* 


0  -'0 


(W) 


Prad  =  1  67351' 10 


(W) 


Directivity  D^,: 

Do  *  12  C12  n  S3L  (dimensionless!  Do2  '=  -  (dimensionless) 

Prad 


Do  =2. 7658  104  (dimensionless)  Do2  =  19.91822  (dimensionless) 
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Axial  Ratio  AR: 


AR  :* 


LX 


jsin(a)  - 


i) 


AR  =  1  17586 


Effective  Isotropic  Radiated  Power  EIRP: 


ElRP^PradDo  (W) 


EIRP2  :  =  Prad  Do2 


EIRP  =46.286  (W) 


EIRP2  =0  03333 


Polarization  Loss  Factor  PLF: 


9p 


( radians) 


0p  =  1.41421  10"J 


( radians ) 


♦P  =  atan|— | 


4p  =0.7854 


Ex  =  E0(0p)  cos(0p)  cos(4p)-  E4(0p)  sin(4p) 


Ey  =E9(0p)  cos(0p)  sin(4p)+  E4(0p)  cos(4p) 


Ez  =E0(0p)  sin(Op)  -  1 


i 

IH 

ca  = 

lEy 

\Ez/ 

/J(|Ex|)24-(|Ey|)2+(|Ez|)2 

O 

L/3 

1 

o 

oa  = 

0  5  +  0.5j 

-9  99999*  10_J 

(dimensionless) 


( dimensionless ) 


<W) 


(W) 


( radians) 


( radians ) 


(V/m) 


(V/m) 


(V/m) 


(dimensionless ) 


(dimensionless) 
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Polarization  Loss  Factor  PLF : 


PLF  ■ 


(dimensionless) 


PLF  =  1 


Radiation  Resistance  Rr: 

Rr  :=  2  — (ft)  Rr  =  3.34702  10'3 

(Uo|)2 


Dual  Parameters 


Axial  Feed 


Peripheral  Feed 


Input  Resistance  R: 
Ra  :=  140  «Jcji  (ft) 

Ra  =  4.42874  102  (ft) 


Rp  =  47  41756 


Voltage  Reflection  Coef fecientj : 
Ra-  Z„ 


fa 


Ra-f  Z„ 


Ta=  0  49399 


(dimensionless) 


(dimensionless) 


r  _RP'Zo 

Tp  .= - 

Rp+Z0 

fp  =-0.51962 


Reflection  Efficiency  gr: 

Era  :=  1  -  ( |ra|  )2  (dimensionless)  crp  :=  1  -  (  |rp|  )2 

tra  =  0.75597  (dimensionless)  erp  =0.72999 


(dimensionless) 


(ft) 


(ft) 


(ft) 


(dimensionless ) 


(dimensionless) 


(dimensionless ) 


(dimensionless) 
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Gain  G 


Ga  :=  era  Do 

(dimensionless ) 

Gp  >  erp  Do 

(dimensionless ) 

Gadb  : s  1 0  log( sra  Do) 

<dB> 

Gpdb  :=  10  log(crp  Do) 

(dB) 

Ga  = 2 09088* 104 

(dimensionless ) 

Gp  =2  01902*  104 

(dimensionless) 

Gadb  =43.20328 

(dB) 

Gpdb  =430514 

(dB) 

Ga2  EraDo2 

(dimensionless ) 

Gp2  :  =  erp  Do2 

(dimensionless ) 

Gadb2  :=  10  log(era  Do2) 

(dB) 

Gpdb2  :=  10  log( Erp  Do2 )  ( dB ) 

Ga2  =  15  05765 

(dimensionless) 

Gp2  =  14  54017 

(dimensionless) 

Gadb2  =  1 1  77757 

(dB) 

Gpdb2  =  11.62569 

(dB) 

Maximum  Effective 

Aperture  A&m: 

.  Era  i2  Do 

Aema  -= -  PLh 

(m2) 

Aemp  :=£P  A-Do  plf 

(m2 ) 

4  it 

4  71 

Aema  =0.16615 

(m2 ) 

Aemp  =  0  1 6044 

(m2 ) 

A^:.ml2D02PLF 

(m2 ) 

.  erp  X.2  Do2  D1  c 

Aemp2  :=— - PLF 

(m2) 

4-71 

4-71 

Aema2  =  11 9657  1 0  4 

(m2 ) 

Aemp2  =  1.15545  10  4 

(m2 ) 
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Maximum  Effective  Height  h^m : 


hema:i2  JRr (m) 


hemp  :=  2  Rr. 


hema  =  2  42912*  10~3 


hemp  =  2.38701*  10 


hema2  :=  2  |Rr  Ae.ma^  (m) 


hemp2:=2  IRr 


Aamp2 


hema2  =6  51873  10 


hemp2  =  6.40574  10 


Bandwidth : 


thigh  :  =  - 


flow 


thigh  =3  44759  10* 


flow  =2.39832*  10* 


BW  :=  thigh  -  flow 


BW=  1  04927  10* 


Acceptable  Conductor  Diameter: 


dmin  :=  005  1 


dmax  05  X 


dmin  =4.9965  10 


dmax  =4  9965  10 
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PARABOLIC  RgrLECTOR  ANALYSIS  WITH  HELICAL  ANTENNA  FEED 

The  helical  antenna  Mathcad  application  (Reflector  Option)  will  compute 
the  following  additional  parameters  (Note:  a  2  or  22  will  follow  most 
secondary  field  parameters  in  this  application,  and  inorder  to  speed  up 
computation  time  corresponding  vectors  are  created  with  subscript  A  &  B 

px, y  =  E-Field  Integral  and  Summation  Equation  Vector  Components 
Eax, ay  =  Aperture  E-Field  Components 

0’,p  =  E-Field  Integral  and  Summation  Equation  Parameters 

Eg  =  Secondary  Field  Theta  E-Field  Component 

E^  =  Secondary  Field  Phi  E-Field  Component 

U  =  Secondary  Field  Radiation  Intensity 

Pra<3  =  Secondary  Field  Radiated  Power 

D0  =  Secondary  Field  Directivity 

G  =  Secondary  Field  Gain 

eap  =  Aperture  Efficiency 

£s  =  Spillover  Efficiency 

£t  =  Taper  Efficiency 

£p  =  Phase  Efficiency 

£x  =  Polarization  Efficiency 

£ohmic  =  Ohmic  Efficiency 

Ae  =  Reflector  Effective  Area 

he  =  Reflector  Effective  Height 

EIRP  =  Reflector  Effective  Radiated  Power 

Ex,y,z  =  Secondary  Field  Cartesian  Components 

oa  =  Reflector  Unit  Polarization  Vector 

PLF  =  Reflector  Polarization  Loss  Factor 

rmin  =  Minimum  Distance  to  Secondary  Far-Field 

0O  =  Reflector  Periphery  Coaltitude 

9p'^p  =  Theta  and  Phi  for  Polarization  Loss 


The  following  additional  data  must  be  inputted  based  on  known  or 
estimated  data: 

EQ  =  Electric  Field  Scale  Factor 

aref  =  Radius  of  the  Mouth  of  the  Reflector 

foe  =  Reflector  Focal  Length 

Ek  =  Blockage  Efficiency 

£Sp  =  Spar  Efficiency 

tl,t2  =  Radiated  Power  Increments  for  Phi  and  Theta 
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il  =  Secondary  Field  Increments 
N,M  =  Summation  Increments 

rff2  =  Secondary  Field  Observation  Distance 
x,  y,  z  =  Polarization  Loss  Factor  Coordinates 
ow  =  Incoming  Wave  Unit  Polarization  Vector 


Enter  additional  input  data  here: 

foe  :=  5  (meters)  aref  :*  5 

(Not*:  (f)  should  b*  greater 
than  minimum  distance  to  the 
far-field  for  the  feed  antenna) 


cb  :  .96 


esp  =  .95  (dimensionless)  tl  =  10 


t2  =  10 


(increments) 


il  =10 


N  =  10 


(increments) 


M  =  10 


rflE2  r  103  (meters ) 


x  =  1000 


ow2  := 


±_ 

/2 

_1_ 

* 

0 


(dimensionless) 


y  =  1000 


z  =1000 


(meters) 


(dimensionless) 


(increments) 


(increments) 


( increments ) 


(m) 


(m) 


(m) 
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Establish  Integral  and  Summation  Increments: 

^1  :*0,^-2..2  h  (radians)  A  •'=  I  N 

tl 


+2  2-x 

il 


( radians ) 


B  :=  1  M 


02  ur6  — - I-  icr6  —  -  io‘‘ 

\36/  il  18  36  36 


ill  *2-*— 
B  M 


rp  :sO,— ..  aref  (m) 

t2 


*22®  *2"m 


rpl  •=  aref — 
A  N 


(m) 


022. 


:=  (.I]  JL 
ll§/  N  36 


022  :  =  -  — 
0  36 


Calculate  reflector  geometric  parameters: 


/  aref1 
0o  .=  2  atan  — — 
,2  focj 


( radians ) 


00  =  0.9273 


( radians ) 


4  foc2+ rp2 

P(rp)-= - 77 -  <m> 

4  foe 


0p(rp)  :=  2  atan ( -*£— 
\2  foe, 

4  foc2+  (rpl^ 


pl 


4  foe 


( radians ) 


(m) 


f— ^ 

U  foe/ 


0pl  :=2atan -  (radians) 

A  12 


r0:=16  1 


r(  :=  10  aref 


rmin  :=  max(r) 


rmin  =  2  0014  104 


(increments) 


(increments ) 


(radians) 


( radians) 


( radians) 


( radians ) 


(m) 


(m) 


(m) 


(m) 


(m) 
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(V/m) 


Aperture  Electric  Fields  Ea _ Ea yj_ 

Eax(rp,41 )  :=  (E0(0p(rp))  -cos(41  )-►  E4(0p(rp))  sin(41  ))-e*  k  < P<n>>- foc 2> 


Eay(rp,41 )  =  (E0(0p(rp))  sin(41  )■»■  E4{0p(rp))  cos(4l ))  e1  k(p<,ll>>  foc2)  (V/m) 


EaxlA  B:  =  [E6(eplJ  (-cos(*l  lj)  +  E^GplJ  s.n(*I  1 J  ]  eJ  + 

EaylA.B  :=  (E0^1/)  Sm(*'  0  +  E*(0p1/)  COS(*'  !B ))  C>  ^  ^'A 


(V/m) 


(v/m) 


Electric  Field  Integral  Equation  Vector  Components  Px,yl 

Iaref 

Eax(rp,41 )  e>  “»< ♦' >  ♦2>^‘in< ♦'»  s,n<*2»  rp drp 

0 

1 2  Jt 

Px(02,42):=  Ax(02,*2>l)d<M 

Jo 

'aref 

Ay( 02, 42,41 )  •-  Eay(rp,41  )  e*  k^n(e2,.(c«»(4.».c«(02)  +  .in(*i)..,n<*2)).ipdrp 

Jo 


Py(02,42):=  Ay(02,42,4Dd41 

Jo 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


AxlA  b  :=  EaxlA  B  ej  -N)  .,pl  ^  (V/m) 


AxlA_Bd4N 


Ayl.n:=  EaylARe 


j  k  rpl  sin^022/)  (cos|fl  IpJ  COS  (♦22^  +  sin 1 1 q)  «n(*22g)) 


(V/m) 


rpl  drp  1  (V/m) 


AyiA.Bd4>1 


(V/m) 
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Secondary  Field  Electric  Field  Components  Ep,  E^: 
i  k  e’J  k  rin 

E02(82,*2)  =2 - (Px(02,*2)  cos(*2)+  Py(02,*2)  sin(*2)) 

2  *  rfi2 


i  ke'J  krfn 

E*2(02,42)  -l - (Px(02,*2)  -sm(*2)  +  Py(02,*2)  cos(*2))  cos(02) 

2  n  rff2 


.  -j  k  rfTl 


.  -j  kriR 

m'-s  " i-*m  lPxl«»(' 4 ^'a.b <H*22J ] HK2J 


Secondary  Field  Radiation  Intensity  U: 

U2(02,*2):=—  [( |E02(02  >2)|  )2+  ( |E42(02.*2)|  )2] 
2 


U3*.B*^[(|E83A.B|)!*(|E*2a.B|)1 


Secondary  Field  Radiated  Power  Prarl: 

PnK!2  5  (ttI  (ts)  Z|LU2K-*22J  -K)) 

Prad2  =0.19318 


Secondary  Field  Directivity  D^: 

Do22  :=  4  n  max^U3) 

Prad2 


Do22  =  1.04762  105 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(W/solid  ang) 


(W/solid  ang) 


(W) 


(W) 


(dimensionless) 


(dimensionless) 
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Spillover  Efficiency  eg: 
■0o 

U(0p)  sin(0p)  d0p 

•0 

es* - 

•n 

U(0p)  sin(0p)  d0p 

•  0 


es=0  19013 


(dimensionless ) 


(dimensionless ) 


Taper  Efficiency  zt  '■ 


8o 

1 


U(fiE!l!  tan  M  dep 

Prad  1 2/ 


12 


2  cot  — 


reo 


U(6p)  4  n 
Prad 


sin(0p)  d0p 


(dimensionless ) 


et  =0  50837 


(dimensionless) 


Phase  Efficiency  c^: 

■  41  6  I0~  5  arcf 

Ep:=e  *  (dimensionless) 


ep  =0  68574 


(dimensionless ) 


Aperture  Efficiency 

ex  :=  98  (dimensionless)  cohmic  :=  98  (dimensionless) 

eap —es  cl  £p  EX-cb  ssp  Eohmic  (dimensionless)  tap  =0.05805  (dimensionless) 
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Effective  Isotropic  Radiated  Power  EIRP: 


EIRP22  :*  Prad2  Do22 


EIRP22  =2  02379  I  O’ 


Gain  G: 


G22  :  =  £ap-Do22  (dimensionless) 


G22dB  :=  10  log(G22)  (dimensionless) 


G22  =608192  10J  (dB) 


G22dB  =37  84041  (dB) 


Effective  Aperture  AP: 


Ae22  :=  eap  n  aref2 


Ae22  =4.55961 


Radiation  Resistance  R, 


Rr22  -2 


(|lo|)2 


Rr22  =0  38636 


Effective  Height  hf 


he22  :=  2  Rr22 


he22  =0.13672 
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Polarization  Loss  Factor  PLF : 

♦p2  :=  atan|^J 
+p2  =0  7854 

E*4  :=  E*2(0p2,*p2) 

Ex2  :=  E04  cos(0p2)  cos($p2)  -  E$4  sin(+p2) 

Ey2  :=  E64  cos(0p2)  sin(^p2)  -  E$4  cos($p2) 

Ez2  :=  •  E04  sin(  0p2 ) 

/Ex2 

oa2  :a- — . . .  * . - . -  |Ev2 

J(|Ex2|)2+(|Ey2|  )2-k|Ez2|)2  \Ez2, 


0p2  :*atanj 


(¥1 


[  radians ) 


0p2  =0  95532 


I  radians) 


E04  :=  E62(0p2,4p2)  (V/m) 


oa2  = 


0.12034  +  0  48533j 
0  12034  +  0  48533j 
-0  60556  -  0.36501  j 


PLF2 


:s  (|ow2  oa2|) 


PLF2  =0.25003 


( radians ) 


(radians) 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


(dimensionless ) 


(dimensionless ) 


( dimensionless ) 


(dimensionless ) 
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PARABOLIC  REFLECTOR  WITH  HELICAL  FEED  FAR-FIELD  ELEVATION  PATTERNS 


E-FIELD  RADIATION  PATTERN: 


ll  =20 


—  -  id6  _1.  io'6 

36  36 


X2(  02 )  = 


cos|  02  -f  — 

*> 


Y2(  62 ) 


E02 


-0  00.- -0  002  -0  002  -0  001  “0  001  0  0  001  0  001  0  002  0  002  0  00? 

X2<«2> 


E  FIELD 
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H-FIELD  RADIATION  PATTERN: 


X3i02  i  :=  |E*2t 02.0  »|  cosfu2  +  -| 


V’i  02 


=  | LI©2(  02 .0  i| 


-0  00? -0  002  -0  002  -0  001  -0  001  0  0  001  0  001  0  002  0  002  0  00? 


—  E  FIELD 


X?(  02 1 
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THE  S  F I  RAT  ANTENNA  'FEFLECTCF  OPTION 
MATH  CAT'  SOFTWARE -CT  SH_S FI  .  MOD 

Spiral  antennas  are  a  family  of  two  or  three  dimensional  devices  that 
possess  frequency  mdependant  parameters  over  a  wide  bandwidth.  Spiral 
antennas  are  commonly  used  for  direction  finding,  satellite  tracking 
and  missile  guidance. 

The  planar  spiral  may  be  of  the  Archemedean,  log-spiral,  or  equiangular 
type.  All  three  radiate  two  main,  circularly  polarized  lobes 
perpendicular  to  the  plane  of  the  antenna.  Additional  gam  for  planar 
spirals  may  be  achieved  by  placing  a  metal  cavity  on  the  side  of  the 
antenna  with  the  unwanted  lobe.  The  cavity  may  be  empty  or  be  filled 
with  electromagnetic  energy  absorbing  material.  These  applications 
pnnciply  examine  the  equiangular  planar  spiral  and  do  not  account  for 
cavity  backed  effects. 

(Note:  Mathcad  equations  cannot  use  symbolic  subscripts.  Therefore, 

symbols  like  /.  will  immediately  follow  the  parameter  m  equations.  ) 

The  spiral  antenna  Mathcad  applications  will  compute  the  following 
parameters  for  equiangular  planar  spirals: 

k  =  Wavenumber 
/.  =  Wavelength 

D0  =  Directivity 

E4  =  Electric  Field  Component 

U  =  Radiation  Intensity 

Umax  =  Maximum  Radiation  Intensity 

prad  =  Radiated  Power 
G  =  Gain 

EIRP  =  Effective  Isotropic  Radiated  Power 
Aem  =  Maximum  Effective  Aperture 
BW  =  Bandwidth 

rnun  =  Mlninium  Distance  to  Far-Field 

Rr  =  Radiation  Resistance 

hem  =  Maximum  Effective  Height 

fhigh  =  Upper  Operating  Frequency 

flow  =  Lower  Operating  Frequency 

rn  =  Any  point  on  the  nth  edge  of  a  spiral 

tex  =  Equiangular  Planar  Spiral  Expansion  Ratio 

Zi  =  Planar  Spiral  Input  Impedance 

T  =  Voltage  Reflection  Coefficient 

crv  =  Reflection  Efficiency 


L  8  8 


fir  =  rclat;:a:;cr.  Less  Fc,:::: 

'•high  =  Upper  Operating  Wavelength 

'•low  =  Lower  Operating  Wavelength 

A  =  Planar  Spiral  Electric  Field  Amplitude 


The  followmo  data  must  be  incut  based  on  known  or  estimated  data: 


M  =  Mode 

N  =  Number  of  Spiral  Arms 
f  =  Frequency  of  Interest 

1  =  Number  of  Increments  for  Far  Field  Radiation  Patterns 
Zff  =  Distance  of  Far-Field  Calculations 
IQ  =  Input  Current  at  Antenna  Terminals 
r0  =  Spiral  Feed  Point 
a  =  Flare  Rate 

6n+i  =  Angular  Arm  Width  of  ntn  Spiral  Arm 
♦r  =  Azimuth  to  Compute  Expansion  Ratio 
R  =  Overall  Radius 

E0  =  Source  Strength  Constant  for  Planar  Spirals 

ow  =  Wave  Unit  Polarization  Vector 

ZQ  =  Characteristic  Impedance  of  Feed  Assembly 

oa  =  Equiangular  Planar  Spiral  Unit  Polarization  Vector 


THE  PLANAR  SPIRAL  ANTENNAS 


Enter  input  data  here: 


N  =2 

( a  rms ) 

f  =  3  5  ID10 

(Hz) 

M  =  1 

(mode } 

lo  :*  1 

(amps ) 

Note:  Mmax 

is  N-l) 

o 

II 

( increments ) 

rfl".=  1  I05 

(meters ) 

OJ 

II 

( dimensionless ) 

ro  =  001 

(m) 

R  =  01 

(m) 

Eo  :=  103 

(V/m) 
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u 


I 


1 

-> 


6  = 


( radians ) 


cm  :*  1 


( dimensionless ) 


3  n 


Zo  100  (O' 


Calculate  planar  spiral  antenna  geometric  parameters  and  define 
constants : 


c  =2  9979  10*  (meters  /  sec) 


H0  :=  120  n 


lQ> 


c 

7 


(meters  /  cycle) 


/.  =  8  56543"  10  ?  (meters  /  cycle) 


s 

o 


(dimensionless ) 


(Farads  /  mu 


»i„  :=4  *  10 


.  7 


(H  /  mi 


Calculate  planar  spiral  antenna  parameters  : 

Define  angular  offset  6 
from  y-z  axis: 

0  :: .  1  _  10'6,-  -  —  -  10'6  4-  I0”6  (radians)  $:=0  ,-2-il  2  n  (radians) 

2  i  2  2  i 
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Distance  to  Far-Fieic 


Hq  ' ' 6  X 


n  j  =  H*  R  i  m  i 


^‘7 


rmin  :smax(rr) 


rmin  =  0  1 


Wavenumber  k: 


Radial  Distance  to  n±“  Spiral  Edge  r: 


»  (»-  Sn-  |) 

rin.fr)  =roe  (m) 


k=  7.33552  10*  (m'1] 


rU  .  2  7i )  =  4  U09 17*  10  ( m ) 


Expansion  Ratio  : 


,  i  .  _  r(n.ir-t-  2  it) 
eex(n,4r)  - - 

r(  n,  $r) 


(dimensionless ) 


eex(  1,2  it)  =4.00917  (dimensionless) 


Bandwidth  BW: 


Equiangular  Spiral 


Log-Periodic  Spiral 


Ihigh  =4  ro 


/.high  I  :=  20 to 


fhigh  :*■ 


fhighl  :  =  • 


/.high  1 
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lhigh  =  ~  4v47 5  |o 


-  Hz 


lhiiili  =  ”  4^)4' 5 


/.low  :=4  R 


( m 


flow  » 


How 


(Hz) 


flow  =7  49475  10 


(HZ  ! 


BW  •=  (high  -  flow 


(Hz)  BW|  fliighl  -  flow 


BW  =6  74528  10 


10 


(Hz,  BW1  =7  49475*10 


Electric  Field  E  (0J)  and  Electric  Field  Amplitude  A ( B ) 


w  =0  i  (increments) 


Eo  k' cos(0)  (  I  4  j  acos(0))  1  tan[-  e 


'  H*fH 


E(0)  :*■ 


sin(0}‘  rff 


A(0)  = 


M  I— I  )Unt  a  c os<y,l 

cos(0)tan|-j  e’ *' 


'  1  4 


sin(0)^l4a  cos(0) 


0)  :*  —  ln(l  4  a2  cos(0)2)  +  atan(a  cos(0) ) 
2  a 


E4(0)  :=  A(0)  e*J 


(Hz) 


(Hz) 


( V/m) 


(V/m) 


(dimensionless ) 


(V/m) 
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(V/m) 


Radiation  Intensity  U (6)  : 

U(0)  :■ — ! —  (A(  0)2) 

2% 


Umax  :=  max(  U I ) 


Umax  =2.2285-10 


Radiated  Power  PraH: 


Prad  :=  4  a 


U(9)  sin(0)  d6 


Prad  =7.94528-10 


-3 


Directivity  Dn: 

^  4  a  Umax  ,  .  . 

D  = -  (dimensionless) 

0  Prad 


DQ  =  3  52462 


(dimensionless) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W  /  solid  angle) 


(W) 


(W) 
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Radiation  Fes.srar.ee  F  r  : 


Rr^lESL  ,u , 


ir.Deaar.ee  ^ i 


(|lo|)2 


Zi  = 


N  30  7. 


Rr  =0  01589 


(fi) 


Zi  =  1  88496* 10*  (Q? 


Voltage  Reflection  Coefficient  T:  Reflecton  Efficiency  crv.: 

cn  =  1  -  (  Irl  ')*  i dimensionless  i 


Zi  -  Zo 

1  ;= -  (dimensionless 


Zi  +  Zo 


r  =  0  306"’5  (dimensionless1 


Ga in  G : 


cn  =0  90591  (dimensionless. 


G:=cr\D  (dimensionless)  GdB  :=  1 0  log(  G )  (dB) 


G =3  19298  (dimensionless) 


GdB  =  5  04 196 


(dB) 


Effective  Isotropic  Radiated  Power  (EIRP) 
Prad  G 


EIRP  =- 


sxv 


EIRP  =0  028 


(W) 


(W) 


Polarization  Loss  Factor  (PLF): 


.=  (UI)2 


PLF 


PLF  =  I 


(dimensionless ) 


(dimensionless ) 
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Kaxirr.urr  Effective  A 


D0  en  PLI 

Aem  * - 

4  n 


Aem  =  1  86416  I0’! 


Maximum  Effective  He 


h«n  *  lSL*2: 

J  % 

hem  =  5  6063*  10  ? 


cert  u  ;  e 


!  rr.- 


im*- ) 


ight  jhemiJ 


(IT.) 


(m) 
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THE  EQUIANGULAR  PLANAR  SPIRAI.  ANTENNA  TAR-FIELD  ELEVATION  PATTERN 


For  the  purpose  of  this  far-fieid  radiation  pattern,  the  spiral  antenna 
lies  parallel  to  the  Ey  =  0  grid  line  and  is  centered  at  the  origin. 

The  magitude  of  the  electric  field  pattern  is  rotationally  symmetric 
with  respect  to  the  Ex=0  grid  line.  The  equiangular  planar  spiral 
antenna  possesses  a  mirror  image  radiator,  pattern  m  the  -y  half  plane. 


E*x(0) : s  | E+( 0 ) | 


E#\(0)  *  |E*(0)|  sin|o-fL 


E«v(0)  uooi 


“0  001  0  0  0  0  “2  711  10  0  0  0  0  0  001 


—  E  FIELD 
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PARABOLIC  REFLECTOR  ANALYSIS  WITH  SPIRAL  ANTENNA  FEED 


The  spiral  antenna  *  ithcad  application  ‘Ref.ector  Option.'  will  compute 
the  following  ad"  -ucnal  parameters  (Note:  a  2  or  22  will  follow  most 
secondary  fiel"1  parameters  ir,  this  application,  and  inoraer  tc  speed  up 
computation  time  cc r responding  vectors  are  created  with  subscript  A  &  B : 

Px  y  =  E-Field  Integral  and  Summation  Equation  Vector  Compc.ner.ts 
Ea>;  ay  =  Aperture  E-Field  Components 

O’.p  =  E-Field  Integral  and  Summation  Equation  Parameters 
E((  =  Secondary  Field  Theta  E-Field  Component 

E^  =  Secondary  Field  Phi  E-Fieid  Component 
U  =  Secondary  Field  Radiation  Intensity 
Prad  =  Secondary  Field  Radiated  Power 
Dc  =  Secondary  Field  Directivity 
G  =  Secondary  Field  Gam 
cap  =  Aperture  Efficienci¬ 
es  =  Spillover  Efficienci¬ 
es  =  Taper  Efficiency 
Cp  =  Phase  Efficiency 
cx  =  Polarization  Efficiency 

cohmic  =  ohm-c  Efficiency 

Ae  =  Reflector  Effective  Area 

he  =  Reflector  Effective  Height 

EIRP  =  Reflector  Effective  Radiated  Power 

Ex,y,z  =  Secondary  Field  Cartesian  Components 

oa  =  Reflector  'Jmt  Polarization  Vector 

PLF  =  Reflector  Polarization  Loss  Factor 

rmm  =  Minimum  Distance  to  Secondary  Far-Field 

60  =  Reflector  Periphery  Coaltitude 

Op,*p  =  Theta  and  Phi  for  Polarization  Loss 


The  following  additional  data  must  be  inputted  based  on  known  or 
estimated  data : 

aret  =  Radius  of  the  Mouth  of  the  Reflector 
foe  =  Reflector  Focal  Length 
Cfc  =  Blockage  Efficiency 

£sp  =  Spar  Efficiency 

tl,t2  =  Radiated  Power  Increments  for  Phi  and  Theta 
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ll  =  Secondary  Field  Increments 
N,M  =  Summation  Increments 

rf£2  *  Secondary  Field  Observation  Distance 
x,  y,  z  *  Polarization  Loss  Factor  Coordinates 
ow  *  Incoming  Wave  Unit  Polarization  Vector 


Enter  additional  input  data  here: 


■ref  *0.32  (meters) 


foe  *0  147  (meters) 


cb  '■*  96 


esp  :*  95  (dimensionless) 


tl  :«  10 


12  :«  10 


( increments ) 


il  :■  10 


N  :«  10 


(increments ) 


M  :«  10 


rff2  :■  103  (meters ) 


x  :*  1000 


ow2  := 


_J_ 

_1_ 

h 

0 


(dimensionless ) 


y  :*  1000 


z  :=  1000 


(dimensionless ) 


( increments ) 


(increments ) 


( increments ) 


(m) 


(m) 


(m) 
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Establish  Integral  and  Summation  Increments: 

*)  « 

41  :*0~—  .2  ii  (radians)  A=1  N 

tl 

^2  :*  Iff3,— -►  Iff  3  2  !t+  Iff3  (radians)  B  =  1  M 

il 


(increments) 


(increments ) 


02  |—j-  Iff3,— 2 - 1-  Iff3  —  -  Iff 

36/  il  18  36  36 


aref 

rp:=0, -  aref  (m) 

l2 


g 

41  L  :«  2  a—  (radians) 

B  M 


( radians ) 


rpl  :=  aref — 
A  N 


(m) 


022  :*  I  —  I  — — —  (radians) 
A  V 18/  N  36 


022  ;*-  — 

0  36 


( radians ) 


Calculate  reflector  geometric  parameters: 


0o  :*  2  atari/-— 

\2  foe/ 

( radians ) 

r0:=1.6>. 

(m) 

0o=  1  65544 

( radians ) 

Tj  :  =  10  aref 

(m) 

4foc2+rp2 
P(,P)  '  4  foe 

(m) 

8  aref2 
r=;'  » 

(m) 

0p(rp)  :=  2  atanf— ] 
\2  foe/ 

(radians) 

rmin  :=  max(  r) 

(m) 

pl  .. <•«*’♦  (■P'J2 
A  4  foe 

(m) 

rmin  =  95.64028 

(m) 

(rp\A\ 

0pl .  *2  atan - 

A  \2  foe/ 

( radians ) 
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Aperture  Electric  Fields  EJV, _ Ea  vj_ 


Eax(rp,41 )  :«E*0p(rp))  sm(*I  )  e1  k'u',rp,-fot'' 


(V/m) 


Eay(rp,*l )  '«  E#0p(rp))  cos(*l )  e>  k  ipi,v>'  foc  2‘ 


=E*(«PlJ  s“(*‘Ib)  ' 

Ea>  >A.B  *E,(ePl<J  “K*1 '»)  '  1  ^ 


(V/m) 


(V/m) 


(V/m) 


Electric  Field  Integral  Equation  Vector  Components  Px,yi 

Iaref 

Eax(rp,+I  )V 

0 

I2" 

Px(02.42).=  Ax(02.*2,*l  )d*l 

Jo 

Iaref 

Eay(rp,*l)e> 

0 


Pv(02,*2)-=  Ay(02,*2.*l)d*l 

Jo 


(V/m) 


(V/m) 


(V/m) 


(V/m) 


Axl  Eaxl  eJ  ^ ""W)  .p,  *pi  (V/m) 

A.B  A.  B 


PX,A.B:*  ^A.B*1' 


(V/m) 


■  kfp|,,n(®2»)  (cotfllln)  co»(*22n)  +  nn($llp)  «n(*22n))  ... 

AylAD:«  Eayl  eJ  s  *  K  v  s  ^  K  *  '^"rpldrpl  (V/m) 

A,B  A.  a 


•2 

:« 

V,B 

Jo 


Ay  i  *  i 


(V/m) 
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(V/m) 


Secondary  Field  Electric  Field  Components  Eft, _ E^j_ 

I  k  e'J  k’*'2 

E02(82,42)  =i - (Px(02.42)  cos(42)+  Pv(02.42)  sm(42)) 

2  *rfl2 


i  k  e‘J  k  rfr2 

£42(02,42)  zL±l - (Px(02.42)  -sin(42)+  Py(02.42)  cos(42))  cos(02) 

2  n  rff2 


E03 


-J  ke 


-j  k  rlT2 


AB  2nrfl2 


(Px|A.Bcos(422B)  +  Py1A,BS.n(422B)) 


E43 


A.B 


.j  k e‘J  krfn 
2  u  rf!2 


Secondary  Field  Radiation  Intensity  U: 


U2(02.42) 


:=—[(|E02(02. 42)|  )2+  (|E42(02. 42)|  )2] 

2% 


U3a.b:~[(IE93..bI)!*(Ie*5*.b|)!] 


2% 


Secondary  Field  Radiated  Power  Praf)' 


Prad2  =  1.84163  10~" 


Secondary  Field  Directivity  Dn: 

Do22:=4"^1(U3) 

Prad2 


Do22  =  4.59318  105 


(V/m) 


(V/m) 


(V/m) 


(W/solid  ang) 


(W/solid  ang) 


(W) 


(W) 


(dimensionless ) 


(dimensionless) 
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Spillover  Efficiency  cc  : 


*0o 

U(0p )  sin(0p)  d0p 
1CT6 


2 

U(0p)  sm(0p)  d0p 

1(T6 


(dimensionless ) 


cs=  1  00042 


(dimensionless ) 


Taper  Efficiency  : 


£t 


d0p 


U(9p)  —  sin(ep)d6p 

Prad 

0 


(dimensionless ) 


tt  =0  45135 


(dimensionless ) 


Phase  Efficiency 


.4*6  10  iref 


cp  :«  e 


(dimensionless ) 


q>  =0.97222 


(dimensionless ) 


Aperture  Efficiency 

ex  =98  (dimensionless)  eohmic  •=  98  (dimensionless) 


sap  •- es  d  cp  oi  tb  csp  cohniic  (dimensionless)  sap  =0  38452  (dimensionless) 
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Effective  Isotropic  Radiated  Power  EIRP: 


EIRP22  •=  Prad2  Do22 

EERP22  =  8.45893- 10-6 

Gain  G: 

G22'»tapDo22  (dimensionless) 

G22  =  1  76615*  105  (dB> 

Effective  Aperture  A„: 

Ae22  '*  cap  n  aref3 

Ae22  =0  1237 

Radiation  Resistance  Rr: 

Rr22:.2™l 

(llol)2 

Rr22  =  3  68326  10~" 

Effective  Height  h0: 

he22  :=2  Irt22 

J  % 

he22  =2  19868  10-7 


IW) 


(W) 


G22dB  -  IO!og(G22)  (dimensionless) 

G22dB  = 52 47028  (dB) 


(m2 ) 


(m2) 


(Q) 


(Q) 


(m) 


(m) 
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Polarization  Loss  Factor  PLF: 


0p2  =  atan 


( radians ) 


l radians ) 


8p2  *0  95532 


( radians ) 


♦p2  =  0  7854 


( radians ) 


E84  :z  E02(6p2,^p2)  (V/m) 


Ef*  :=E*2(ep2,*p2)  (V/m) 


Ev2  :=E04  cos(0p2)  cos(4p2)-  E+4  sin(4p2) 


(V/m) 


Ey 2  :=  E04  cos(  0p2 )  sin(  $p2 )  -  E+4  cos(  4p2 ) 


(V/m) 


Ez2  :=-E64  sin(9p2) 


(V/m) 


fex2\ 

oa2  . .  1  - |E v2  1 

J(  |Ex2|  )24  ( |Ey2|  )24(|Ez2|)2  \ez2 / 


(dimensionless ) 


oa2  = 


-0.40672  +  0.0353 1 j 
-0.40672  +  0.0353  lj 
,0.81344-  0  07063j 


(dimensionless ) 


PLF2  := 


(dimensionless ) 


PLF2  -  0  16667 


(dimensionless ) 
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s 


H- PLANE  RADIATION  PATTERN: 


Y2(62) 


X2( 62 )  :■  !e42(02.O)| 


Y2(02 )  =  |E*2<  62,0)1  sin|o2-*-^J 


1?  10 


118  10~* 


1  06*10 


94*10 


8  2  10 


7*10_7 


5  8*10 


4.6  10 


3.4  10 


2.2*10 


1  10 


r  7 
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